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1.0  INTRODUCTION 


This  report  describes  the  third  year  of  work  on  Contract  Number 
N0001 4-80-001 7.  The  work  has  the  following  general  goals: 

Demonstration  of  a  superconducting  300  megasample  per  second 
A/D  converter  within  the  first  three  years. 

Demonstration  of  a  superconducting  A/D  converter  operating 
well  into  the  gigasample  per  second  region  within  five  years. 

During  the  second  year  the  goal  for  the  third  year  was  completed. 

During  the  third  year  the  goal  for  the  fifth  year  was  essentially  accomplished 
and  described  in  a  publication  in  Electron  Device  Letters. 

Future  efforts  under  this  contract  will  involve  a  continued  effort  to 

*  raise  the  sample  rate  and  to  accurately  characterize  the  device  at  GHz 
sampling  rates.  In  addition  work  will  begin  on  efforts  to  incorporate  the 
analog-to-digital  converter  into  a  transient  event  recorder.  This  entire 
device  will  be  fabricated  in  superconducting  integrated  circuit  technology. 

2.0  A/D  CONVERTER 
f 

2.1  2  GHz  A/D  Sample  Rate 

During  the  contract  period  the  performance  of  the  A/D  converter  has  been 
improved  from  4-bit  accuracy  at  a  200  MHz  sampling  rate  (see  Appendix  A1 )  to 
6-bit  accuracy  at  a  2  GHz  sampling  rate.  This  improvement  has  been  made 

*  possible  by  a  new  design  which  allows  a  substantial  size  reduction  and  a  new 
sample  holder  with  reduced  crosstalk.  The  new  design  is  based  on  2-junction 
interferometers  which  are  smaller  and  have  much  more  periodic  threshold 
curves  than  the  3-junction  devices  used  previously.  The  2- junction  inter¬ 
ferometers  also  improve  the  operating  margins  and  are  less  sensitive  to 
trapped  flux. 

The  sample  holder  has  been  modified  to  increase  the  number  of  contacts 
from  24  to  32  thus  allowing  a  separate  ground  wire  for  each  of  the  eight 
high  speed  contacts.  The  resulting  reduction  in  crosstalk  makes  possible 
operation  of  the  converter  above  1  GHz.  Even  with  this  Improvement,  sample 
holder  crosstalk  Is  the  limiting  factor  at  the  highest  operating  speeds.  The 
1  new  design  and  experimental  results  are  described  In  detail  in  a  publication 

reproduced  In  Appendix  A2:  "A  Superconducting  6-B1 t  Analog- to-Digital  Con¬ 
verter  with  Operation  to  2  x  10s  Samples/Second." 

2.2  Direct  Injection  A/D  Converter 

«  An  alternative  to  the  magnetically-coupled  signal  line  of  previous  designs 

Is  the  direct  Injection  of  the  signal  current  into  the  interferometer 
comparator.  An  A/D  converter  based  on  this  idea  has  been  designed  and  suc¬ 
cessfully  tested.  Its  schematic  is  shown  in  Fig.  1A  and  a  circuit  photo  is 
shown  in  Fig.  IB.  A  resistor  ladder  is  used  to  divide  the  signal  in  a  binary 
sequence.  Each  signal  component  is  then  Injected  through  an  interferometer 
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inductance  to  ground.  The  use  of  identical  interferometers  and  the 
symmetry  of  the  resistor  ladder  assures  the  achievement  of  the  desired 
binary  coupling  ratios— independent  of  processing  parameters.  The  absence  of 
insulated  control  lines  and  a  reduced  number  of  crossovers  makes  the  circuit 
simpler  to  fabricate  and  more  reliable.  It  is  also  expected  that  crosstalk 
between  interferometers  along  the  signal  line  will  be  greatly  reduced.  The 
principal  disadvantage  of  this  design  is  that  the  signal  must  be  divided  and 
applied  to  the  comparators  in  parallel  rather  than  the  series  technique  used 
previously.  This  results  in  reduced  sensitivity.  Two  changes  in  the  interfe¬ 
rometers  have  been  made  to  help  restore  the  sensitivity:  (1)  The  ground  point 
has  been  moved  producing  an  asymmetric  interferometer  in  which  the  injected 
current  flows  through  about  90%  (rather  than  50%)  of  the  total  inductance,  and 
(2)  the  interferometer  inductance  has  been  increased  from  8  to  12  pH.  This 
second  change  reduces  the  depth  of  modulation  of  the  threshold  curve  with  a 
slight  reduction  in  operating  margin.  The  present  full-scale  sensitivity  for 
the  direct  injection  design  is  5.7  mA  compared  to  3.5  mA  for  the  previous 
magnetic  coupling  design. 

2.3  Sample  Holder  Improvements 

Since  crosstalk  in  the  sample  holder  remains  a  limiting  problem  for 
high  speed  measurements,  we  are  making  continuing  efforts  to  design  a  better 
sample  holder.  The  difficulty  of  the  crosstalk  problem  is  illustrated  by 
the  fact  that  the  present  A/0  converter  operates  at  clock  rates  approximately 
three  times  the  rate  being  used  in  prototype  Josephson  computers. 

One  of  the  causes  of  crosstalk  is  that  metal  parts  in  the  sample  holder 
cannot  be  solidly  grounded  to  the  chip  ground  plane.  Stray  capacitance 
between  the  signal  lines  and  these  metal  parts  provides  a  crosstalk  path. 

The  most  recent  sample  holder  approaches  this  problem  by  using  plastic  parts 
to  reduce  stray  capacitance  and  by  using  ground-signal -ground  coplanar  strip 
lines  to  bring  signals  from  coaxial  cables  to  the  chip.  This  results  in  a 
38  pin  chip  with  14  slow  lines,  and  8  fast  lines  with  two  ground  contacts 
each.  A  photo  of  the  new  sample  holder  is  shown  in  Fig.  2.  Evaluation  of 
this  device  is  currently  In  progress. 

2.4  A/D  Converter  Test  Apparatus 

Performance  measurements  on  the  A/D  converter  have  in  the  past  been 
made  with  two  different  experimental  arrangements.  The  first  operates  the 
converter  at  modest  speed  with  a  slowly  varying  signal.  This  mode  of  operation 
is  appropriate  for  making  the  required  bias  adjustments  and  for  observing  the 
staircase  function.  For  high  speed  measurements  the  converter  must  be  connected 
to  high  frequency  clock  and  signal  generators  and  a  sampling  oscilloscope  for 
observing  the  output  data.  During  this  reconnection  process  it  is  likely  that 
the  bias  adjustments  will  drift  thus  rendering  the  output  data  inaccurate. 
Clearly,  a  method  is  required  for  quickly  making  and  maintaining  correct  adjust¬ 
ment  during  high  speed  operation.  An  algorithm  for  automatic  adjustment  is 
described  in  Appendix  A2.  The  implementation  of  this  algorithm  in  an  all-purpose 
test  fixture  is  shown  in  the  block  diagram  of  Fig.  3  and  photo  of  Fig.  4.  The 
test  fixture  will  (hopefully)  make  and  maintain  the  bias  adjustments,  display 
the  bit  curves  and  staircase  function  and  allow  high  speed  operation  at  any 
time. 
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This  is  accomplished  by  tapping  the  data  transmission  lines  with 
small  IK  resistors  which  allow  the  measurement  and  injection  of  low 
frequency  signals  without  disturbing  the  high  frequency  signals.  A  nest 
of  12  interacting  servo  loops  operates  on  these  low  frequency  signals  to 
make  the  required  adjustments.  A  multiplexer  and  D/A  converter  provide 
scope  displays  which  completely  characterize  the  low-frequency  behavior 
of  the  A/D  converter. 

2.5  Circuit  Simulations 

In  the  past  most  of  our  circuit  simulations  have  been  performed  using 
the  SUPER*SCEPTRE  network  analysis  program.  Because  of  the  size  of  this 
program  it  requires  a  large  computer,  response  is  slow,  and  it  is  costly 
to  run.  Many  cases  can  be  solved  using  a  much  simpler  approach  which  operates 
as  follows:  over  small  time  increments,  capacitors  are  assumed  to  be  con¬ 
stant  voltage  sources,  and  inductors  and  junctions  are  assumed  to  be  constant 
current  sources.  The  resulting  simplified  circuit  equations  are  then  solved 
for  the  currents  into  capacitors  and  the  voltages  across  inductors.  A  simple 
numerical  integration  then  updates  the  capacitor  voltages,  and  junction  and 
inductor  currents.  Repeating  this  process  for  each  time  increment  produces 
the  complete  time  dependent  behavior  of  the  circuit. 

A  general  version  of  the  program  (which  we  call  micro*SCEPTRE)  for  a 
circuit  with  an  arbitrary  number  of  components  requires  only  about  30  lines 
of  BASIC  or  FORTRAN  code.  This  program  runs  easily  on  any  desktop  computer 
and  has  duplicated  SUPER*SCEPTRE  results  even  on  moderately  complex  circuits 
(15  elements).  When  programmed  in  SASIC,  running  times  are  in  minutes  to 
tens  of  minutes.  In  either  case,  results  are  usually  obtained  more  quickly 
than  when  using  SUPER*SCEPTRE.  An  example  of  an  application  for  this  program 
is  described  in  the  following  section  on  aperture  time. 

One  shortcoming  of  this  approach  is  that  the  user  must  learn  a  procedure 
for  generating  equations  to  be  used  in  programming  the  method.  It  thus  lacks 
the  topological  input  description  of  the  circuit  found  in  SUPER*SCEPTRE.  It 
appears  that  a  combination  of  approaches  will  be  required  to  rapidly  obtain 
circuit  simulations.  These  involve  changing  the  method  of  simulation  from 
SUPER*$CEPTRE  to  another  more  efficient  method  and  acquisition  of  a  highly 
responsive  time-sharing  computer.  Ideally,  an  engineer  should  be  able  to  sit 
down  at  a  computer  terminal,  enter  the  elements  of  the  circuit,  and  see  the 
simulated  response  of  the  circuit  on  his  screen  in  seconds  or  minutes.  With 
such  highly  responsive  simulation,  he  should  be  able  to  adjust  component 
values,  or  other  circuit  details,  interactively,  much  in  the  way  one  develops 
a  breadboard  circuit  with  slower  components. 

3.0  APERTURE  TIME 

The  sampling  theorem  demonstrates  that  a  signal  can  be  reconstructed 
from  its  values  at  a  number  of  sampling  instants  if  the  length  of  time  over 
which  the  sample  is  taken  is  infinitely  short.  In  practice  the  length  of 
time  must  be  short  enough  that  one  makes  an  error  in  the  value  of  the  signal 
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of  less  than  one  least  significant  bit.  This  length  of  time  is  called  the 
"aperture  time".  For  a  bipolar  signal  the  maximum  aperture  time  is  given 
by  the  following: 


where  f  is  the  sampling  frequency  and  N  is  the  number  of  bits.  Thus  for 
a  6-bitsconverter,  with  a  2  GHz  sampling  rate,  the  aperture  time  must  be 
less  than  5  ps  and  for  a  4-bit  converter,  less  than  20  ps.  In  this 
section  we  discuss  several  approaches  to  limiting  the  aperture  time  to 
such  short  intervals. 

The  reader  should  note  that  the  sampling  theorem  also  allows  one  to 
measure  the  integral  of  the  signal  over  a  longer  interval  whose  end  points 
are  accurately  known.  Although  this  approach  may  prove  useful  in  the 
future,  it  has  not  yet  been  considered  for  implementation  in  superconducting 
technology. 

3.1  Pulsed  Bias  Supply 

The  required  short  aperture  time  can  be  provided  by  an  additional 
circuit,  or  by  building  it  intrinsically  into  the  converter. 

An  approach  to  the  latter  is  the  use  of  self-gating  "AND"  comparators 
as  proposed  by  Van  Duzer's  group  at  Berkeley.  These  comparators  are 
sensitive  to  the  signal  only  during  the  rising  edge  of  the  power  supply 
waveform,  thus  yielding  an  aperture  time  equal  to  the  supply  rise  time. 

This  approach  doubles  the  complexity  of  the  circuit  and  may  have  rather 
small  margins  for  error-free  operation. 

Another  approach  suggested  by  Dennis  Herrell  of  IBM  is  to  drive  the 
interferometer  comparators  with  a  power  supply  pulse  which  makes  a  single 
very  short  excursion  to  the  normal  level  for  comparator  operation  (about 
0.75  Ic)  and  then  falls  back  to  a  point  below  the  "floor"  of  the  interfe¬ 
rometer  gain  curve  for  the  remainder  of  the  cycle  (see  Fig.  5).  In  this  way 
the  signal  is  sampled  only  for  a  very  short  time  at  the  beginning  of  the 
cycle.  During  the  rest  of  the  cycle  the  comparators  are  insensitive  to 
signal  changes  but  retain  the  digitized  data  for  transfer  to  succeeding 
logic  stages. 

A  method  to  generate  the  required  pulse  is  to  differentiate  the  switching 
transition  of  a  single  junction  by  using  a  set  of  small  capacitors  which 
couple  the  junction  to  each  of  the  interferometers.  The  circuit  (with  only 
one  comparator  interferometer  shown)  is  illustrated  in  Fig.  6. 

It  operates  as  follows:  the  sine  wave  power  supply  applies  bias  to 
the  comparator  interferometer  (J2,  J3)  and  to  the  pulse  generating  junction 
Jl .  Near  the. peak  of  the  sine  wave,  Jl  switches,  thus  delivering  a  pulse  of 
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current  through  capacitor  C  to  the  comparator.  It  will  switch  or  not  depending 
on  the  signal  level.  For  the  remainder  of  the  cycle  the  comparator  remains  in 
the  same  state  regardless  of  the  signal  level.  R3  damps  the  principal  resonance 
of  the  circuit.  If  the  capacitance  of  Jl  is  large  compared  to  C,  the  necessary 
isolation  between  comparators  can  be  achieved. 

This  circuit  was  simulated  using  our  micro*SCEPTRE  analysis  program.  Figure 
7  is  a  plot  of  the  current  I  to  the  interferometer.  The  effective  aperture  time 
is  the  width  of  the  spike  which  occurs  near  the  peak  of  the  sine  wave. 

This  circuitry  has  been  added  to  the  latest  version  of  our  direct  injection 
A/D  converter  and  is  shown  at  the  top  of  Figs.  1A  and  IB.  The  timing  of  the 
sampling  pulse  can  be  adjusted  by  applying  a  dc  bias  to  the  pulse  generating 
junction  via  pin  5.  The  amplitude  of  the  sampling  pulse  can  be  measured  by 
observing  the  apparent  decrease  in  the  critical  currents  of  the  interferometers 
when  the  pulse  is  present.  These  measurements  show  the  pulse  amplitude  to  be 
about  30%  less  than  the  design  value.  When  this  defect  is  corrected  in  the 
next  mask  set  it  is  expected  that  operation  of  the  sample-and-hold  can  be 
demonstrated. 

3.2  Persistent  Current  Sample-and-Hold 


Another  approach  to  achieving  a  short  aperture  time  is  to  add  a  separate 
sample  and  hold  circuit  to  the  input  of  the  A/D  converter.  This  can  be  done  by 
trapping  a  persistent  current  proportional  to  the  signal  in  a  superconducting 
loop  as  illustrated  in  Fig.  8.  In  this  circuit  a  control  current  reduces  the 
critical  current  in  the  junction  J  to  a  low  value.  As  a  result  the  junction 
switches  into  the  voltage  state  and  the  current  is  transferred  into  tne  super¬ 
conducting  loop  having  inductance  L.  As  long  as  the  signal  current  remains 
above  the  reduced  critical  current  of  the  junction,  the  signal  current  will 
flow  predominantly  in  L.  When  a  sample  is  to  be  taken  the  control  current 
above  J  is  reduced  allowing  the  critical  current  to  rise.  The  signal  current  is 
such  that  it  always  remains  below  the  unreduced  critical  current  of  J.  When 
the  critical  current  of  J  is  increased,  the  signal  current  is  trapped  in  L, 
and  any  further  changes  in  signal  current  flow  through  the  junction,  rather 
than  L,  because  L  has  a  much  larger  inductance  than  that  of  J. 

Simulations  of  this  device  show  that  the  current  in  L  is  indeed  trapped 
when  the  critical  current  rises,  and  further  changes  in  the  signal  flow  through 
J.  The  loop  L  is  used  as  the  control  current  for  the  interferometers  in  the 
A/D  converter.  Thus  with  proper  timing,  the  A/D  always  samples  a  static 
signal.  The  aperture  time  is  just  the  time  required  to  switch  J  into  the 
superconducting  state. 

The  simulations  were  done  using  a  single  junction  with  controllable 
critical  current.  The  circuit  looks  a  good  deal  less  ideal  when  the  single 
junction  is  replaced  with  a  2-junction  interferometer.  Several  complications 
should  now  be  discussed. 


Figure  7.  The  results  of  a  circuit  simulation  of  Figure  5.  The  A/D 
aperture  time  is  the  width  of  the  spike  which  occurs  at  the 
peak  of  the  sine  wave  supply  current. 
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Figure  8.  A  schematic  diagram  for  a  persistent  current  sample  and  hold 
circuit. 
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1.  The  magnetic  flux  corresponding  to  the  persistent  current  in 
inductance  L  is  quantized.  Thus  L  must  be  large  enough  that  a  number  of 
flux  quanta  greater  than  the  number  of  resolvable  levels  of  the  converter 
can  be  trapped.  Thus 


L  > 


Wx 


where  i  is  the  current  corresponding  to  the  maximum  digital  output  of  the 
converterxand  <j>  =  2.07  x  10"  3  weber.  Since  1  is  3.5  mA  for  our  most 

sensitive  converter,  L  must  be  at  least  40  pH.  inis  is  not  difficult 
and  indeed  L  may  be  several  hundred  picohenries  if  the  crossing  inductances 
of  the  interferometers  are  included  as  they  must  be. 

2.  Unfortunately,  it  requires  a  large  voltage  to  change  the  current  in 
a  large  inductance  quickly.  In  our  case  the  voltage  is  limited  to  the  gap 
voltage  of  the  junction,  about  2.9  mV.  The  response  time  of  the  loop  is  thus 
about  50  ps.  This  is  a  reasonable  speed,  since  the  sampling  time  of  our 
present  converters  is  500  ps.  Furthermore,  if  the  sample-and-hold  is  operated 
in  the  tracking  mode  indicated  above,  this  long  response  time  is  only  evident 
after  a  sample  when  the  signal  current  is  again  steered  through  the  loop  L. 

3.  Most  difficult  of  all,  however,  is  the  necessity  to  reduce  the 
critical  current  of  the  interferometer  below  the  current  corresponding  to  a 
least  significant  bit  of  signal  current.  This  may  impose  a  serious  limitation 
on  the  dynamic  range  of  this  sample  and  hold  circuit. 


4.0  OVERALL  SYSTEM  CONCEPTS 

Work  on  this  contract  has  progressed  to  the  point  that  thought  is  being 
given  to  the  use  of  the  A/D  converter  in  a  complete  system.  Not  only  is  a 
complete  system  necessary  for  the  practical  use  of  the  converter,  it  also 
appears  to  be  imperative  in  order  to  test  the  device.  Data  is  produced  so 
fast  by  the  device  that  no  conventional  method  of  storing  it  exists.  Thus 
the  initial  system  will  Include  fast  superconducting  memory  to  store  the  data 
and  later  allow  it  to  be  read  out  more  slowly  into  a  conventional  computer 
for  analysis.  The  system  will  in  fact  be  a  transient  event  recorder. 

One  possibility  for  the  system  Is  Illustrated  in  Fig.  9.  One  bit 
of  the  A/D  converter  is  shown  being  driven  from  a  gated  sine  wave  supply 
together  with  a  dc  offset.  The  same  supply  is  split  into  three  phases 
by  delay  lines  and  operates  a  shift  register  consisting  of  a  long  string 
of  AND  gates.  Thus  a  burst  of  M  cycles  from  the  sine  wave  supply  converts 
and  stores  M  samples  of  the  Input  signal.  The  dc  offset  in  the  supply 
causes  the  data  to  be  retained  after  shifting  stops.  A  low  frequency 
supply  can  then  be  used  to  shift  the  data  out  for  further  processing. 
Preliminary  design  is  underway  for  such  a  shift  register  memory. 


Figure  9.  A  schematic  diagram  for  an  A/D  converter  (only  1  bit  shown) 
with  shift  register  memory. 


Figure  10.  Typical  I-V  curve  for  a  junction  with  a  Pb-Bi  upper  electrode. 
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5.0  FABRICATION  CHANGES 

This  section  describes  several  changes  which  have  been  made  in  the 
fabrication  process  during  the  past  year.  In  addition  one  anticipated 
improvement  is  discussed. 

5.1  New  Mask  AT igner 

A  new  Cobilt  2020H  mask  aligner  is  now  in  routine  use.  This  contact 
printer  may  be  the  most  sophisticated  available  and  enables  alignment  to 
within  about  one  micron.  In  addition  to  its  improved  optics  and  alignment 
capability,  it  also  offers  greatly  improved  exposure  control,  uniformity, 
and  reliability. 

5.2  Pb-Bi  Upper  Junction  Electrodes 

A  Lead-Bismuth  alloy  is  now  used  for  junction  upper  electrodes. 

Junction  electrical  quality  is  extremely  high  and  quite  reproducible. 

An  I-V  curve  from  such  a  junction  is  shown  in  Fig.  10.  A  significant 
difficulty  remains,  however,  with  this  alloy.  Films  have  very  large 
grain  size  (about  one  micron)  and  edge  coverage  is  unreliable.  Evidently 
this  problem  can  be  alleviated  by  cooling  the  substrate  during  deposition. 
See  the  section  below  on  temperature  control. 

5.3  New  Nb  Deposition  Capability 

Early  in  this  program  niobium  ground  planes  were  deposited  in  an  ultra- 
high  vacuum  system  with  an  electron  beam  evaporator.  A  catastrophic 
failure  of  the  evaporator  later  caused  an  outside  company  to  be  used  for 
Nb  deposition. 

In  the  last  two  months,  however,  a  new  Nb  deposition  system  has 
been  installed.  It  consists  of  a  small  dc  magnetron  sputter  gun  installed 
in  a  conventional  diffusion  pumped,  nitrogen  trapped  vacuum  system.  The 
source  is  only  large  enough  to  coat  one  wafer  at  a  time,  so  a  rotating 
sample  holder  was  constructed  to  enable  coating  multiple  wafers  in  one 
pumpdown.  The  system  produces  niobium  films  having  essentially  bulk 
transition  temperature  with  thicknesses  of  a  few  thousand  angstroms. 

We  have  done  intial  experiments  in  which  the  Nb  films  are  patterned 
using  our  usual  liftoff  technique.  Edge  quality  has  been  satisfactory. 

The  new  deposition  system  has  enabled  us  to  begin  making  Nb-Pb  alloy 
junctions  which  have  future  potential  for  a  high  reliability  technology. 

5.4  Future  Need  for  Wide-Range  Temperature  Control  of  Substrates 

At  present  the  substrate  temperature  control  in  our  deposition- 
oxidation  system  covers  a  range  of  about  10  to  50°  Celsius.  Two  processes 
which  we  may  wish  to  use  in  the  future  require  the  range  to  be  extended. 
First,  as  mentioned  above,  Pb-Bi  films  exhibit  very  large  grain  size 
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which  can  be  reduced  by  deposition  onto  0°  Celsius  substrates.  Second,  a 
new  technique  is  about  to  appear  In  the  literature  for  fabricating  very 
fine  grained  lower  electrodes.  Again  the  grain  size  is  reduced  using  a 
low  temperature  substrate,  in  this  case  77  Kelvin. 

Our  present  vacuum  system  is  designed  such  that  achieving  this 
range  of  temperature  will  involve  considerable  work.  This  work  will  be 
done  during  the  next  year. 


6.0  RELATED  WORK 

Additional  work  done  by  the  Cryoelectronic  Group  of  NBS  adds  to  the 
total  expertise  which  can  be  brought  to  bear  on  the  A/D  converter  problem. 

In  this  section  we  discuss  four  such  publications,  some  of  which  were 
supported  by  ONR  under  other  contracts.  The  publications  are  included  in 
this  report  as  appendices. 

6.1  Turn-On  Delay  in  Junctions  and  Interferometers 

R.  L.  Peterson  has  simplified  previous  results  on  turn-on  delay  and 
considerably  extended  their  range  of  validity.  His  results  are  directly 
applicable  to  interferometers  used  in  the  A/D  converter.  They  will  play  an 
important  part  in  eventual  estimates  of  the  limiting  speed  of  operation  of 
the  converter.  The  results  for  single  junctions  will  appear  in  the  Proceedings 
of  IC  SQUID,  Berlin,  May  1,  1980,  (Appendix  A3). 

6.2  Review  of  Picosecond  Response  of  Josephson  Devices 

A  number  of  members  of  the  Cryoelectronic  Group  have  reviewed  the 
expected  limiting  response  of  Josephson  devices  in  a  paper  to  appear  in 
IEEE  Transactions  on  Devices,  October,  1980,  (Appendix  A4). 

6.3  Microwave  Mixers  and  Detectors 

Collaborative  work  with  Professor  P.  L.  Richards  and  T.  M.  Shen  has  led 
to  the  experimental  observation  of  conversion  gain  in  SIS  mixers.  These 
mixers,  produced  at  NBS,  exhibit  the  best  signal-to-noise  ratios  of  any 
mixers  ever  tested  at  36  GHz,  (Appendix  A5). 

6.4  Analog  Sampling 

An  analog  sampling  device  has  been  constructed  which  has  been  used  to 
observe  the  shortest  transition  time  ever  seen  In  an  electric  circuit,  9  ps. 
Further  Improvements  in  this  device  should  lead  to  even  better  sampling  times. 
An  effort  Is  now  underway  to  couple  the  device  to  signals  generated  at  room 
temperature,  (Appendix  A6). 
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7.0  INFORMATION  EXCHANGE 


Leadership 

R.  E.  Harris  served  as  vice  chairman  and  local  arrangements  chairman  of 
the  Workshop  on  Josephson  Dfgital  Devices,  Circuits,  and  Systems,  Copper 
Mountain,  Colorado. 

C.  A.  Hamilton  served  on  the  organizing  committee  for  the  Workshop  on  High 
Speed  A/D  Conversion  -  1980,  Portland,  Oregon. 


Invited  Talks 


1.  "Superconducting  Electronics",  seminar  at  Tektronix,  Inc.,  Beaverton, 

Oregon,  C.  A.  Hamilton,  June  4,  1979. 

2.  "A  High  Speed  Superconducting  A/D  Converter",  37th  Annual  Device 
Research  Conference,  Boulder,  Colorado,  C.  A.  Hamilton,  June  25,  1979. 

3.  "Introduction  to  Josephson  Devices",  Workshop  on  Josephson  Digital  Devices, 
Circuits,  and  Systems",  Copper  Mountain,  Colorado,  R.  E.  Harris,  August  21,  1979. 

4.  "Superconducting  Striplines",  Workshop  on  Josephson  Digital  Devices, 

Circuits,  and  Systems",  Copper  Mountain,  Colorado,  R.  L.  Kautz,  August  21,  1979. 

5.  "An  A/D  Converter  Based  on  Quantum  Interference",  Workshop  on  Josephson 
Digital  Devices,  Circuits,  and  Systems",  Copper  Mountain,  Colorado,  C.  A. 
Hamilton,  August  22,  1979. 

6.  "A/D  Conversion  Using  Superconducting  Quantum  Interference",  Workshop  on 
High  Speed  A/D  Conversion  -  1980,  Portland,  Oregon,  C.  A.  Hamilton,  February 
12,  1980. 

7.  Panel  on  "Competing  Technologies  for  Gb  Logic",  International  Solid-State 
Circuits  Conference,  San  Francisco,  California,  R.  E.  Harris,  with  others, 
February  13,  1980. 
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Multiple-quantum  interference  superconducting  analog-to-digitai 
converter** 

Richard  E.  Harris,  C.  A.  Hamilton,  and  Frances  L.  Lloyd 

Electromagnetic  Technology  Division.  National  Bureau  of  Standards,  Boulder,  Colorado  80303 
(Received  16  March  1979;  accepted  for  publication  1 1  May  1979) 

Multiple-quantum  interference  in  a  superconducting  interferometer  is  used  for  analog-to-digital 
conversion.  The  simple  fully  parallel  four-bit  converter  which  is  described  is  the  first  known  use 
of  this  effect  in  a  digital  circuit.  Sampling  rates  of  2X  10*  per  second  were  achieved,  and  much 
higher  rates  appear  possible. 

PACS  numbers:  06.60.Jn,  74.50.  +  r,  06.70.Dn,  06.70.Ep 


This  letter  describes  the  first  known  use  in  a  digital  cir¬ 
cuit  of  the  multiple  lobes  of  the  threshold  curve  of  a  super¬ 
conducting  interferometer.  Although  digital  applications  of 
such  interferometers  have  been  known  for  several  years,12 
only  the  first  lobe  has  been  used  previously.  An  analog-to- 
digital  (A/D)  converter  using  multiple  lobes  has  been  pro¬ 
posed.3  4  We  demonstrate  a  realization  of  this  concept  in  a 
four-bit  A/D  converter  operating  at  a  rate  of  200  megasam¬ 
ples  per  second.  An  increase  to  rates  well  above  one  gigasam- 
ple  per  second  may  be  possible. 

This  experiment  makes  use  of  four  current  compara¬ 
tors,  each  consisting  of  a  Josephson  three-junction  interfer¬ 
ometer  which  has  the  junctions  connected  in  parallel  by  su¬ 
perconducting  loops.’  A  photograph  of  one  of  the 
interferometers  is  shown  in  Fig.  1  (a).  The  center  junction 
has  twice  the  critical  current  of  each  of  the  symmetrically 
arranged  outer  junctions.  Such  an  interferometer  has  a  cur¬ 
rent-voltage  characteristic  like  that  of  a  single  junction.  Its 
maximum  0-V  current  Im  (threshold  current)  is  modulated 
by  an  applied  magnetic  field  in  a  way  determined  by  the 
junction  critical  currents  and  the  loop  inductances.  The 
junctions  are  sufficiently  small  (about  4^m  square)  that  the 
magnetic  fields  used  have  negligible  effect  on  them  directly. 
By  choosing  Li0  s  i<A0  • a  threshold  curve  like  that  shown  in 
Fig.  2(a)  is  obtained.  Here  L  is  the  inductance  of  one  of  the 
superconducting  loops,  i0  is  the  critical  current  of  one  of  the 
outer  junctions,  and  <bn  is  the  flux  quantum  h  /Te. 

The  periodicity  of  the  threshold  curve  is  a  manifestation 
of  quantum  interference  occurring  in  the  interferometer. 
The  net  phase  change  of  the  wave  function  describing  the 
superconducting  loops  in  the  interferometer  increases  in 
proportion  to  the  amount  of  magnetic  flux  applied  through 
the  loops  of  the  interferometer.  A  net  phase  change  of  2rr 
corresponds  to  one  additional  flux  quantum  penetrating  the 
loop.  The  threshold  curve  repeats  whenever  two  additional 
flux  quanta  are  added.  In  a  two-junction  interferometer,  the 
curve  would  repeat  for  each  additional  flux  quantum. 

To  use  this  interferometer  as  a  current  comparator,  a 
magnetic  field  is  applied  equally  through  both  loops.  The 
magnetic  field  results  from  passing  a  current  lc  through  a 
control  line  above  the  interferometer.  Its  effect  is  illustrated 
by  the  horizontal  arrow  in  Fig.  2(a).  The  current  through  the 


interferometer  It  is  then  raised  from  zero  to  a  predetermined 

level  Is  mtl  indicated  by  the  vertical  arrow.  If  It is  above 

the  threshold  curve  the  interferometer  switches  to 
the  voltage  state;  otherwise  it  does  not. 

A  single  interferometer  is  used  to  generate  a  single  bit  of 
the  output  of  an  A/D  converter,  by  choosing  so  that 
the  interferometer  switches,  or  does  not  switch,  for  equal 
intervals  of  control  current.  By  making  the  control  current 
the  analog  signal  to  be  digitized,  the  voltage  output  of  the 
gate  Vt  is  one  bit  of  the  digital  value  of  the  signal  as  illustrat¬ 
ed  in  Fig.  2(b).  To  take  repeated  samples,  lt  is  cyclically 
varied  between  0  and  It  mtx  ■ 

The  next  most  significant  bit  can  be  generated  by  ex¬ 
panding  the  bit  pattern  of  the  least  significant  bit  by  a  factor 
of  2.  The  expansion  is  produced  by  reducing  the  coupling  of  a 
second  interferometer  to  the  control  line  by  a  factor  of  2. 
Similarly,  the  third  and  fourth  bits  are  generated  by  further 
expanding  the  bit  patterns  of  third  and  fourth  interferomet¬ 
ers.  Thus  only  four  interferometers  are  required  for  A/D 
conversion  of  four-bit  accuracy.  Moreover,  since  these  inter¬ 
ferometers  operate  in  parallel,  maximum  speed  is  achieved 
with  all  four  bits  being  produced  essentially  simultaneously. 
In  general,  only  N  interferometers  are  required  to  digitize  to 
AT-bit  accuracy,  with  little  decrease  in  speed  needed  for  addi¬ 
tional  bits. 
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•’Work  supported  in  put  by  the  Office  of  Neva)  Research  under  contract  FIG.  1  (a)  Photograph  of  one  of  the  four  interferometers  in  the  A/D  con- 
number  NOOOI4-77-F-48.  verter.  (b)  Photograph  of  a  complete  four-bit  converter. 
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FIG.  2.  (a)  Threshold  curve  showing  maximum  zero  voltage  current  Im  as 
a  funcr  -n  of  control  current  Ic.  Ideal  maximum  curve  current  is 
shown  by  dashed  line,  (b)  Output  voltage  for  ideal  as  function  of 
control  current,  (c)  Actual  bit  pattern  for  the  four-bit  A/D  converter,  (d) 
Value  of  digital  output  of  converter  versus  analog  signal  (traced). 


Figure  1  (b)  shows  a  complete  four-bit  converter.  The 
observed  bit  patterns  for  this  converter  are  shown  in  Fig. 
2(c).  This  pattern  is  a  Gray  code  representation  of  the  analog 
signal.  Gray  code  is  frequently  used  in  A/D  conversion  since 
only  one  bit  changes  as  the  signal  crosses  from  one  digital 
level  to  the  next.  Another  significant  feature  of  this  A/D 
converter  can  be  realized  by  adding  a  second  control  line  to 
each  interferometer.  By  applying  a  dc  bias  through  one  of 
these  control  lines,  the  bit  pattern  from  that  gate  can  be  shift¬ 
ed  to  values  of  higher  or  lower  current.  Using  this  feature,  a 
variety  of  unique  codes  can  be  produced. 

Superconducting  A/D  converters  of  the  type  described 
have  been  fabricated  with  two  converters  on  each  6.4-mm- 
square  chip.  The  interferometers  are  much  like  those  de¬ 
scribed  in  Ref.  5,  although  they  are  scaled  up  in  linear  size  by 
a  factor  of  4  to  make  fabrication  simpler.  The  circuits  con¬ 
tain  resistors  for  damping  the  interferometer  resonances  and 
for  terminating  transmission  lines.  The  impedances  of  the 
striplines  on  the  chip  are  chosen  for  proper  matching  to  the 
interferometer  impedances.  The  coupling  of  the  signal  to  the 
interferometers  is  varied  by  adjusting  the  length  of  the  con¬ 
trol  line  above  the  interferometers  in  the  ratios  1:2:4:8.  As 
can  be  seen  in  Fig.  1,  the  control  line  is  diverted  around  the 
portion  of  the  wide  counterelectrode  passing  over  the  inter¬ 
ferometer.  In  this  region  the  counterelectrode  provides 
shielding  which  is  not  readily  calculable.  Furthermore,  this 
control-line  arrangement  causes  the  coupling  to  each  half  of 


each  interferometer  to  be  insensitive  to  slight  misalignment 
during  fabrication. 

The  interferometers  have  been  tested  with  a  gate  cur¬ 
rent  having  a  frequency  (i.e.,  sampling  rate)  of  200  MHz. 
Due  to  equipment  limitations  only  two  of  the  four  bits  were 
tested  at  this  speed.  They  were  found  to  produce  outputs 
which  repeat  synchronously  with  a  nominally  50-MHz  sig¬ 
nal.  In  this  sense,  the  converter  has  been  demonstrated  to 
operate  at  200  megasamples  per  second.  Based  on  the  ob¬ 
served  switching  times,  we  believe  that  the  devices  fabricated 
will  operate  properly  at  rates  up  to  about  700  megasamples 
per  second.  Improving  the  ground  plane  contacts  should  al¬ 
low  the  converter  to  operate  at  much  higher  rates.  Reported 
switching  times®  of  42  ps  for  these  devices  suggest  rates  well 
above  1  gigasample  per  second. 

To  assess  the  accuracy  of  the  full  device,  conventional 
circuitry  was  used  to  convert  from  Gray  code  back  to  an 
analog  signal.  The  resulting  output-vs-input  curve  [Fig. 
2(d)]  is  monotonic  with  a  maximum  error  in  step  width  of  0. 7 
least  significant  bits.  These  measurements  were  performed 
at  about  10  megasamples  per  second  with  the  signal  swept  at 
a  few  hundred  hertz. 

The  coupling  of  t’.e  signal  to  the  interferometers  was 
determined  to  be  in  tht  ratio  1:2:4:8  within  about  1%,  indi¬ 
cating  that  this  design  may  be  suitable  for  converters  of  6  to  8 
bits. 

The  initial  tests  of  this  A/D  converter  design  suffer 
from  variations  in  interferometer  characteristics,  trapped 
flux,  and  the  general  difficulties  of  high-speed  testing.  In 
addition  the  ultimate  accuracy  of  the  device  may  be  limited 
by  dynamic  effects  including  subthreshold  switching  and 
signal  line  transients  caused  by  switching.  Nevertheless,  the 
design  appears  to  offer  great  promise  for  significantly  higher 
speeds  in  the  future. 

It  is  a  pleasure  to  acknowledge  helpful  conversations 
with  T.  Fulton,  L.  Geppert,  D.  Herrell,  and  H.  Zappe.  R. 
Peterson  assisted  with  the  understanding  of  the  threshold 
curves,  D.  McDonald  provided  constant  encouragement, 
and  L.  Mullen  assisted  with  the  fabrication. 
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A  Superconducting  6-Bit  Analog-to-Digital  Converter 
with  Operation  to  2  x  109  Samples /Secondf 

C.A.  HAMILTON,  member,  ieee,  and  FRANCES  L.  LLOYD 


Abstract  —  The  deufn  tad  pcrfornaacc  of  a  6-bit  sapmoaducting 
A/D  coamter  are  described.  The  coavefter  is  baaed  oa  doable  junctioa 
iaterferometera  ased  as  carnal  companion.  The  uoiqac  periodic 
response  of  these  companion  makes  possible  a  fully  parallel  N-bit  con* 
verier  requiring  only  N  companion.  Conversion  rales  up  to  2  x  10’ 
samples  per  second  have  been  demonslnted. 

This  letter  describes  a  superconducting  A/D  converter 
with  6-bit  accuracy  and  operating  speeds  up  to  2  x  109 
samples  per  second  [  1  ] ,  [2] .  The  A  /  D  converter  is  based  on  a 
set  of  superconducting  interferometers  [3]  which  operate  as 
current  comparators.  Each  of  these  interferometers  consists 
of  two  or  more  Josephson  junctions  connected  by  supercon¬ 
ducting  loops.  The  circuit  for  an  interferometer  is  shown  in 
Fig.  lc  and  its  I-V  curve  in  Fig.  la.  For  currents  below  the 
maximum  threshold  current,  I«„,  the  device  has  zero  voltage 
drop.  When  I™,  is  exceeded  the  interferometer  switches  very 
rapidly  (10-30  ps  depending  on  overdrive)  along  the  load  line 
to  the  energy  gap  voltage  (Vc  s  2.9  mV).  As  shown  in  Fig. 
lb,  the  threshold  current  I„  has  a  maximum  value  Im  and  a 
periodic  dependence  on  a  control  current  L  which  is 
magnetically  coupled  to  or  directly  injected  into  the  in¬ 
terferometer.  Thus,  by  applying  a  signal,  L,  to  the  control 
input  and  pulsing  the  interferometer  current  If  to  a  pre¬ 
determined  level,  (see  Fig.  lb),  the  interferometer  can  be 
used  as  a  comparator.  As  shown  in  Fig.  Id,  this  comparator 
has  the  unique  property  that  its  digital  output,  V,  is  a 
periodic  function  of  the  analog  signal  input,  L.  This  is  ex¬ 
actly  the  property  required  for  an  A/D  converter  since  each 
output  of  an  A/D  converter  is  a  periodic  function  of  the 
signal  input  with  the  period  changing  by  a  factor  of  two 
from  one  bit  to  the  next. 

Figure  2a  is  the  circuit  diagram  for  a  6-bit  converter  using 
this  idea.  The  required  variation  in  periodicity  is  achieved  by 
factor  of  two  changes  in  the  mutual  inductance  between  the 
common  signal  line  and  each  successive  interferometer. 
When  a  conversion  is  to  be  made,  a  current  pulse,  I„  is  ap¬ 
plied  to  the  power  supply  line  so  that  all  interferometers  are 
biased  to  the  level  shown  in  Fig.  lb.  Each  interferometer  for 
which  this  bias  level  exceeds  the  threshold  curve  will  switch 
to  the  voltage  state  producing  a  2.9  mV  pulse  on  its  output 
line.  These  pulses  are  a  digital  representation  of  the  analog 
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signal.  For  continuous  operation  the  supply  line  is  driven  by 
a  pulse  train,  or  at  the  highest  speeds  by  an  offset 
(monopolar)  sine  wave.  This  A/D  converter  operates  in  a 
fully  parallel  fashion  and  requires  only  one  comparator  for 
each  bit  of  accuracy.  Figure  2b  is  a  photo  of  the  circuit 
which  was  made  using  an  eight  level  fabrication  procedure 
described  elsewhere  [4]. 

Figure  3  shows  the  six  bit  patterns  as  a  function  of  analog 
signal  level.  Note  that  the  output  is  in  Gray  code  so  that 
perfect  alignment  of  the  threshold  levels  is  not  required.  In 
Fig.  3b  the  Gray  code  is  converted  back  to  an  analog  signal. 
The  monotonic  staircase  verifies  that  the  accuracy  of  the 
converter  is  ±  1  LSB  (least  significant  bit). 

Crosstalk  between  the  chip  leads  causes  considerable  dif¬ 
ficulty  in  testing  this  device  at  very  high  speed.  For  this 
reason  a  mount  has  been  designed  in  which  the  pads  on  the 
chip  are  pressed  against  the  ends  of  8  differential  microstrip 
line  pairs  which  are  etched  on  a  G-10  circuit  board.  Six 
millimeters  away  from  the  chip,  each  stripline  pair  converts 
to  SO  Q  coax.  The  coax  shields  are  tied  together  only  at  the 
superconducting  ground  plane  on  the  chip. 

Even  with  these  measures,  crosstalk  from  the  high  level  ac 
power  supply  (30  mV)  exceeds  the  output  signal  levels  (2.9 
mV)  at  sample  rates  above  2  x  109  per  second.  We  are  cur¬ 
rently  not  able  to  fully  verify  the  operation  of  the  converter 
at  this  rate.  However,  by  synchronizing  a  signal  waveform 
with  the  2  GHz  power  supply,  it  is  possible  to  use  a  sampling 
oscilloscope  to  observe  the  digital  data  on  the  output  lines. 
Figure  3c  shows  the  data  on  one  of  the  six  output  lines  for 
eight  successive  samples  of  a  6  ns  wide  signal  waveform.  The 
crosstalk  has  been  electronically  subtracted  from  this  data. 
It  is  clear  that  the  full  utilization  of  this  A/D  converter  will 
require  either  a  greatly  improved  sample  holder  and  inter¬ 
face  electronics  or  some  data  rate  reduction  on  the  chip.  The 
ultimate  speed  achievable  with  this  type  of  converter  is 
determined  by  the  tum-on-delay  of  the  interferometers  and 
the  time  required  to  reset  them  to  the  zero  voltage  state  at 
the  end  of  each  cycle. 

The  Fust  A/D  converter  of  this  type  [1]  used  3-junction 
interferometers  similar  to  those  developed  for  logic  applica¬ 
tions  [5].  Small  fabrication  variations  can  cause  these 
devices  to  have  nonperiodic  threshold  curves.  For  this 
reason  the  new  design  uses  2-junction  interferometers.  The 
desired  mutual  inductance  ratios  of  1:2:4:8:16:32  are 
achieved  by  varying  the  length  of  signal  line  which  passes 
over  each  interferometer.  As  shown  in  Fig.  2b,  these  lengths 
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vary  from  a  small  fraction  of  a  turn  (MSB)  to  2  full  turns 
(LSB). 

In  order  to  ensure  large  operating  margins,  it  is  desirable 
to  have  the  depth  of  modulation  of  the  threshold  current 
(Fig.  lb)  as  large  as  possible.  This  modulation  depth  is  given 
by: 

Wmin)/U  a  *s-  >  I*. 

<bo  L 

where  <t>.  =  2.07  x  10',5Wb  =  2.07  pH  mA.  In  order  to 
drive  SO  0  output  lines,  I„.  has  a  minimum  value  of  about 
100  nA.  There  is  thus  a  trade  off  between  operating  margin 
which  improves  for  small  L  and  sensitivity  to  the  signal 
which  improves  for  large  L.  We  have  chosen  L  =  6.5  pH 
which  gives  I„  (min)/ 1™.  -  0.3.  The  resulting  full  stale 
signal  current  is  3.5  mA. 

Various  types  of  on-chip  crosstalk,  signal  reflections,  and 
spurious  resonances  must  be  eliminated.  For  this  reason  the 
entire  circuit  is  built  on  a  superconducting  ground  plane  and 
all  lines  are  matched  and  terminated.  The  damping  resistor 
Ro  *  2  0  minimizes  a  resonance  between  the  interferometer 
inductance  and  the  junction  capacitance  [6].  The  power  sup¬ 
ply  isolation  resistor  R  must  also  be  large  enough  to  prevent 
a  transition  in  one  interferometer  from  significantly  affec¬ 
ting  its  neighbors.  The  minimum  value  of  R  to  prevent  LSB 
errors  is  found  to  be 

R  2  2pk  £  1  +  V  1  T j-‘ 

where  N  is  the  number  of  bits  and  Z,  is  the  power  supply 
source  impedance.  For  I_  =  100  pA,  Vc  =  2.9  mV,  N  = 
6.  and  Z,  =  50  Q,  R  a  400  Q.  Since  the  interferometer  bias 
levels  are  approximately  75  V#  of  I»,  the  peak  supply  voltage 
is  30  mV  and  the  average  supply  dissipation  is  about  9  pW 
for  the  entire  converter. 

In  order  to  achieve  monotonic  operation  with  6-bit  ac¬ 
curacy  it  has  proven  necessary  to  compensate  for  circuit 
variations  by  making  adjustments  to  both  the  on-off  ratio 
(signal  range  for  which  data  =  1  divided  by  signal  range  for 
which  data  »  0)  and  “phase”  of  each  of  thee  6-bit  curves 
(Fig.  3a).  These  12  adjustable  parameters  can  be 
automatically  set  by  simple  room  temperature  electronics  as 
follows:  With  the  signal  line  driven  by  a  full  scale  triangle 
wave,  the  on-off  ratio  of  each  bit  curve  is  measured  by  sim¬ 
ple  averaging.  Any  deviation  from  unity  is  used  to  control  a 
dc  current  which  is  injected  back  into  the  output  line.  This 
changes  the  peak  supply  level  for  that  bit  (i.e.,  the  intersec¬ 
tion  of  the  dashed  line  and  the  threshold  curve.  Fig.  lb)  and 
thus  corrects  the  bit  curve  on-off  ratio.  Six  of  these  servo 
loops  are  used  to  achieve  unity  on-off  ratio  for  each  in¬ 
terferometer.  It  can  be  shown  that  when  the  Gray  code  out¬ 
put  is  converted  to  binary  code  a  “phase”  error  in  Gray 
code  becomes  an  on-off  ratio  error  in  binary  code.  Thus,  a 
similar  :et  of  averaging  circuits  on  the  binary  code  outputs 
can  be  used  to  drive  dc  currents  into  an  additional  bias  con¬ 
trol  line  (see  Fig.  2)  which  overlaps  each  interferometer.  The 


currents  in  these  lines  shift  the  bit  curves  to  the  right  or  left 
to  achieve  the  proper  phase.  This  automatic  circuitry  can 
make  the  12  adjustments  in  a  fraction  of  a  second.  The  bias 
levels  are  then  held  constant  while  data  conversion  takes 
place.  At  present,  DC  drifts,  temperature  variations  or 
changes  in  local  magnetic  field  require  the  adjustment  to  be 
repeated  several  times  per  hour. 
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Fig.  1  (a)  The  l-V  curve  of  a  superconducting  interferometer,  (b)  the 
dependence  of  the  critical  current  L.  on  a  control  current  L,  (c)  the  cir¬ 
cuit  diagram  of  a  superconducting  interferometer,  and  (d)  the  input- 
output  characteristic  of  an  interferometer  comparator. 
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Fig.  2  (a)  The  circuit  diagram  for  a  superconducting  A/D  converter  and 
(b)  a  photograph  of  the  circuit. 


Fig.  3  The  six  digital  outputs  as  a  function  of  a  signal  output,  (b)  the 
digital  value  of  the  output  as  a  function  of  signal  level  (staircase  func¬ 
tion)  and  (c)  digital  output  data  on  one  of  the  six  output  lines  as  a  func¬ 
tion  of  time.  Eight  successive  samples  are  shown. 
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TURN-ON  DELAYS  IN  SINGLE  JOSEPHSON  JUNCTION  DEVICES 
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Introduction 

When  a  current  applied  to  a  device  containing  one  or  more 
Josephson  junctions  exceeds  the  critical  value  for  the 
development  of  a  voltage,  there  is  a  characteristic  delay 
before  the  onset  of  switching  to  the  energy-gap  voltage. 

This  is  called  the  turn-on  delay,  and  was  first  studied 
theoretically  by  E.  P.  Harris  (1).  Figure  1  is  a  simulation 
illustrating  a  typical  response  of  a  single  Josephson  junction 
to  a  step  increase  of  a  bias  current  Ife  through  the  junction 
critical  current  I  .  The  risetime  t_  and  turn-on  delay  t_ 
are  Shown.  Turn-on  delay  is  an  important  consideration  in 
the  design  of  high-speed  electronic  circuits  which  are 
required  to  switch,  as  in  A/D  converters  (2,  3) . 

Because  the  risetime  is  only  of  peripheral  concern  in  this 
paper,  we  shall  devote  just  a  few  comments  to  it  here,  and 
then  turn  our  attention  principally  to  turn-on  delay  in  the 
remainder  of  the  paper.  If  the  bias  current  does  not  change 
appreciably  during  the  switching,  an  expression  for  xR  is 
easily  obtained  (4)  from  the  recognition  that,  on  average, 
all  of  the  current  is  diverted  from  the  junction  supercur¬ 
rent  term  Icsin<t>  to  the  junction  capacitance  C.  The  risetime 
is  thus  approximately  the  time  required  to  charge  the  capaci¬ 
tance  to  the  gap  voltage  V  ,  and  is  about  CV^/Ic,  which 
reduces  to  4RnC/it  if  the  relation  Vg  =  4IcRNA  for  weak- 
coupling  superconductors  is  used.  Notice  that  it  is  the 


[ 


normal  state  resistance  R^,  rather  than  the  actual  resistance 
shunting  the  junction,  which  appears  in  this  expression. 

The  expression  is  valid  for  a  constant  Ifa  >  lc.  if  ifa 
changes  during  the  switching,  the  risetime  will  be  affected. 
An  expression  for  the  ramp-rate  dependence  of  tr  will  be 
given  in  the  last  section. 

The  turn-on  delay  time  is  defined  here  by  extrapolating  the 

<  rapidly  rising  portion  of  the  switching  curve  back  to  the 

time  axis;  the  time  interval  between  this  intersection  and 

the  time  at  which  I.  crosses  I  defines  i_.  The  extra- 

b  c  D 

polation  is  not  very  precise  because  of  the  shape  of  the 
rising  curve  (see  Fig.  1)  ,  and  this  makes  the  measured  value 
of  td  uncertain  typically  by  several  percent.  Hence  any 
expression  for  rD  need  not  be  very  precise.  This  is  true 
also  because  one  is  not  ordinarily  interested  in  a  precise 
;  measurement  of  it,  but  only  a  good  estimate.  This  is  the 

spirit  of  the  work  described  here. 


In  the  following  section,  we  analyze  the  turn-on  delay 
resulting  from  abruptly  stepping  the  bias  current  from  below 
to  above  the  critical  current  of  the  Josephson  device.  We 
consider  both  a  "bare"  junction,  i.e.,  a  resistively  and 
capacitively  shunted  junction  (the  RCSJ  model) ,  and  the  RCSJ 
shunted  by  an  inductor,  i.e.,  a  single- junction  SQUID. 

Exact  expressions  for  the  turn-on  delay  time  are  derived  in 
the  limit  of  no  damping,  and  compared  to  approximate  formulas 
derived  earlier  (4) .  In  the  third  section  we  develop  approxi¬ 
mate  expressions  for  the  ramp-rate  dependence  of  both  rD  and 
t  .  Since  this  case  more  closely  approximates  the  actual 

A 

operation  of  Josephson  switching  devices,  the  results  should 
be  of  some  practical  importance. 


A.  Stepped  Current 


Figure  2a  shows  the  circuits  we  consider.  The  indicated  re¬ 
sistance  R (V)  may  be  the  parallel  combination  of  an  external 
shunt  resistance  and  the  quasiparticle  resistance  illustrated 
in  Fig.  2b.  The  latter  figure  also  shows  a  load  line  inter¬ 
secting  the  quasiparticle  curve* above  the  lower  corner,  the 
voltage  at  this  corner  being  labelled  V  .  If  the  resistance 
is  so  low  that  the  load  line  intersects  below  V  ,  the  junction 
will  not  switch  to  the  gap  voltage.  This  case  is  not  of  con¬ 
cern  to  us  here.  For  higher  resistance,  it  turns  out  that  the 
resistive  term  has  relatively  little  effect  on  the  turn-on 
delay  time.  This  is  because  the  resistive  term  is  analogous 
to  viscous  drag,  which  is  proportional  to  velocity.  As'  long 
as  the  velocity  (voltage)  is  sufficiently  small,  this  force 
is  relatively  small.  By  definition,  the  voltage  is  small 
during  the  turn-on  delay  period,  and  thus  the  resistive  term 


is  not  very  important  for  turn-on  delay. 


a) 


Fig.  2a.  The  two  models 
considered  are  the  re- 
sistively  and  capaci- 
tively  shunted  junction 
(RCSJ) ,  and  the  RCSJ 
with  an  inductive  shunt, 
forming  a  SQUID. 


Fig.  2b.  The  quasi¬ 
particle  I-V  curve  is 
approximated  by  the 
three  legged  model.  A 
load  line  is  shown 
intersecting  the  I-V 
curve  in  the  region  of 
the  gap. 


The  equations  describing  the  RCSJ  model  are 


Ib  "  CV  + 


V/R  +  Icsin 4>, 


(1) 


V 


(2) 


where  the  dot  means  d/dt,  and  T  =  $  /2tt  =  )A/2e.  Earlier,  we  (4) 

o  o 

and  Harris  (1)  used  these  full  equations,  approximated  the  sin<t> 
term,  and  solved  the  resulting  equation  for  the  time  required 
for  the  phase  to  advance  by  a  specified  amount.  We  showed  (4) 
that  an  advance  of  1/2  radian  beyond  tt/2  was  generally  an  ap¬ 
propriate  amount.  Our  result  (4)  for  a  bias  current  of 

I  (1-6.)  for  t  <  0,  stepped  to  I  (1  +  6  )  at  t  =  0  and 
ci  c  ° 

thereafter  held  constant,  was 

1/2 


x  fl  +  2  (tt/2  -  4>±)  * 

td  =  aT  L  60  +  . 


(3) 


where 


(4) 


The  initial  value  of  the  phase,  4>^,  is  related  to  the  bias 
current  at  t  <  0  by 


sin  <|k  =  1  -  6^.  (5) 

Equation  (3)  was  obtained  after  dropping  all  terms  containing 
R,  which  are  typically  unimportant. 

We  have  since  realized  that  the  V/R  term  in  Eq.  (1)  can  be 
dropped  at  the  outset,  and  the  resulting  equation  solved 


exactly,  provided  1^  is  constant.  Thus,  multiplying  through 
by  <p  allows  us  to  write 


d 

Ht 


2  2 

-  ai  cos <b  -  oi 

P  P 


0, 


(6) 


or,  since  <t>  is  initially  zero  for  the  stepped-current  case. 


.2 

$ 


(7) 


where 


=  cos  $+  (1  +  6q)4> 


(8) 


Hence 


If,  as  in  the  earlier  derivation  of  the  approximate  formula 
(3),  we  choose  the  phase  at  which  the  voltage  begins  ac¬ 
celerating  significantly  to  be  1/2  radian  beyond  tt/2,  we  get 


T 


D 


2-l/2 


cu 

P 


+ 


1 

1 


-1/2 


d$ ' 


')-  f  (* 


i>] 


(10) 


We  have  not  found  a  way  to  evaluate  the  integral  analytically. 
However,  since  it  is  independent  of  all  junction  parameters,  it 
need  be  numerically  evaluated  only  once  for  a  given  overdrive 
or  "underdrive"  6^.  In  Figure  3  we  show  as  a  function 
of  SQ  for  two  values  of  namely  0.0  and  0.1.  The  solid 
lines  are  the  results  of  Eq.  (10),  and  the  dashed  lines 
represent  the  approximate  result  of  Eq.  (3) .  The  circles 
are  taken  from  simulations  as  shown  in  Fig.  1  and  discussed 
above.  Agreement  is  generally  excellent  for  the  smaller 
overdrives.  At  the  larger  overdrives,  the  overestimate  is 
probably  because  the  final  phase  <t> ^  used  to  evaluate  td 
should  not  be  a  fixed  quantity.  However,  the  inaccuracies 
are  not  serious.  The  simple  expression  (3)  is  probably 
sufficient  for  all  practical  purposes.  It  can  be  obtained 
from  Eq.  (10)  by  replacing  cos<*>  by  the  first  two  terms  of 
its  Taylor  expansion  cos  <|k  -  ( (f>  -  (fc^sinijK  over  the  entire 
region  of  integration. 

Next  we  turn  our  attention  to  a  single- junction  SQUID,  illus¬ 
trated  by  the  addition  of  the  inductor  shown  in  Fig.  2a.  This 
SQUID  is  described  by  Eq.  (1)  with  $  4>/L  added  to  the  right 
side.  In  its  usual  mode  of  operation  such  a  SQUID  does  not 
switch  to  a  finite  voltage  state.  There  will  be,  however,  a 
turn-on  delay  before  the  entry  of  a  flux  quantum.  An  A/D 
converter  based  on  this  type  of  SQUID  has  been  proposed  (5, 

6) .  This  turn-on  delay  is  easily  evaluated  along  the  same 
lines  as  above. 

The  "critical"  value  of  the  bias  current  before  flux  can 
enter  the  SQUID  is  not  that  of  the  junction,  because  of  the 
shunting  inductance,  but  rather  \ 


D1/2]  , 


cos-1  (3"1)+  (g2 


(11) 


FRACTIONAL  OVERDRIVE,  8Q 


Curves  of  <o  td  versus  fractional  overdrive  <S  for 
two  values  of  the  fractional  underdrive  Sefined 
hY  Ijj  “  IcCl-5iI.  The  solid  curves  are  fjom  Eq. 
CiOI/  and  the  clashed  curves  are  from  the  approximate 
Eq.  (31.  The  circles  are  values  of  x_  taken  from 
simulations,  as  in  Pig.  1.  u 


where 


6 


LW 


(12) 


(Mathematical  details  may  be  found  in  ref.  (7).)  We  assume 
that  1^  is  stepped  from  1^  (1  -  6^)  to  1^  (1  +  <SQ)  .  Dropping 
V/R,  we  may  cast  Eq.  (1)  (with  4>Q<fi/L  included)  into  the  form 

2 

LC  $  =  g(4>)  -  g  (<i> . )  ,  (13) 

2  1 

where 

g(4>)  »  6cos0  -  ■j  ($  -  8^)  ^  ,  .(14) 


8, 


LIb(t  >0)  LI®  (1  +  SQ) 


(15) 


The  turn-on  delay  is  then  given  by 

.♦f 


♦f  -1/2 

T  =  1/2  J  d*[g(*)  -  g  (*..)]  ,  (16) 


where  is  yet  to  be  specified.  The  initial  phase  is 
determined  from 


4k  +  $sin4>i  =  LIb  (1  -  Si)/^>Q. 


(17) 


To  choose  the  final  phase  4>f,  we  must  decide  first  how  to 
define  for  a  SQUID.  Figure  4  shows  the  response  of  a  6  - 

20  SQUID  to  a  bias  current  stepped  from  l£  to  1.01  l£.  The 
six  small  voltage  oscillations  in  each  large  voltage  excursion 
seen  in  the  figure,  corresponds  to  the  entry  (exit)  of  six 
fluxoids  into  (out  of)  the  SQUID.  The  slight  dampling  seen  is 
because  of  the  200  ft  subgap  resistance  assumed.  Perhaps  the 
only  logical  way  to  define  td  is  to  choose  it  at  the  maximum 
voltage  corresponding  to  the  first  flux  quantum  entry.  Now, 
dV/dt  =  0  implies  d<fi/d<p  =  0,  which  in  turn  requires  dg/dtp  =  0 
from  Eq.  (13)  .  Thus,  with  Eq.  (14)  ,  is  the  solution  of 

<J>0  +  Bsin  4>f  =  Ll£(l  +  6q)  .  (18) 


As  before,  the  integral  in  Eq.  (16)  cannot  be  determined 
analytically.  But  whereas  the  earlier  integral  was  indepen¬ 
dent  of  junction  parameters,  this  integral  depends  upon  the 
SQUID  6.  Figure  5  shows  the  results  of  Eq.  (16)  for  6^  =  0 
and  0.1  for  8  =  20.  Agreement  with  values  of  td  measured 
from  simulations  such  as  in  Fig.  4  is  now  almost  exact 
because  the  criterion  for  is  fixed  to  a  definite  feature 
on  the  voltage  curves.  The  result  of  neglecting  the  damping 
term  in  the  formula  for  is  not  observable. 

Analogously  to  the  preceding  case,  we  can  extract  an  ap¬ 
proximate  expression  from  Eq.  (16).  By  expanding  g(<p), 
defined  in  Eq.  (14),  about  its  initial  value,  one  obtains 


1/2 


(19) 


VOLTAGE  (mV)  ACROSS  JUNCTION 


Fig.  4 


Response  of  a  single  junction  SQUID.  The  junction 
of  Fig.  1  is  placed  in  a  superconducting  loop  of 
L  =  19.75  pH.  making  8  =  20.  I.  is  stepped  from 
1°  to  1.01  I®.  D 


8i=o.o 


8i=o.i 
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Pig.  5.  Curves  of  co  t_  versus  overdrive  6  for  6.  =  0.0  and 
Q.l.  The  solid  curves  are  from  Eq.  tl6)^  and  the 
dashed  curves  from  the  approximation,  Eq.  (19) .  The 
circles  are  the  results  of  simulations.  The  dis¬ 
continuities  reflect  the  vanishing  of  the  first  local 
voltage  maximum,  its  place  being  taken  by  the 
(until  then)  second  local  maximum.  (An  increasing 
oyerdriye  distorts  the  potential  energy  curve  in 
such  a  way  that  the  earlier  local  maxima  become 
less  pronounced,  eventually  becoming  only  ripples 
on  the  envelope  without  being  true  maxima.) 
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B.  Ramped  Current 


In  the  foregoing,  we  have  assumed  that  the  bias  current  was 
stepped  instantaneously  through  the  critical  current.  Although 
physically  one  can  probably  effectively  achieve  this  (e.g., 
one  should  be  able  to  develop  a  1  mA  current  in  a  picosecond 
by  a  number  of  methods) ,  there  will  be  cases  in  which  one 
would  ordinarily  be  ramping  more  slowly.  We  now  examine  the 
dependence  of  turn-on  delay  upon  the  ramp  rate  of  1^. 

Let  1^  vary  according  to 


Xb(tl  =  IcC1  +  cot)  '  (20) 

so  that  »  Ic  at  t  =  0.  We  again  omit  the  V/R  term  in  Eq. 

(11.  If  <*>  (t=Ql  is  near  enough  ir/2  that  sin$  (o)  =  1,  Eq. 

Cl  I  may  be  written  j 

=  to2  cot  .  (21) 

r 

This  is  integrated  twice  to  give 

<Ktl  =  $(0)  +  i  COlt  +  co 2  cot3/6  .  (22) 

F 

When  CO  1  and  $C0)  are  known,  and  a  suitable  final  rhase  <p^  for 
evaluating  T^is  picked,  can  be  found  by  solving  the  cubic  Eq. 

C22L. 

At  this  point  we  need  to  distinguish  three  cases.  (1)  The 
ramp-rate  is  suitably  high.  By  this  we  mean  that  the  in¬ 
tended  overdrive  has  been  reached  before  the  system  has 
sensibly  responded.  As  will  be  seen  below,  or  as  may  be 
surmised  from  Eq.  (22),  the  criterion  for  this  case  is 


oi  >>  oip.  (2)  The  ramp-rate  is  suitably  small.  By  this  we 

mean  not  only  a  <<  a  ,  but  also  that  the  switching  is  com- 

P 

pleted  before  the  ramping  stops.  (3)  The  intermediate  case, 
in  which  the  ramp  rate  cannot  be  considered  a  step  function, 
but  in  which  the  intended  overdrive  is  reached  before  switching 
is  finished. 


For  case  (11,  the  formulas  of  the  preceding  section  apply,  and 
we  proceed  to  case  (2) .  Since  here  the  ramp  rate  is  small,  we 
can  expect  <K0)a  n/2,  i.e.,  that  sin<J>  should  track  well  with 
I./I  .  For  t  <  0,  while  4>  is  small,  $  ^  o>;  as  t  approaches 
zero,  <t>  increases  but  (0)  will  still  be  of  the  order  of  o>. 

Thus  Eq.  (221  becomes 

_ 

at  (a+4wV)  =  <P  -  <H0),  (23) 

6&T 


where  a  is  somewhat  larger  than  unity.  Since  oi  >>  os  in  the 

2  2^.  3 

present  case,  it  will  be  true  that  cu  /6co  >>  a,  and  the  t 

P  n 

term -dominates  for  4>  — 4>  CO  1  <  1 .  Finally,  if  we  pick  + 

i,  as  the  preceding  section,  as  the  phase  at  which  to  evaluate 
td,  we  obtain 


T 


(24) 


Equation  (241  overestimates  xD  since  the  linear  term  in  Eq. 
C23L  was  dropped.  Keeping  the  linear  term  and  setting 
a  =  1,  one  gets 


T 


G1/3 


[(S+l)1/3-(S-l)1/3], 


(25) 


where 


G  =  3cop/2co  ,  C261 

2  1/2 

S  =  (1  +  8/G  J  .  (.27) 

In  Fig.  6  we  compare  Eqs.  (24)  and  (25)  with  the  results  of 
simulations.  Equation  (25)  is  seen  to  be  better,  but  the 
simple  Eq.  (24)  is  probably  good  enough  for  most  purposes. 

The  fact  that  the  turn-on  delays  for  the  slowest  ramp  rates 
do  not  converge  toward  Eqs.  (24)  or  (25)  very  well  is  probably 
a  reflection  of  the  fact  that  a  fixed  is  not  the  best  way 
for  evaluating  xD,  as  mentioned  also  in  the  preceding  section. 

To  complete  case  (2) ,  we  examine  the  corresponding  risetime 
tr.  For  t  >  td,  sin<f>  oscillates  more  and  more  rapidly,  so  that 
on  average  it  should  be  negligible.  Then  Eq.  (1)  becomes 


CV  -  ICU 


+  cot) 


(28) 


Integrating  this  for  t  >  xD,  setting  V(td)  =  0,  and  defining 
xR  by  the  time  after  xQ  at  which  V  reaches  V  ,  we  get 


TR  = 


1  +C0X- 


CO 


2coCV  /I  \ 
_ SL _ £_  J 

(l+coxD>2  / 


1/2 


(29) 


This  formula  is  found  to  be  within  10%  of  the  estimated  Xj 

from  simulations  over  the  ramp-rate  range  co  =  0.005  cop  to 

0.36  co  . 

P 


Comparison  of  the  turn-on  delay  time  x_  from  Eqs. 

C2 4 }  and  (251 ,  with  that  estimated  from  simulations 
Ccircles) ,  for  current  ramped  up  from  zero  at  dif¬ 
ferent  rates.  The  junction  is  the  same  as  that  in 
Fig.  1.  It  is  assumed  that  the  overdrive  is  not 
reached  before  the  completion  of  the  turn-on  delay 
period. 


Case  (3) ,  the  intermediate  case,  is  not  so  easily  evaluated, 
and  probably  not  worth  a  detailed  examination.  This  is  the 
case  in  which  the  ramp  rate  is  high  enough  that  4>{0)  in  Eq. 
(22)  is  not  well  estimated  by  H/2,  nor  <j>(0)  by  w.  Further 
this  is  the  case  in  which,  because  the  ramp  rate  is  in  the 
intermediate  range,  the  intended  overdrive  is  likely  to  be 
reached  in  the  midst  of  the  switching  process.  Formulas  can 
indeed  be  worked  out  but  they  are  not  simple  and  each  is  ap¬ 
plicable  only  for  a  narrow  range  of  overdrives  and  ramp 
rates.  Therefore  we  will  not  present  any  additional 
formulas.  It  is  sufficient  to  observe  that  once  the  ramping 
has  stopped,  the  system  is  no  longer  being  "pushed",  and  the 
turn-on  delay  time,  or  risetime,  whichever  has  been  interupted 
by  the  halted  reunping,  will  be  lengthened  over  what  the  con¬ 
tinuous  ramp  rate  formulas  would  indicate.  Since  the  simple 
ramping  result  of  Eq.  (24)  is  an  overestimate,  its  use  in 
the  intermediate  cases  is  usually  satisfactory,  as  simulations 
show. 

In  summary,  the  and  tr  expressions  of  the  stepped  current 
section  may  be  used  for  fast  ramp  rates,  i.e.,  for  co  >>  co^. 

For  slow  rates,  m  <<  co^,  Eqs.  (24)  or  (25),  and  (29)  are  good. 
In  the  intermediate  cases,  Eq.  (24)  will  be  satisfactory  for 
td,  and  Eq.  C291  will  overestimate  tr  by  an  acceptable  amount. 
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Picosecond  Applications  of  Josephson  Junctions* 
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Abstract 

The  behavior  of  simple  superconducting  circuits  In  the  picosecond 
regime  Is  described  In  a  comprehensive  way,  with  primary  emphasis  being 
given  to  the  step  function  and  pulse  responses  of  these  circuits.  Topics 
receiving  detailed  discussion  Include  Josephson  junction  modeling  with 
both  the  microscopic  and  shunted  junction  models.  Limitations  of  the 
shunted  junction  model  are  explored  by  comparing  It  with  experimental 
results  and  with  the  microscopic  model.  An  approximate  evaluation  is 
given  of  the  Important  dynamical  properties  of  junctions  made  with  the 
dominant  fabrication  technology  (Pb  alloy  systems),  as  a  function  of 
tunneling  barrier  thickness.  Rounding  out  the  device  aspects  of  the 
discussion,  we  describe  in  detail  the  properties  of  superconducting  micro¬ 
strip  transmission  lines,  with  an  emphasis  on  their  high  speed  behavior. 
Turning  to  simple  circuits  we  review  experimental  results  on  the  measurement 
of  picosecond  regime  transient  signals.  The  concept  of  turn-on  delay  is 
analyzed  anew,  providing  simplified  and  extended  results.  Details  of  concepts 
for  pulse  height  and  pulse  width  measurements  are  explored,  leading  to  the 
conclusion  that  the  time  resolution  of  superconducting  circuits  is  limited 
to  approximately  the  period  of  one  plasma  oscillation.  With  present  Pb  alloy 
fabrication  technology  this  limit  Is  2  ps. 
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Picosecond  Applications  of  Josephson  Junctions 

I,  Introduction 

Josephson  junctions  are  now  being  used  for  the  generation,  detection, 
and  measurement  of  pulses  and  switching  transients  with  durations  as 
short  as  about  10  ps.  Our  purpose  here  Is  to  present  a  unified  picture 
of  the  more  basic  developments  In  this  area,  including  some  new  results. 

The  earliest  developments  were  the  microscopic  theories  of  superconducting 
tunnel  junctions  by  Riedel  (1)  and  Werthamer  (2),  which  predicted  essentially 
a  flat  spectral  response  out  to  the  frequency  of  the  energy  gap,  where  a 
somewhat  enhanced  response  was  expected.  At  frequencies  above  the  gap 
the  superconducting  response  goes  monotonically  to  zero,  but  with 
measureable  effects  extending  out  to  several  times  the  gap.  Experimental 
investigations  of  the  bandwidth  of  Josephson  junctions  were  first  carried 
out  In  the  frequency  domain,  perhaps  most  Importantly  at  891  GHz  (30% 
above  the  gapl,  where  It  was  found  that  the  response  was  Indeed  close  to 
that  at  low  frequencies.  These  experiments  were  done  with  point  contact 
devtces  for  which  It  was  known,  from  of  the  analyses  of  Stewart  (3)  and 
McCumber  (4),  that  the  capacitance  was  quite  small. 

In  the  early  years  of  the  Josephson  effect  it  was  commonly  thought 
that  the  response  time  of  tunnel  junctions  would  be  severely  limited  by 
the  capacitance  associated  with  the  exceedingly  thin  tunneling  barrier. 
Basavatah,  Eldrldge,  and  Matlsoo  (5)  corrected  this  misconception  by 
experimentally  evaluating  phenomenological  tunneling  theory  and  producing 
for  the  first  time  a  clear  picture  of  the  relationship  between  tunneling 
barrier  thickness  and  current  density.  From  this  It  was  only  a  small 


step  to  relate  current  density  to  junction  response  time,  from  which  It 
became  clear  that  the  picosecond  regime  should  be  accessible  with  tunnel 
junctions. 

The  evolving  development  of  high  quality  superconducting  integrated 
circuits,  discussed  In  detail  In  this  special  Issue,  creates  the  possibility 
of  practical  utilization  of  nearly  1  THz  of  bandwidth  In  electronic  systems. 
In  the  following  we  will  more  thoroughly  examine  the  experimental  and 
theoretical  basis  for  the  quantitative  views  expressed  here.  We  will  also 
present  new  results  regarding  turn-on  delay,  pulse  detection,  and  pulse  width 
measurement,  with  an  emphasis  on  performance  at  the  highest  possible  speeds. 
The  order  of  discussion  which  follows  Is:  II.  EXPERIMENTAL  RESULTS, 

III.  THEORETICAL  BASIS,  IV.  SWITCHING  TIME,  V.  PULSE  GENERATION, 

DETECTION,  AND  MEASUREMENT,  VI.  TRANSMISSION  LINES,  and  VII.  DISCUSSION. 

Our  discussion  emphasizes  single  junction  devices  rather  than  multi junction 
SQUIDs  both  because  the  former  are  simpler  to  analyze  and  to  fabricate  and 
because  they  probably  give  the  fastest  response,  though  the  latter  point 
has  not  been  proven.  We  have  also  chosen  to  largely  omit  discussion  of 
nonlatching  devices  (5]  and  of  the  resetting  (or  punchthrough]  problems 
{7L  Finally,  discussion  of  noise  was  omitted  stnce  fundamental  physical 
noise  will  not  ordinarily  dominate  the  circuits  of  Interest.  More  germane 
to  the  Irregularities  of  practical  ctrcuits  Is  the  problem  of  reflections 
from  tmpedance  mismatches,  which  are  nearly  always  present. 

II.  EXPERIMENTAL  RESULTS 

Although  most  of  our  discussion  will  be  about  the  time  domain  response 
of  superconducting  circuits.  It  is  appropriate  to  begin  with  a  brief  review 
of  experiments  in  the  frequency  domain,  both  because  they  gave  the  first 
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experimental  evidence  of  high  frequency  response  and  because  time  domain 
measurements  are  still  limited  to  times  greater  than  about  10  ps,  whereas 
frequency  domain  measurements  using  far  infrared  lasers  probe  the  sub¬ 
picosecond  regime. 

A.  Bandwidth  Measurements  In  the  Frequency  Domain 

The  most  direct  proof  of  the  large  bandwidth  of  Josephson  Junctions 
comes  from  a  rather  simple  point  of  view.  We  take  the  low  frequency 
description  of  the  Josephson  effect  and  ask  at  what  frequency  does  It 
break  down?  Fortunately  Josephson  junctions  have  a  special  property 
which  allows  this  to  be  done  with  little  uncertainty. 

The  low  frequency  form  of  the  Josephson  equations  relating  current 
I  to  voltage  V  ts: 

I  -  I0  *1n  ♦  ,  (1) 

♦  •  V  ,  (2) 

where  the  dot  means  d/dt  and  e  and  H  are  the  usual  fundamental  constants. 
These  equations  predict  that  if  a  superconducting  junction  is  dc  biased 
at  voltage  Y&,  the  phase  $  will  advance  uniformly  In  time  and  the  current 
I  will  oscillate  with  amplitude  IQ  at  frequency 

fj  "  Hi”  Vb  "  484  SH2/mV)-  (3) 

A  more  complete  theory  of  Junctions,  which  will  be  discussed  later, 
makes  tt  clear  that  It  Is  the  quantity  IQ,  the  junction  critical  current. 


which  eventually  begins  to  decrease  with  frequency.  That  Is  not  at  all 
surprising  since  one  would  not  expect  to  be  able  to  produce  significant 
optical  frequency  current,  for  example,  from  a  state  of  matter  which  has 
a  binding  energy  of  only  a  few  mV. 

The  experimental  problem  then  Is  to  measure  IQ  as  a  function  of 
frequency.  Since  the  frequency  range  of  Interest  exceeds  1  THz  It  Is  a 
difficult  practical  problem. 

Brian  Josephson,  In  his  original  paper  on  superconducting  tunneling  (8), 
suggested  an  indirect  method  for  observing  the  supercurrent  I  at  microwave 
frequencies.  This  method,  first  successfully  used  by  Shapiro  (9),  Is  based 
on  the  simultaneous  application  to  a  junction  of  dc  bias  and  microwave 
radiation.  The  basic  Idea  is  that  for  proper  adjustment  of  Vfa  the 
frequency  of  the  Josephson  oscillation  fj  will  be  equal  to  the  frequency 
f.  of  the  applied  radiation.  Since  the  Josephson  junction  Is  a  nonlinear 

a 

device.  It  will  produce  a  beat  frequency  current  at  the  difference 
frequency  f,  -  f  ,  which  in  this  case  Is  at  zero  frequency.  The  result 
Is  a  current  step  on  the  I-V  characteristic.  This  situation  can  be 
analyzed  using  Eqs.  (1)  and  (2)  with 

V  ■  Vb  +  v!  cos(2irfat  +  0^.  (4) 

From  this  analysis  it  is  found  that  the  Induced  dc  Is: 

'<lcCt0  ■  fa>  •  2Io  V2,Vf4>'  <5> 
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where  J-j  is  the  first  order  Bessel  function  of  the  first  kind.  Thus  Idc> 

considered  as  a  function  of  V,,  has  a  maximum  value  of  1.16  I  .  The 

i  o 

anaylsis  given  here  omits  the  shunt  capacitance  and  shunt  resistance  of 
real  junctions.  Simulations  by  Russer  (10)  and  Kautz  (11)  have  shown  that 
the  induced  dc  given  by  £q.  (5)  Is  an  upper  bound  which  Is  closely 
approached  In  the  right  circumstances.  The  upper  bound  suffices  for  our 
purposes  so  we  will  not  consider  In  detail  the  effects  of  the  shunts. 

Experiments  of  this  type  specifically  designed  to  explore  the 
bandwtdth  of  the  Josephson  effect  were  reported  several  years  ago  by 
McDonald,  et  al.,  (12,  13)  and  more  recently  by  Weitz,  Skocpol ,  and  Tlnkham  (14). 
Perhaps  the  most  direct  evidence  for  large  amplitude  Josephson  oscillations 
at  a  frequency  In  the  vicinity  of  the  energy  gap  Is  Illustrated  (12)  In 
fig.  1.  The  experimental  technique  which  was  used  for  acquiring  this  data 
largely  removed  the  background  I-V  curve  from  the  data  and  emphasized 
the  Induced  dc.  This  was  an  experiment  In  which  a  Josephson  junction 
was  simultaneously  Irradiated  at  891  and  805  GHz.  As  the  voltage  across 
the  junction  was  slowly  Increased,  the  frequency  of  the  Josephson  oscillation 
Increased  according  to  Eq.  (3)  until,  at  a  voltage  of  1.66  mV,  this 
osctllatton  phase  locked  with  the  applied  far  Infrared  laser  radiation 
at  805  GHz  and  produced  a  zero  frequency  beat  current  which  Is  Illustrated 
In  the  figure.  As  the  voltage  was  Increased  further,  the  Josephson 
oscillation  phase  locked  with  the  891  GHz  radiation  at  1.84  mV  and  again 
produced  a  zero  frequency  beat,  but  this  time  of  greater  amplitude 
because  of  htgher  Intensity  radiation  (larger  V-|).  In  effect  the  Junction 
acts  as  a  spectrum  analyzer  of  unusually  large  bandwidth.  Our  main 
Interest  however  Is  In  the  result  that  the  junction  gives  a  large  response 
at  thts  extremely  high  frequency;  in  fact  the  magnitude  of  the  response  Is 


0.8  IQ.  Considering  this  result,  as  well  as  other  more  extensive  data  04), 
we  conclude  that  the  Josephson  response  is  <\*  705!  of  the  upper  bound  of 
Eq.  C5). 

The  principal  frequency  of  this  experiment,  891  GHz,  Is  somewhat 
above  the  energy  gap  of  675  GHz  of  the  superconductor  (Nb)  used  In  the 
experiment.  Thus  tn  an  approximate  way,  considering  the  possible  deleterious 
effects  of  shunt  impedances,  heating  (15),  and  noise  rounding  (16,  15),  we 
conclude  that  the  Josephson  effect  operates  at  nearly  full  amplitude  In 
the  vicinity  of  the  energy  gap,  a  very  high  frequency  for  electronics. 

More  generally  it  has  been  found  experimentally  that  I  decreases  above 
the  gap  frequency  such  that  by  approximately  12  times  the  gap  it  is  no 
longer  observable  above  the  noise  (12-14). 

These  high  frequency  experiments  were  done  with  point  contact  junctions 
to  accommodate  the  far  Infrared  radiation.  Closely  related  experiments  with 
tunnel  junctions  lead  to  the  same  conclusions  (18,19). 

Also  consistent  with  the  discussion  given  here  Is  the  study  of  self 
resonances  In  Josephson  interferometers  by  Zappe  and  Landman  (20).  They 
conclude  that  the  simple  model  of  the  Josephson  effect  Is  applicable  up  to 
frequencies  at  least  as  high  as  300  GHz. 


8.  Time  Domain  Experiments 

The  early  high  frequency  experiments  with  Josephson  devices  suggested 

their  use  as  very  fast  logic  switches.  Matlsoo  (21)  demonstrated  this 

potential  in  the  first  important  time  domain  experiment.  This  experiment 

placed  an  upper  limit  of  800  ps  on  the  transition  time  of  a  Josephson  tunnel 

junction  for  switching  from  the  zero  voltage  state  to  the  energy  gap  voltage. 

After  many  improvements  in  fabrication  and  testing,  Jutzi  (22)  published 

data  In  1972  which  put  an  upper  limit  of  34  ps  on  the  switching  time  of  a 
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small  area  junction  with  a  critical  current  density  of  1000  A/cm  .  The 
most  recent  published  data  for  Josephson  logic  devices  Indicate  switching 
times  of  about  7  ps  and  total  logic  delay  of  as  little  as  13  ps  per  gate 
(23). 

The  fact  that  Josephson  devices  can  produce  and  detect  these  exceedingly 
fast  signals  gives  them  the  potential  for  improving  the  state  of  the  art  In 
very  fast  waveform  analysis.  To  understand  how  this  might  be  done  we  can 
consider,  as  a  first  approximation,  a  Josephson  junction  to  be  a  device 
which  switches  from  zero  voltage  to  the  energy  gap  voltage  (^  3  mV)  at 
the  Instant  its  critical  current  is  exceeded.  This  model  is  discussed  more 
fully  in  Section  IV.  Since  the  transition  is  very  rapid,  a  Josephson 
junction  can  act  as  a  comparator  to  determine  the  time  at  which  some  threshold 
Is  crossed.  Also,  Josephson  devices  which  operate  In  a  latching  mode  can 
be  used  to  detect  and  remember  an  extremely  short  excursion  above  an 
arbitrary  threshold  level.  Both  of  these  properties  are  used  to  measure 
the  waveforms  of  repetitive  analog  signals. 

Figure  2(a)  shows  one  common  way  this  is  done  (24-26).  The  unknown 
repetitive  signal  Is,  together  with  a  bias  current  1^,  fs  applied  to  a 
tunnel  junction  Jl .  When  the  sum  of  these  two  currents  exceeds  the  critical 


current  IQ  of  Jl ,  It  will  switch  causing  a  voltage  transition  to  be 
transmitted  to  room  temperature  apparatus  where  its  time  of  arrival  is 
measured.  Since  the  bias  and  the  critical  current  of  Jl  are  known,  the 
signal  level  at  the  time  the  detector  switches  can  be  estimated.  Thus,  by 
measuring  the  delay  In  the  Jl  transition  as  a  function  of  bias,  rising  edges 
on  the  signal  may  be  mapped  out.  In  practice  this  may  be  automated  so  that 
the  unknown  signal  Is  traced  out  directly  on  an  X-Y  recorder.  Hamilton,  et  al . , 
(26)  used  this  method  to  observe  the  risetime  of  a  4200  A/cm  junction 
adjacent  to  the  sampler  (see  Fig.  2(b)).  The  two  traces,  recorded  about 
one  minute  apart,  are  for  Increasing  and  decreasing  bias.  They  are 
different  because  time  drift  In  the  room  temperature  electronics  decreases 
the  apparent  risetime  in  one  case  and  Increases  It  in  the  other.  Assuming 
constant  time  drift  the  correct  risetime  is  the  average  or  about  9  ps. 

This  result  represents  the  combination  of  the  actual  signal  risetime  and 
the  time-jitter  limited  system  refponse  of  about  7  ps.  Simulations  of 
this  circuit  show  that  there  Is  a  turn-on  delay  between  the  time  that  I 

o 

Is  exceeded  and  the  switching  transition  of  Jl .  This  turn-on  delay 
depends  on  the  degree  of  overdrive  of  Jl ,  among  other  factors  (see  Section 
IV).  When  Ifa  Is  high,  the  overdrive  is  high  and  the  turn-on  delay  is 
minimal.  As  1^  decreases  the  turn-on  delay  Increases.  The  result  Is 
that  signals  observed  with  this  technique  are  distorted,  with  the  observed 
signal  risetime  being  a  few  picoseconds  greater  than  the  actual  risetime. 

The  character  of  the  distortion  of  the  signal  Is  consistent  with  the 
theoretical  turn-on  delay  formula  Eq.  (23).  Whether  the  agreement  is 
quantitatively  accurate  Is  Impossible  to  say  because  of  the  amount  of  jitter 
present.  In  addition  to  the  jitter  problem  this  technique  is  further 
limited  in  that  it  can  only  observe  rising  edges  of  the  signal. 
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Both  of  these  limitations  may  be  removed  by  another  method  (27,  28) 
depicted  In  Fig.  3.  In  this  circuit  the  current  through  Jl  Is  the  sum  of  a 

dc  bias  and  a  fast  ramp  Ip  *  at.  Jl  will  switch  to  the  voltage  state 

when  *bl  +  *r  *  *ol *  producing  a  voltage  transition  at  a  time  T  *  ^ol^bl^0. 
This  voltage  transition  Is  differentiated  by  a  capacitor  C  to  produce  a  sampling 
pulse  of  maximum  value  1^.  This  pulse  Is  applied  to  junction  J2.  J2  also 
receives  current  form  the  signal,  I$ Ct ) ,  and  a  second  bias  source  Ib2 •  The 
maximum  values  of  Is  and  Ib2  are  arranged  so  that  Iq2  can  never  be  exceeded 
except  at  the  Instant  of  the  sampling  pulse.  Thus  if  Ip-|  +  Is (T )  +  >  IQ2» 

J2  will  switch  producing  an  output  voltage  VQ  3  Vg.  If  Ipl  +  I$(T)  +  Ib2  <  IQ2 

the  output  voltage  is  zero. 

The  output  voltage,  averaged  over  many  repetitions  of  the  signal,  is 
compared  with  a  value  Vg/2  and  the  amplified  difference  is  used  to  drive  Ib2. 
This  servo  loop  automatically  adjusts  Ip2  to  the  threshold  value  where  J2  is 
switching  half  the  time  and  Ip^  +  I$(T)  +  Ib2  *  I02*  Since  Ip1 ,  Ib2  and  Iq2 
are  known,  the  signal  IS(T)  Is  determined.  The  complete  signal  waveform  may 
be  recovered  by  sweeping  Ipl  and  plotting  the  computed  signal  (Iq2  -  Ipl  -  Ib2) 
versus  equivalent  time  (Iol  -  Ib|)/o.  Note  that  this  circuit  works  in  a 
repetitive  mode  where  Jl  and  J2  are  reset  to  the  zero  voltage  state  between 
each  signal  recurrence.  The  reset  circuits  are  not  shown  In  Fig.  3.  In 
practice  the  sampling  pulse  may  be  created  In  several  different  ways  and 
the  detector  junction  J2  may  be  a  Josephson  Interferometer.  Faris  (28)  and 
Tuckerman  (29)  have  used  this  principle  In  a  general  waveform  sampling 
circuit  with  an  estimated  time  resolution  of  6  ps. 

In  order  to  produce  an  undistorted  picture  of  the  signal  with  this 
technique,  the  maximum  value  of  IsCtl.+  Ip]  must  occur  at  the  same  time 
as  the  maximum  of  Ip^.  Since  Ip^  Is  rounded  on  top,  any  variation  of  I$(t) 


9 


across  the  rounded  top  of  Ip^  will  cause  the  maxima  not  to  coincide  and  thus 
lead  to  distortion  In  the  observed  signal.  The  magnitude  of  this  distortion 
depends  on  the  width  of  Ip1  and  the  amplitude  of  I$(t)  relative  to  Ipl .  In 
practice  this  error  can  be  reduced  to  a  few  picoseconds  (28).  A  second  source 
of  distortion  results  from  the  finite  response  time  of  the  detector  junction  J2. 
Discussion  of  this  point  will  be  left  to  Sections  IV  and  V  which  consider  in 
detail  the  dynamics  of  Josephson  junctions  used  as  latching  detectors. 
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III.  THEORETICAL  BASIS 


The  design  of  circuits  Involving  Josephson  junctions  is  facilitated  by 
the  availability  of  accurate  models  for  the  junctions.  The  most  commonly  used 
model  is  a  phenomenological  description  known  as  the  RCSJ  model,  or  reslstlvely 
and  capacltlvely  shunted  junction  model  first  described  by  Stewart  (3)  and 
McCumber  (4).  Following  a  brief  description  of  this  model  we  discuss  certain 
features  of  the  more  complex  microscopic  theory  on  which  it  Is  based. 

The  microscopic  theory  is  of  practical  Importance  because  It  provides  a 
basis  for  testing  the  phenomenological  model  in  areas  that  are  experimentally 
difficult,  such  as  the  picosecond  regime. 

A.  RCSJ  Model 

Stewart  and  McCumber  used  Eqs.  (1)  and  (2)  to  describe  supercurrent 
flow  in  a  Josephson  junction  and  added  a  resistive  shunt  to  account  for 
normal  current  at  finite  voltages.  A  capacitive  shunt  was  also  needed 
because  of  the  physical  nature  of  a  tunnel  junction.  Their  circuit 
model  is  Illustrated  In  Fig.  4  (a).  In  order  to  encompass  the  Idea  of 
switching  from  the  zero  voltage  state  to  the  energy  gap  voltage  V  , 
the  shunting  resistance  can  be  made  voltage  dependent,  as  illustrated  In 
Fig.  4  (b).  We  designate  this  junction  resistance  as  Rj. 

If  the  RCSJ  model  Is  driven  by  a  current  souce  I(t)  the  circuit 
equation  for  the  case  with  no  load  conductance  Is: 

I(t)  -  Cj$  +  V  Rj"1  +  I0sin<t>  .  (6) 
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Much  of  the  analysis  to  be.  discussed  here  will  involve  junction  behavior  at 


voltages  below  the  energy  gap,  in  which  case  Rj  *  R1  (see  Fig.  4  (b)). 
Restricting  analysis  to  this  region,  one  may  write  Eq.  (6)  in  the  form: 

^  *  “p  {rncj)2  ”£  +  S”  He  +  s1n**  (7) 

where  R^  is  the  junction  normal  state  resistance  and  wp  is  the  plasma 
frequency  (30): 


“p  5  <VCJ  *o>1/2  •  (8) 

where  5  2e.  The  dimensionless  time  9  =  t,  “g  3  2  ^oW^o' 

Equation  (7)  contains  the  three  time  scales  of  a  Josephson  junction  related 
to  oip,  RnCj,  and  to  tug.  We  note  incidentally  that  these  time  scales  are 
not  Independent,  In  fact  uip  R^Cj  *  oig. 

The  aig  time  scale  is  based  on  the  quantity  IQR^.  As  will  be  seen 
later  when  discussing  the  microscopic  theory  of  a  Josephson  junction, 
the  quantity  IQRN  *  ira/2e  where  t  is  the  energy  gap.  Thus  IQRN  is 
determined  by  a  microscopic  quantity  and  we  shall  therefore  refer  to  9  as 
the  microscopic  time.  On  the  other  hand  the  quantities  u>p  and  R^Cj  both 
depend  on  capacitance,  or  in  other  words,  on  the  presence  of  a  macroscopic 
dielectric  layer,  so  we  shall  refer  to  the  corresponding  time  scales  as 
the  macroscopic  scales.  These  three  time  scales  will  be  examined  in  more 
detail  in  the  following  two  sections. 
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B.  Microscopic  Theory 

Josephson,  in  his  original  work  on  superconducting  tunneling  (6), 
described  the  consequences  of  Eqs.  0)  and  C2).  It  remained  for  others 
to  evaluate  the  theory  for  frequencies  in  the  vicinity  of  the  energy  gap 
of  the  superconductors  making  up  the  Junction.  That  was  subsequently  done 
by  Riedel  (1),  and  more  completely  by  Werthamer  (2),  who  were  in  turn 
followed  by  other  significant  contributors  (31-34).  Because  the  microscopic 
theory  includes  detailed  effects  which  are  not  present  in  the  phenomenological 
theory  It  is  more  complex  and  less  straightforward  to  Interpret.  Here  we 
will  discuss  only  those  aspects  of  the  microscopic  theory  which  are  con¬ 
ceptually  relevant  to  picosecond  response.  Further  details  are  presented 
in  Appendix  A. 

Initially  the  microscopic  theory  was  presented  in  a  frequency  domain 
formulation.  This  theory  describes  only  the  tunneling  current  but  the 
capacitive  current  can  be  added  exactly  as  for  the  phenomenological 
model.  The  tunneling  current  Is  the  sum  of  two  currents:  one  of  normal 
electrons  called  the  quasiparticle  current  and  described  by  I^Cai)  and 
one  of  paired  electrons  called  the  supercurrent  and  described  by  IjC<*>). 
Although  both  I  and  Ij  have  real  and  Imaginary  parts,  which  are  all 
required  for  the  complete  theory,  it  is  Im  I  and  Re  Ij  which  provide 
the  theoretical  basis  for  the  phenomenological  theory.  These  functions 
are  shown  In  Fig.  5.  The  junction  critical  current  of  the  phenomenological 
model  Is  the  low  frequency  limit  of  Re  I ji 

I0  *  ‘  11m  Re  *  2efc  •  C9) 

u  -  o  N 

The  right  hand  equation  Is  correct  for  superconductors  with  weak  electron- 
phonon  coupling  but  requires  a  correction  factor  for  strongly  coupled 


materials  (35,  361.  The  frequency  dependence  of  IjGu)  is  the  frequency 
dependence  of  IQ,  which  was  discussed  in  Section  II.  The  peak  in  Re  Ij 
at  the  energy  gap  describes  an  internal  microscopic  resonance  called  the 
Riedel  peak  which  has  been  experimentally  verified  (18,  191. 

The  phenomenological  resistive  current  VRj1  in  Eq.  (6)  appears  in  the 
microscopic  theory  as  Im  I^fa).  When  the  full  microscopic  theory  is  used 
to  calculate  an  I-V. curve  (37)  the  abrupt  rise  in  Im  I  (w)  of  Fig.  5  occurs 
at  a  voltage  V  »  =  2A/e,  the  energy  gap  voltage.  Thus  the  phenomenological 

resistive  current  is  Im  I  (V).  Indeed  this  function  provides  a  highly  ac¬ 
curate  description  of  experimentally  measured  current-voltage  characteristics  (38). 

Tn  our  later  discussions  the  plasma  frequency  plays  an  important  role. 
Physically  this  frequency  represents  the  resonance  between  the  junction 
capacitance  Cj  and  the  Josephson  inductance  which  is  derivable  from  Eqs. 

(1)  and  (2): 


Lj  *♦„/!„  •  CIO) 

In  the  context  of  the  microscopic  theory  one  should  note  first  of  all  that 
the  analog  of  IQ,  l.e.,  IjU),  is  frequency  dependent  so  the  inductance 
should  also  be  frequency  dependent.  Secondly,  I  («)  has  a  reactive  part 
which  will  affect  the  resonance.  As  will  become  especially  clear  when  we 
discuss  the  theory  of  pulse  detection,  the  plasma  frequency  is  the  natural 
frequency  for  small  oscillations  of  junction  voltage  about  the  zero  voltage 
state.  Consequently  we  use  small  signal  theory  (34)  to  calculate  the 
effective  inductance  In  the  microscopic  model  Including  the  above  effects, 
with  the  result: 
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[L.f?r  • } "  {•  i>o<°>  *  ♦  tiqp(o) .  . 


Cn) 


In  the  limit  as  ai  goes  to  zero  the  I  terms  cancel  and,  using  Eq.  (9), 
t-eff  becomes  Lj  of  £q.  (10).  If  is  evaluated  at  0  K  It  Is  found 
that  the  singularity  In  Re  I j (u» )  at  the  energy  gap  is  canceled  by  a 
singularity  in  Re  I^Cm)  at  the  same  frequency.  Increases  monotonlcally 

from  Its  minimum  value  lj  at  u  *  0,  reaches  a  value  of  1.57  Lj  at  the 
energy  gap  frequency,  and  asymptotically  approaches  a  linear  dependence 
on  id,  Leff  3  L  j)W&.  The  increase  of  inductance  with  frequency  simply 
means  that  the  additional  parallel  resistive  and  capacitive  current  paths 

dominate  the  circuit  at  high  frequencies. 

Including  high  frequency  effects  the  plasma  frequency  is  given  by 
-1  /2 

(LgfrfC)"  '  .  Therefore  as  C  is  reduced  so  as  to  make  approach  the 
energy  gap  frequency  ug,  one  sees  that  Leff,  which  Is  always  >_  Lj,  makes 
Up  smaller  than  expected  without  high  frequency  corrections.  This  result 
has  Implications  for  the  later  discussion  In  Section  V  where  It  Is  found  that 
the  plasma  frequency  is  directly  related  to  the  time  resolution  of  superconducting 
circuits. 

Returning  to  the  general  subject  of  the  microscopic  theory,  we 
note  that  only  a  few  evaluations  of  It  have  been  made.  These  Include 
detailed  evaluations  of  (a)  I-V  curves  as  a  function  of  R^Cj  (37,39,40), 

(b)  a  proposed  technique  for  producing  picosecond  pulses  (41),  and  (c) 
the  Impulse  response  of  a  Josephson  junction  (42).  Since  the  last 
reveals  the  microscopic  limitations  of  the  pulse  response  of  junctions 
most  clearly,  we  confine  our  attention  to  it  here. 

A  Josephson  junction,  or  any  physical  system  for  that  matter,  does  not 
respond  instantaneously  to  an  applied  stimulus.  Delays  arise  from  the  Inter¬ 
actions  of  electrons  within  the  superconducting  electrodes  of  a  junction. 

These  microscopic  delays  represent  a  fundamental  limitation  to  device  response. 
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To  Illustrate  the  intrinsic  delay  In  a  Josephson  junction,  R.  E.  Harris  (42) 
has  calculated  the  idealized  Impulse  response  of  a  hypothetical  junction  which 
has  no  capacitance.  The  calculation  was  based  on  a  reformulation  of  the 
microscopic  theory  in  the  time  domain  which  Is  discussed  further  in  Appendix 
A.  In  this  calculation  the  stimulus  is  a  voltage  pulse  described  as  a  delta 
function.  The  amplitude  of  the  pulse  is  adjusted  to  produce  a  phase  shift 
of  ir/2  in  $.  The  resulting  current  in  the  junction  is  illustrated  In  Fig.  6. 

The  calculation  reveals  three  significant  aspects  of  the  response.  First 
the  Instantaneous  response  at  time  zero  is  linearly  related  to  the  driving 
voltage  and  is  a  result  of  the  resistive  nature  of  the  I-V  curve  for  large 
bias  voltages.  Second,  after  this  response  comes  a  rapid  rise,  on  a  time 
scale  of  )»/2a  which  Is  about  0.23  ps  for  the  Pb  alloy  junctions  we  will 
discuss  In  the  next  section.  This  rise  reveals  the  Intrinsic  response 
time  of  the  junction.  It  Is  a  very  fast  response,  about  a  factor  of  25 
faster  than  any  observed  to  date.  As  shown  in  Appendix  A  this  delay 
arises  from  both  the  supercurrent  Ij  and  the  quasiparticle  current  I  . 

Third,  at  later  times  the  current  oscillates  at  the  gap  frequency  2A/h. 

This  is  just  the  appearance  in  the  time  domain  of  the  previously-mentioned 
Internal  resonance  named  after  Riedel.  An  important  question  is  now 
apparent:  how  close  can  real  devices,  which  suffer  from  macroscopic 
limitations,  come  to  this  limiting  response  time? 


C.  Macroscopic  Junction  Properties 

Several  intrinsic  macroscopic  parameters  are  needed  for  evaluating 

the  behavior  of  real  Josephson  junctions.  We  will  summarize  the  general 

formulas  for  these  quantities  and  evaluate  them  for  Pb  alloy  technology. 

Following  that  we  return  to  a  discussion  of  the  macroscopic  time  scale 

suggested  by  the  differential  Eq.  (7)  of  the  RCSJ  model. 

The  macroscopic  properties  of  interest  are  the  response  time  RnCj, 

the  plasma  frequency  <op,  the  critical  current  density  Jc,  and  the  Josephson 

penetration  depth  Xj,  all  of  which  are  interrelated.  These  quantities  are 

independent  of  junction  area  except  for  the  plasma  frequency,  which  is 

1/2 

proportional  to  (I „/C.)  '  ,  according  to  Eq.  (8).  This  ratio  is  independent 
0  u 

of  area  if  the  current  density  is  uniform.  Another  point  of  interest,  which 
we  discuss  in  due  course,  is  that  the  plasma  frequency  is  maximized,  for 
a  given  capacitance,  if  the  current  density  is  uniform.  The  Josephson 
penetration  depth  is  a  measure  of  the  maximum  si^e  a  junction  can  have  and 
still  retain  uniform  current  density  (43).  Consequently  we  use  Eq.  (8) 
combined  with  Eq.  (12)  below  to  compute  the  plasma  frequency  with  the  priviso 
that  the  junction  size  not  exceed  Xj  by  Xj. 

Basavalah,  Eldrldge,  and  Matlsoo  (5)  empirically  related  critical 
current  density  to  tunneling  barrier  thickness  d  using  the  following 
relationship: 

jc  -4exP  (’c2d)‘  02) 

where  c^  and  Cg  are  constants  to  be  determined. 

The  next  relationship  of  Interest  Is  between  j  and  the  junction  normal  state 
resistance  «N.  For  a  junction  using  electrodes  of  Identical  superconductors  at 
zero  temperature  we  have 
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where  A  Is  the  junction  area.  This  formula  Is  similar  to  Eq.  (9)  but 
Includes  the  strong  coupling  (29,  30)  correction  Cj.  If  T  >  0,  c3  also 
includes  a  minor  temperature  correction.  For  dissimilar  electrode  materials 
we  use  an  approximate  expression  obtained  by  replacing  2a  with  a-j  +  a2>  rather 
than  the  exact  but  more  complicated  result  of  Ambegaokar  and  Baratoff  (44). 

Combining  Eq.  (13)  with  the  usual  expression  for  capacitance,  C  *  ep  eQA/d, 
the  expression  for  R^Cj  is  obtained: 


Wj 
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The  Josephson  penetration  depth  is  given  by: 


To 
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1+X2 


(14) 


(15) 


where  uQ  Is  the  magnetic  permeability  of  vacuum,  and  and  x2  are  the  magnetic 
penetration  depths  of  the  base  and  upper  superconducting  electrodes  of  the 
Josephson  junction. 

For  evaluation  of  the  foregoing  formulas,  the  constants  c-| ,  c2,  c^, 

A,  ep,  X-j  and  x2  must  be  determined.  In  general  one  would  expect  these 
constants  to  depend  on  the  alloy  composition  of  the  base  and  upper 
electrodes  of  the  junction  and  the  details  of  the  procedure  for  producing 
the  tunneling  barrier.  Recent  work  of  Greiner,  Klrcher,  and  Ames  (45), 
and  of  Baker,  Klrcher,  and  Matthews  (46)  has  begun  to  clarify  what 
factors  are  Important  in  the  fabrication  procedure.  Unfortunately  not  all 
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of  the  material  constants  we  need  have  been  determined  for  the  alloys  of 
greatest  interest  (45)  so  we  are  forced  to  assume  typical  values  taken  from 
several  different  sources. 

We  begin  with  the  constants  in  the  relationship  between  j  and  d  of 
Eq.  (12).  For  this  we  use  tunneling  barrier  data  from  the  work  of 
Eldridge  and  Matisoo  (47)  which  used  Pb-In  alloy  for  the  base  electrode 
and  Pb  for  the  upper  electrode.  In  this  work  the  barrier  is  described 
in  terms  of  various  fractional  concentrations  of  In  and  Pb  in  the  oxide 
and  we  somewhat  arbitrarily  choose  the  data  for  6.5%  In  and  93.5%  Pb. 
Equation  (12)  has  been  fitted  by  a  least  squares  procedure  to  this 

5 

data  to  yield  c^  *  2.207  x  10  ,  »  1.57872  for  d  and  jc  in  units 

2 

of  nm  and  A/cm  ,  respectively.  With  these  constants  Eq.  (12)  gives  jc 

2 

to  within  5%  of  the  experimental  data  (40)  over  the  range  from  1  A/cm 
to  1000  A/cm2. 

The  dielectric  constant  of  the  tunneling  barrier  can  be  evaluated 

from  the  recent  data  of  Gheewala  (48).  He  reports  a  critical  current  of 

50  iiA  for  a  junction  with  2.5  pm  diameter  and  0.2  pF  capacitance.  Using 

2 

the  corresponding  j  of  1019  A/cm  ,  the  barrier  thickness  Is  evaluated 

V» 

from  Eq.  (12),  and  the  dielectric  constant  er  *  12.7  comes  immediately 
from  the  standard  capacitance  formula.  (Note  that  this  value  of  er  is 
obtained  by  combining  data  for  two  different  alloy  systems  and  consequently 
has  considerable,  but  unknown,  uncertainty  associated  with  it.)  A  working 
rule  of  thumb  for  capacitance  is  4  pF  per  raA  of  critical  current  at  1000 
A/ cm2. 

An  I-V  curve  for  the  most  recent  alloy  system  of  Greiner,  et  al .  (45 ) 
allows  the  evaluation  of  additional  constants.  From  the  voltage  of  the 
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steep  rise  of  the  quasiparticle  current,  our  own  data  gives  2a/e  *  2.9  mV. 
For  a  critical  current  of  0.177  mA  we  measure  a  normal  state  resistance 
Rn  *  8.0  n,  which  from  Eq.  (14)  gives  a  value  of  c^  =  0.62.  Values  of 
the  magnetic  penetration  depths  are  taken  from  ref.  (45):  *  137  nm, 

Xg  *  202  nm. 

These  constants  allow  computation  of  the  graphs  shown  in  Fig.  7. 

It  should  be  understood  that  these  graphs  are  only  approximate  since 
they  are  based  on  data  from  more  than  one  alloy  system.  Nevertheless  we 
take  them  to  be  representative  of  the  properties  of  Pb  alloy  junctions 
and  find  them  to  be  useful  for  design  analysis. 

Most  of  the  detailed  analysis  in  subsequent  sections  Involves 
junction  behavior  at  voltages  below  the  steep  rise  of  the  quasiparticle 
current  at  the  energy  gap.  For  this  case  the  appropriate  RC  time  is  not 
Rj^Cj  but  rather  R1 C j ,  where  R-j  Is  the  resistance  below  the  energy  gap  of 
Fig.  4  (b).  Thus  to  obtain  the  physical  RC  time  the  graphical  data  for 
RnCj  of  Fig.  7  must  be  multiplied  by  R-j *  a  number  whl’cl1  depends  on  the 
fabrication  process.  For  example.  In  our  laboratory  R1/RN  was  typically  5 
when  a  PbAu  upper  electrode  was  used  but  a  PbBI  upper  electrode  at  the 
present  time  gives  a  representative  value  of  17. 

Using  data  from  Fig.  7  one  can  readily  establish  that  the  quantity 

cu  Rt  C ,  >  >  1  for  the  entire  range  of  the  figure.  That  means  that  Junctions 
p  i  j 

operating  in  the  subgap  resistance  region  have  plasma  oscillations  that  are 
only  lightly  damped,  a  fact  that  has  significant  consequences  in  the  study 
of  pulse  coincidences  in  a  detector  In  Section  V. 

He  observe  from  Fig.  7  that  R^Cj  and  refer  to  substantially  dif¬ 
ferent  numerical  time  scales.  Consequently  we  must  consider  which  is 
physically  more  Important. 


In  the  RCSJ  model.  Fig.  4  (a),  there  are  three  parallel  current  paths, 

namely  those  through  ,  Cj,  and  the  Josephson  Inductance  Lj  given  by 

Eq.  (10).  For  effective  high  speed  action  of  a  Josephson  junction  In  a 

switching  circuit  Lj  must  not  be  heavily  shunted  by  R-j  or  Cj.  The 

plasma  frequency,  by  definition,  is  the  frequency  above  which  the  Impedance 

of  Cj  is  less  than  that  of  Lj.  Using  the  graphical  data  of  Fig.  7  It 

can  readily  be  determined  that  at  the  plasma  frequency  the  capacitance 

heavily  shunts  R-j .  Thus  R^  plays  only  a  small  role  In  the  high  speed 

dynamics.  From  these  considerations  It  becomes  clear  that  w  and  not 

RnCj  should  be  the  basis  for  the  physically  more  important  time  scale  for 

Josephson  junctions  without  external  shunts  of  low  Impedance. 

Earlier  we  found  the  microscopic  response  time  of  a  Pb  alloy  junction 

to  be  J1/2A  *  0.23  ps.  The  corresponding  limiting  macroscopic  time  Is 

approximately  rp  =  2ir/o»p.  In  order  to  have  the  macroscopic  and  microscopic 

response  times  approximately  equal,  we  find  that  the  plasma  frequency 

must  be  approximately  equal  to  the  energy  gap  frequency  2  2a/J4 .  In  this 

regime  the  RCSJ  model  has  Inaccuracies  of  about  30S  which  we  will  largely 

Ignore.  More  importantly  we  note  that  the  current  density  required  for 

5  2 

such  a  high  u  would  be  well  In  excess  of  10  A/cm  for  Pb  alloy  technology. 

P 

Such  high  Jc  means  that  nonequilibrium  and  self  heating  effects  would  be 

5  2 

significant.  Even  though  current  densities  In  excess  of  10  A/cra  have 
been  achieved  (49),  little  Is  known  about  the  dynamical  properties  of 
such  junctions.  Consequently,  It  must  be  left  as  an  open  question 
whether  or  not  the  macroscopic  and  microscopic  speeds  can  be  made  comparable. 

An  example  of  what  Is.  achleveable  comes  from  the  work  of  Broom, 

Jutzl,  and  Mohr  who  used  junctions  with  current  densities  of  28000 
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2 

A/cm  in  a  memory  cell  (.501.  Heating  effects  were  sufficiently  large  in 

this  case  that  the  I-V  curve  was  traced  by  a  pulsed  method  to  reduce 

errors  from  self  heating  of  the  junction.  Talcing  this  example  as  a 

rough  estimate  of  the  maximum  j  that  is  practical  without  using  extra 

thick  electrodes  or  special  substrates  (51),  Fig.  7  gives  the  corresponding 
+1 2  -1 

(Dp  *  3  x  10  s  .  This  is  a  factor  of  ten  below  the  plasma  frequency 
needed  to  make  the  macroscopic  and  microscopic  response  times  equal. 


IV.  Switching  Time 


The  switching  time  of  a  junction  from  the  zero  voltage  state  to  the 
energy  gap  voltage  divides  roughly  Into  two  parts,  a  turn-on  delay  time 
during  which  the  voltage  remains  small,  and  a  rise  time  during  which  the 
voltage  is  rapidly  Increasing,  as  Illustrated  in  Fig.  9.  The  analysis 
of  these  given  here  Is  based  on  the  RCSJ  model.  For  the  period  during 
which  the  voltage  remains  below  the  gap  voltage,  one  may  use  Eq.  (6)  with 
R.  replaced  by  R  where  R  Is  the  parallel  combination  of  the  subgap  resistance 
R.|  and  a  possible  load  resistance  R^.  Setting  I(t)  *  Ijj(t)  (6)» 

we  note  that  if  I^Ct )  changes  sufficiently  slowly  then  Eq.  (6)  can  be 
written  in  the  form  of  an  energy  equation 

fa  (T  +  V)  a  '  fej  *  °6) 

where  the  kinetic  and  potential  energies,  T  and  V,  are  respectively 
T  *  $  and 


p  *  -  dip  (cos$  +  (I{/Iq)  ♦  !  •  ^7) 

These  equations  suggest  the  helpful  analogy  of  a  particle  of  unit  mass 

subject  to  a  viscous  retarding  force,  sliding  down  a  potential  energy 

"washboard*.  Figure  8  Illustrates  the  potential  energy  as  a  function 

of  4  and  Ib.  Clearly  there  Is  no  stable  position  for  the  particle  if 

I  >  I  .  Consideration  of  the  washboard  picture  helps  one  to  understand 
b  o 

qualitatively  many  features  of  the  RCSJ  model,  including  the  switching 
behavior  of  a  junction.  For  example,  the  particle  initially  moves  slowly 


until  the  downward  curvature  of  the  potential  is  reached.  This  corresponds 
to  the  turn-on  delay  period,  which  was  first  discussed  by  E.  P.  Harris  (52). 
Subsequently  the  particle  accelerates  until  either  its  terminal  velocity 
is  reached,  corresponding  to  a  voltage  IbR,  or  the  velocity  corresponding 
to  the  gap  voltage  is  reached.  This  defines  the  rise  time  of  the  switching 
process.  Figure  9  shows  a  simulation  which  illustrates  this  behavior. 

We  will  develop  formulas  for  the  turn-on  delay  and  the  rise  time  tR 
for  lightly  and  moderately  loaded  circuits.  The  sum  of  these  is  the 
approximate  switching  time. 

A.  Turn-on  Delay  of  a  Josephson  Junction 

Suppose  that  the  current  Ib(t)  is  increased  slowly  from  zero  so  that 
the  junction  phase  $  will  be  about  ir/2  when  Ib  reaches  IQ.  At  that  point 
Ib  is  then  abruptly  stepped  to  IQ(1  +  5Q)  where  5Q  is  the  fractional  over¬ 
drive.  This  is  a  bit  artificial  in  that  in  practice  I&  typically  ramps 
up  through  I  at  a  finite  rate.  One  technique  which  approximates  this 
case  is  that  In  which  a  junction  Is  biased  at  a  constant  current  below 
IQ,  and  then  triggered  by  a  current  pulse  resulting  from  the  switching 
of  another  junction,  as  discussed  In  Section  II. B.  We  will  shortly 
generalize  to  the  case  where  the  initial  Ib  is  below  IQ.  The  effects  of 
a  finite  ramp  rate  will,  however,  not  be  taken  up  here. 

In  Eq.  (6)  we  thus  use  Ib  ®  IQ  (1  +  $Q)  for  t  >  0,  and  have 
<j>(0)*ir/2.  We  set  sin$  equal  to  unity  and  consider  the  evolution 
of  V  until  $  advances  by  one-half  radian.  The  supercurrent  term  is  in 
error  by  just  12%  at  the  end  of  this  period,  so  that  the  estimated  time 
for  4  to  get  to  ir/2  +  1/2  Is  quite  accurate.  Fortunately,  $  *  ir/2  +  1/2 
Is  also  about  the  place  where  $  begins  Increasing  significantly,  as 
simulations  show.  We  thus  take  this  time  as  an  estimation  to  tQ.  From 
Eq.  (6)  we  then  easily  obtain,  provided  tq  «  RCj, 
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where  Is  given  In  Eq.  (8). 

E.  P.  Harris  (52)  has  published  an  expression  for  the  turn-on 
delay  time  for  this  case.  He  approximated  sin<j>  by  the  first  two  terms  in  Its 
Taylor  expansion,  Integrated  the  resulting  expression  exactly  in  terms  of 
elliptic  Integrals,  and  approximated  the  latter  by  simpler  functions.  From 
this  result  he  calculated  the  time  required  for  the  phase  to  advance  to  3w/2. 
He  subtracted  from  this  the  time  required  for  the  phase  to  advance  the  same 
amount  if  only  the  junction  capacitance  were  being  charged.  His  result, 
again  for  tq  <  <  RCj,  is 


Figure  10  compares  Eqs.  (18)  and  (19).  Also  shown  are  values  of  Tp 
taken  from  numerical  simulations.  The  latter  have  been  determined  by 
extrapolating  the  rise  time  portion  of  the  curve  down  to  the  time  axis. 
Typically  there  Is  a  few-percent  uncertainty  In  this,  both  because  of 
the  superimposed  Josephson  oscillations,  and  because  the  rising  portion 
Is  not  a  straight  line  even  with  the  oscillations  removed.  The  points 
shown  In  Fig.  10  were  determined  from  the  midpoint  of  two  straight  lines 
drawn  tangent  to  the  first  two  maxima  and  minima,  as  shown  In  Fig.  9. 

As  is  seen  in  Fig.  10,  Eq.  (18)  Is  more  accurate  than  Eq.  (19)  In 
the  region  of  about  1 -50%  overdrive.  Outside  this  region,  Harris’ 
expression  Is  more  accurate.  Thus,  from  1 -50*  overdrive,  Eq.  (18)  would 


be  preferred  for  reasons  of  both  accuracy  and  simplicity.  Below  1%, 
turn-on  delay  becomes  a  more  serious  problem,  and  one  may  prefer  to 
use  Eq.  (19);  the  simple  formula  overestimates  tg  there.  Above  about 
50%,  Tg  is  typically  very  small  compared  to  the  risetirae  and  therefore 
of  decreased  significance. 

Now  we  remove  the  restriction  that  I.  is  initally  at  I  .  Let 
Ib(t  <  0)  =  1^  <_  IQ,  and  $(0)  =  *  sln-^  (Ifa.j/I0) •  There  are  two  new 

effects  to  consider  now.  One  is  that  the  particle  in  the  washboard 
model  must  travel  a  longer  distance  if  the  same  "final”  position  is  used 
to  define  the  delay  time.  The  other  is  that  the  initial  potential 
energy  is  now  greater,  giving  the  particle  more  velocity  after  the  step, 
thus  decreasing  the  delay  time.  This  Is  illustrated  in  Fig.  8  (b).  The 
latter  effect  usually  dominates. 

Me  write  1^  *  IQ  (1  -  5^)  and  Ib  (t  >  0)  *  lQ  (1  +  SQ]  where  «^  >  0 
and  SQ  >  0.  Our  procedure  is,  in  effect,  to  approximate  s1n<j>  by  its 
Initial  value  and  to  drop  the  damping  term  V/R  from  Eq.  (16),  which  thus 
gives  a  constant  acceleration  to  the  system  particle  and  an  easy  determination 
of  Tg.  Lest  this  seem  too  crude,  we  shall  show  that  the  same  result  is 
obtained  by  expanding  s1n$  through  the  term  linear  term  In  <|>  ,  followed 

by  some  usually  excellent  approximations  and  the  assumption  that  tg  «  RCj 
as  In  the  preceding  case.  When  sin*  is  thus  expanded,  the  potential  energy  of 
Eq.  (17)  is  changed  to  a  typically  shallow  parabola.  Equation  (6)  then 
describes  a  damped  harmonic  oscillator: 

I  (1«  )  -  Cj*  +  V/R  +  IjjSln^+C*  -  ♦i)I0cos*i.  (20) 

This  has  the  solution 


25 


(21) 


♦(t>0)  »  +  ,  +  ^l^-1  |l-e't/2RCJ(cosait  +  (ZuRCjl^sIn  wt)j, 

where 

0)  *  {4R2C^i>2  cos^  -1  )1/2  /2RCj  *  ^(cos*,)1'2.  (22) 


Thus  &»  Is  essentially  the  phase-dependent  plasma  frequency.  The  approximate 

equality  in  Eq.  (22)  Is  quite  typically  valid:  If  we  take  Vg/ IQ  as  a 

minimum  value  for  R  (to  ensure  switching  to  the  gap),  the  condition  for 

_  2 

the  validity  of  the  approximate  equality  is  cos$.j>>4>0I0/4CjVg.  For  the  junction 

of  Fig.  9  this  becomes  cos^  >>  10"2,  or  6^  >  >  10*®.  Clearly,  R  may  be 

considerably  less  than  Vg/IQ,  as  in  high  speed  latching  logic,  without 

practically  invalidating  the  approximation  of  Eq.  (22). 

The  coefficient  of  sin  ut  in  Eq.  (21)  is  very  small  under  the  same 

condition.  Thus  the  Initial  motion  of  4  results  principally  from  the 

2 

cos  nit  term.  Expanding  cos  wt  ■  1  -  (ait)  /2,  and  defining  td  as  the 
time  at  which  $  reaches  w/2  +  1/2,  as  before,  we  find 


p  +  2(tt/2 
l  fio+6i 


(23) 


This  is  a  simple  generalization  of  Eq.  (18)  to  which  it  reduces  as  1^ 
goes  to  IQ.  Unless  ^  is  far  from  ir/2,  the  numerator  is  not  very  sensitive 
to  the  precise  value  of  but  the  denominator  can  be  very  sensitive  to  6^ 
confirming  our  earlier  remarks  concerning  the  relative  effects  of  the 
initial  position  and  potential  energy. 

Equation  (23)  underestimates  the  delay  time  for  small  overdrives,  as 
determined  from  simulations.  The  reason  is  that  the  approximations 


used  Tn  the  above  analysis  are  equivalent  to  taking  the  acceleration  of 

♦  as  unchan9ed  from  its  Initial  value  during  the  motion.  For  «f  »  0.1,  the 

result  of  Eq.  (23)  for  the  junction  of  Fig.  9  is  about  30%  low  for  6  » 

o 

0.01,  correct  for  <5q  about  0.2,  and  progressively  too  high  beyond  that. 

It  is  not  difficult  to  refine  the  approximation  at  the  expense  of  a  more 
complicated  formula.  However,  since  different  methods  for  defining  r0 
can  give  values  differing  by  several  percent,  and  since  one  is  ordinarily 
only  interested  in  estimating  tq,  the  estimate  of  Eq.  (23)  is  surely 
adequate  for  most  purposes. 

Thus,  in  summary,  the  total  size  of  the  step  rather  than  just  the 
amount  of  overdrive  is  the  important  quantity  in  the  turn-on  delay  resulting 
from  a  step  increase  in  1^.  Further,  is  modified  by  the  initial  phase. 
Discussions  of  the  turn-on  delays  in  single- junction  and  two- junction 


SQUIDs  are  given  in  Appendices  B  and  C. 


The  rise  time  tr  refers  to  the  time  duration  of  the  rapidly  rising 
portion  of  the  switching  transition  Illustrated  In  Fig.  9.  We  shall  show 
that  the  rise  time  of  the  RCSJ  model  Is  approximately  R^Cj 'where  R^  Is 
the  normal  state  resistance  of  the  junction.  Although  this  result  Is 
wtdely  known.  It  Is  often  misquoted.  It  needs  to  be  emphasized  that 
the  resistance  here,  R^,  Is  not  the  actual  resistance  R  shunting  the 
Junction. 

During  the  steep  rise  of  the  transition  the  sin#  term  of  Eq.  (6) 

Is  rapidly  oscillating  so  that  the  current  I0s1n<(i  averages  nearly  to 
zero.  Equation  (6)  with  the  s1n$  term  neglected  describes  the  charging 
of  the  Junction  capacitance.  The  solution  to  this  equation  is 

t 

‘  V(t)  -  IjjRfl-expC-t/RCjl] , 

tf  the  bias  current  has  the  value  Ijj  which  does  not  change  significantly 
during  the  rapid  rise.  We  choose  t  »  0  as  the  time  when  the  significant 
acceleration  begins.  If  IbR  Is  less  than  the  gap  voltage  Vg,  the  voltage 
rises  only  to  I^R  with  time  constant  RCj.  However,  if  Vg  Is  less  than 
IbR,  the  rtsetlme  xR  at  which  the  gap  is  reached  Is 


•  -RCjfnCl-Y  /IbR). 


C25) 


ff,  moreover,  Vg  <  <  I^R  we»tee..that 


I. 


CA  4 

■  ?Vj  \ 


where  we  have  used  the  relation  V-  ■  ^rom  ^us 


or 
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,f  h  -  'o  and  vg  «  *0  *.  10- 90S  rise tine  .ill  be  close  to  «Bc,. 

This  mey  be  verified  In  Fig.  9.  For  strong  coupling  superconductors, 

the  10-902  risetime  will  be  somewhat  longer  than  R  C 

N  J" 


V.  Pulse  Generation,  Detection,  and  Measurement 

In  this  section  we  adopt  the  primitive  view  that  pulses  have  only  two 
characteristics  of  Interest,  a  height  and  a  width.  As  a  reasonable  ap¬ 
proximation  to  practical  pulses  we  use  the  Gaussian  pulse  shape  for  most 
of  the  analysis.  The  main  problem  of  Interest  Is  to  develop  the  concepts 
needed  for  measurements  of  pulses  that  are  available  on  a  repeating  basis. 

A.  Pulse  Generators 

% 

Probably  the  simplest  means  of  producing  pulses  In  superconducting 
electronics  Is  the  circuit  of  Fig.  2  (a).  The  junction  J1  is  assumed  to 
be  Initially  In  the  zero  voltage  state  with  I$(t)  +  1^  <  IQ.  The  switching 
transition  of  the  hysteretic  Junction  Jl  occurs  very  shortly  after  I$(t) 
exceeds  IQ  -  1^.  The  resulting  Increased  impedance  of  Jl  diverts  the  bias 
current  down  the  transmission  line  producing  the  rising  edge  of  the  pulse. 

If  the  shunt  Impedance  across  the  Junction  Is  low  enough  (53)  the  junction 
will  not  stay  In  the  switched  position,  l.e..  It  will  not  latch,  but  will 
go  back  to  the  zero  voltage  state  producing  the  trailing  edge  of  the  pulse. 

A  related  technique  for  producing  especially  narrow  pulses  Is  to  pass 
the  rising  edge  of  the  switching  transition  through  a  differentiating  network, 
as  Indicated  for  example  as  part  of  Fig.  3. 


Paris  (28)  and  Tuckerman  (29)  have  recently  demonstrated  a  pulse 
generator,  based  on  a  SQUID  Interacting  with  another  junction,  which 
produced  12  ps  wide  pulses. 

Several  years  ago  Peterson  and  McDonald  (54)  proposed  a  conceptually 
different  form  of  pulse  generator  based  on  a  Josephson  junction  driven  only 
by  a  microwave  current  source  and  loaded  by  a  shunt  resistor.  In  this 
circuit  the  junction  may  be  either  hysteretic  or  nonhysteretlc  without 
changing  the  qualitative  behavior  of  the  circuit.  A  train  of  alternating 
positive  and  negative  pulses  is  produced  when  the  drive  current  amplitude 
just  exceeds  IQ  and  encounters  a  strong  nonlinearity.  One  cycle  of  the 
pulse  train  Is  shown  In  Fig.  11.  This  design  is  notable  in  part  because 
it  has  been  analyzed  using  the  microscopic  theory  of  Junctions  (35)  to 
calculate  the  resulting  pulse  widths.  A  comparison  of  the  waveforms 
calculated  using  the  microscopic  theory  and  the  RCSJ  model  is  one  of  the 

few  benchmark  comparisons  between  the  two  theories.  The  differences 

/ 

between  the  waveforms  are  approximately  1 5%  for  pulse  widths  of  4.7  ps  and 
about  35%  for  the  narrowest  pulse  widths  of  0.8  ps. 

The  original  analysis  did  not  consider  the  details  of  pulse 
detection  and  measurement.  The  need  to  understand  these  processes 
provided  the  primary  motivation  for  the  analysis  which  follows.  In  all 
of  the  simulations  which  follow,  the  affects  that  would  arise  In  practice 
from  parasitic  Inductance  and  finite  source  and  load  Impedances  are 
neglected.  The  goal  is  to  describe  the  main  qualitative  features 
of  pulse  Interactions  with  detectors. 


B.  Pulse  Detectors 


The  elementary  problem  considered  here  is  represented  in  Fig.  13  (a),  in 
which  a  current  source  Ip(t)  provides  a  single  pulse  to  a  Josephson 
junction  detector,  and  a  current  source  Ip  provides  dc  bias.  We  then 
ask  under  what  circumstances  does  the  junction  voltage  V  switch  from 
zero  to  the  energy  gap  voltage.  The  basis  for  our  analysis  is  Eq.  (6) 
with  the  driving  current  I{t)  *  Ip(t)  +  Ip. 

If  the  pulse  risetime  is  slow  enough,  the  detector  will  switch 
as  soon  as  Ip{t)  +  Ip  «  IQ,  assuming 


Jb  41  V  However*  1f  th*  Pu1se  risetime  Is  too  fast  or  If  the  pulse 
•width  Is  too  narrow,  large  deviations  from  this  limiting  case  are  sug¬ 
gested  by  the  qualitative  Ideas  of  the  washboard  model.  The  emphasis 
In  the  present  discussion  will  be  on  determining  the  minimum  value  of  1^ 
for  which  switching  occurs,  i.e.,  the  threshold  condition,  with  a  particular 
pulse  amplitude  and  pulse  width.  In  a  subsequent  section  the  measurement 
of  pulse  width  Is  explored.. 

Initially  we  consider  simulations  In  which  the  applied  pulse  has  a 
Gaussian  shape  described  as: 


ytj » i 


pi 


exp  {-(4£/i2) 


9 


(27) 


where  1 p(t)  a  Ip(t)/I0,  1pl  =  Ipl/10,  *s  the  of  arrival  of  the 
center  of  the  pulse  at  the  detector,  and  W  Is  the  full  width  of  the 
pulse  at  a  current  level  of  j  ip1. 

Figure  1 2  Is  a  composite  Illustration  of  the  calculated  responses  of 
a  detector  junction  to  pulses  of  four  different  widths.  The  driving 
pulses  are  shown  In  the  first  row  of  illustrations.  The  junction  properties 
and  the  pulse  amplitudes,  1p-|  ■  0.5,  are  the  same  for  all  cases.  Various 
junction  responses  are  shown  for  biases  1^  just  below  and  just  above 
threshold  for  switching.  Although  both  bias  levels  are  very  close  to 
threshold,  the  corresponding  graphs  qualitatively  represent  the  general 
behavior  for  any  values  of  1b  below  and  above  threshold,  respectively. 

In  the  second  row  of  Illustrations  It  Is  seen  that  for  ip  below 
threshold  the  voltage  stays  near  zero  and  undergoes  damped  plasma 
oscillations  after  excitation  by  the  applied  current  pulse.  For  1b 
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above  threshold  the  voltage  rises  rapidly  from  zero  to  the  energy  gap 
value  Vg.  Josephson  oscillations  driven  by  are  evident  after  the 
junctions  switch  in  all  cases,  but  the  time  scale  Is  too  coarse  for 
their  resolution  In  the  two  left-most  illustrations  of  this  group. 

These  figures  clearly  Illustrate  why  a  Josephson  junction  is  useful  as  a 
pulse  detector:  The  device  permanently  changes  its  voltage  state  after 
the  passage  of  a  brief  transient  event  if  the  pulse  height,  pulse  width, 
and  bias  have  the  proper  relationship.  Of  course  this  Is  much  the  same 
as  the  operation  of  a  Josephson  latching  logic  gate. 

In  the  third  row  of  Fig.  12  the  time  dependence  of  the  phase  $  of 
the  junction  is  shown  for  each  of  the  four  driving  pulse  widths.  For 
the  curves  corresponding  to  bias  values  below  threshold, the  phase  is 
driven  to  some  maximum  positive  value  by  the  applied  pulse  and  then 
undergoes  damped  plasma  oscillations  in  returning  to  its  equilibrium 
position.  The  plasma  oscillations,  which  seem  like  an  unimportant 
detail,  assume  major  Importance* when  we  consider  the  detection  of  pulse 
coincidences  In  a  later  section.  An  extensive  numerical  study  of  the 
phase  supports  the  general  idea  that  a  junction  will  switch  If  the  phase 
reaches  a  certain  universal  critical  value.  Thus  If  one  could  watch  the 
evolution  of  the  phase  In  time  using  this  critical  phase  Idea,  one  would 
know  well  In  advance  whether  or  not  the  voltage  was  committed  to  a 
trajectory  to  Vg.  The  critical  phase  Is  ir-<*>b,  where  $b  Is  the  phase  due 
to  the  steady  bias: 


s1n”] 1 b* 


(28) 
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This  idea  is  central  to  the  detailed  models  that  will  be  used  to  analyze 
threshold  switching  conditions.  Zappe  has  used  a  similar  argument  when 
analyzing  response  to  a  current  step  (53).  We  note  that  the  critical  phase 
idea  originates  from  considerations  with  ideal  pulses.  If  a  pulse  is  complex, 
i.e. ,  if  it  goes  alternately  both  positive  and  negative  for  example,  the 
critical  phase  idea  may  not  apply.  This  is  only  to  say  that  the  direction  of 
motion  of  a  pendulum  can  always  be  reversed  regardless  of  its  position  relative 
to  a  critical  angle. 

The  switching  process  for  pulse  detectors  involves  a  turn-on  delay 
and  a  rise  time  which  are  qualitatively  the  same  as  discussed  in  Section 
IV.  These  delays  are  particularly  evident  for  the  switching  transition 
of  Fig.  12  (k),  the  simulation  with  the  narrowest  driving  pulse.  The 
formulas  given  in  Section  IV  for  the  turn-on  delay  are  derived  for  the 
case  in  which  the  applied  current  remains  at  a  constant  level  above  the 
critical  .current,  thus  continually  feeding  energy  into  the  system  and 
shortening  the  turn-on  delay.  Turn-on  delay  expressions  for  the  pulsed 
case  can  be  derived.  However,  they  are  not  simple  and  we  will  not 
present  them  here. 

In  summary  Fig,  12  Illustrates  the  behavior  of  the  detector  biased 
near  threshold  for  a  variety  of  pulse  widths.  The  relationship  between 
pulse  height,  pulse  width,  threshold  bias,  and  the  properties  of  the 
detector  junction  will  now  be  analyzed  by  considering  three  different 
cases:  (1)  Static  Limit,  (2)  Impulse  Approximation,  and  (3}  Square 
Pulse  Approximation. 

1 .  Static  Limit 

This  limit  applies  to  very  broad  pulses,  broad  in  comparison  with 
a  junction  plasma  period  tp  s  2-x/u^.  With  the  pulse  current  rising 

32 


/ 


sufficiently  slowly  It  is  obvious  that  the  junction  will  switch  very 
soon  after  ip(t)  +  1^  exceeds  unity.  The  pendulum  analog  (55)  provides 
a  simple  interpretation  for  this  case.  This  analog  is  based  on  the 
observation  that  the  differential  equation  for  the  angle  of  deflection 
of  a  simple  pendulum  is  similar  to  Eq.  (7)  for  the  phase  of  a  Josephson 
junction  If  one  makes  a  suitable  Identification  of  the  analogous 
variables.  The  torque  ^analogous  to  the 


! 
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driving  currentjis  balanced  by  gravitational  torque  so  that  the  bias 
angle  In  the  analogous  cases  Is: 

♦b*510’1  (tSt)  "  sin‘lfb-  (29) 

where  mg  Is  the  gravitational  force  on  the  pendulum  of  length  i.  Since 
the  torque  required  to  hold  the  pendulum  statically  at  angle  Is 
proportional  to  s1nq>b ,  It  Is  clear  that  the  torque  Increases  for  angles  from 
zero  up  to  ir/2  and  then  decreases  for  larger  angles.  Thus  for  any  torque 
t.  <  mgz,  there  are  two  angular  positions  where  that  torque  can  hold  the 
pendulum:  ^  and  -nr-<j>b .  However  the  upper  position  Is  unstable  and  random 
fluctuations  In  position  will  cause  the  pendulum  to  go  "over  the  top".  Thus 
our  hypothesis  Is  that  If  the  pendulum  can  ever  reach  it  -$b>  the  mirror  Image 
position  of  for  reflection  In  a  horizontal  plane,  then  continuous  rotation 
of  the  pendulum  will  be  Initiated  (assuming  no  damping),  which  Is  analogous 
to  switching  to  the  energy  gap  In  a  Josephson  junction. 

In  the  static  limit  the  V  and  V1  terms  (which  are  equivalent  to 
|  and  $  terms)  In  Eq.  (6)  are  negligible.  An  examination  of  this 
equation  leads  to  the  conclusion  that  In  the  static  limit  the  critical 
angle  Idea  has  a  special  twist.  It  becomes  obvious  that  In  order  for  a 
slow  pulse  to  carry  the  pendulum  to  the  mirror  position  with  the  lowest 
possible  (threshold)  bias  level,  the  pendulum  must  reach  *  »  w/2  In 
coincidence  with  the  pulse  maximum. 

The  simulations  In  the  first  column  of  Fig.  12,  with  a  pulse  width  of 
324  ps,  are  close  to  the  static  limit.  Threshold  bias  Is  between  the  two 
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illustrated  bias  levels  of  0.500214  and  0.500215.  Thus  we  see  immediately 
that  the  threshold  bias  1b(Th]  Is  very  near  the  simple  relation  ib(Th)  + 
ipl  3  s*nce  *P1  3  °-5- 


2.  Impulse  Approximation 

Another  limiting  case  which  can  be  evaluated  in  detail  is  for  extremely 
narrow  pulses,  the  Impulse  approximation.  In  this  case  the  current  very 
rapidly  rises  and  *alls  back  to  1b  before  the  junction  can  respond.  All 

of  the  charge  qp  in  the  impulse  goes  onto  the  capacitor,  charging  it 

1  2 

to  voltage  Vm  and  providing  an  amount  of  energy  ^CjV^  for  driving  the 


system. 


i: 


Ip(t)dt 


C.V_. 

j  m 
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In  the  pendulum  analog  an  impulse  of  torque  is  applied  but  it  Is  over 
before  the  pendulum  has  time  to  move.  In  the  limiting  case  the  pendulum 
undergoes  a  step-function  change  In  angular  velocity.  Based  on  our 
previous  discussion  of  the  analog,  which  ignores  dissipation,  we  see  that 
the  system  will  go  into  a  state  of  steady  rotation  if  the  impulse  causes 
the  pendulum  to  reach  the  mirror  Image  position. 

Another  source  of  energy  is  the  bias  supply.  To  calculate  this  energy, 
it  is  appropriate  to  integrate  from  the  initial  phase  <j> ^  to  the  phase  of  the 
image  position  tt  -  <j>b: 


r^b 

V(t)  dt-  *0Ib  \  d<j> 


♦0  Ib(ir  -  2*b).  (31) 


The  only  other  energy  In  the  system  is  the  analog  of  the  potential 
energy  of  the  pendulum  (42),  which  Is  given  by:  IQcosij>.  As  the  pendulum 

Is  raised  the  potential  energy  continuously  Increases.  The  change  In 
potential  energy  between  and  ir-$b  Is  2JoIo  cos$b. 

Ignoring  dissipation,  conservatic  r  energy  requires  that  the 
energy  in  the  Impulse  plus  the  energy  supplied  by  the  bias  equals  the 
change  in  potential  energy: 


1 

7 


qn2  - 

— Ib(Th){ir-  2$b(Th)} 
J 


2*o  !o  cos*bCTh)*  (32) 


at  threshold. 
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For  a  Gaussian  pulse  described  by  Eq.  (27): 


qp  -  1.064  iplW 


(33) 


Combining  (32)  and  (33)  and  changing  to  reduced  variables  gives: 


1-pl2  "  ^f7  {2  c°s^b(Th)  -  1b(Th)  t*-2+b(Th)]  } 


(34) 


This  is  the  desired  relationship  between  pulse  height  1p^  and  bias  at  threshold 
i‘b(Th)  for  a  narrow  pulse  of  width  W.  A  clearer  idea  of  what  "narrow"  pulse 
means  will  be  obtained  by  comparing  this  formula  with  detailed  simulations. 

Figures  12  (j)  -  (l)  illustrate  an  example  calculation  for  which  the  impul 
approximation  is  reasonably  accurate.  For  this  case  the  pulse  width  W  * 

0.063  Tp  *  0.324  ps,  and  the  pulse  height  ip-|  •  0.5.  Figure  12  (k) 
illustrates  the  time  dependence  of  the  detector  voltage  for  biases  of 
0.95527  and  0.95528,  which  are  just  below  and  just  above  threshold. 

Setting  1 b(Th)  *  0.95528,  the  value  for  ip^  at  threshold  can  be 
calculated  from  the  Impulse  approximation  formula,  Eq.  (34).  The  result 
Is  1pi  *  0.4505,  which  Is  about  10%  off  the  simulation  value  of  0.5. 

This  simulation  differs  from  the  static  limit  In  that  the  sum  1 fa(Th)  + 

1  pi  «  1.455  Instead  of  simply  1.  Greater  variations  of  the  sum  (1b  +  1p1)Th 
are  found  if  1b  Is  held  fixed  and  the  threshold  pulse  height  Is  evaluated. 

The  required  pulse  height  diverges  as  W’^,  according  to  Eq.  (34). 
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A  physically  more  Interesting  difference  from  the  static  limit  is 
illustrated  In  Fig.  12  0 i)  where  the  time  development  of  the  phase  Is 
shown.  For  ib  just  barely  below  threshold,  $  advances  to  a  maximum  of 
1.865  radians  and  then  falls  back  and  undergoes  damped  plasma  oscillations. 

The  maximum  advance  of  the  phase  Is  very  close  to  the  image  position  w  -  $b  ■ 

1.872  radians,  validating  the  physical  significance  of  that  position. 

The  plasma  oscillations  are  much  more  evident  below  threshold  In 
the  impulse  limit  than  in  the  static  limit.  Fig.  12  (c).  This  behavior 
is  readily  understood  from  a  Fourier  analysis  of  the  two  different  input 
pulses.  The  bandwidth  of  the  pulses  (defined  as  the  frequency  at  which 
the  frequency  distribution  has  fallen  to  1/2  its  maximum  value)  is: 


ZLni 

fl/2  *  Tw“  ‘ 


(35) 


In  approximating  the  static  limit  we  had  W  *  324  ps,  implying  f^2  » 

1.36  GHz.  The  pulse  used  In  the  Impulse  limit  has  VI  *  0.324  ps,  and 
f^2  *  1360  GHz.  For  both  of  these  simulations,  which  differ  only  in 
the  applied  pulse  width,  the  plasma  frequency  is  194  GHz.  Comparing  the 
two  values  of  f-j^  the  P^sma  frequency  It  Is  obvious  that  the  wide 
pulse  has  very  little  energy  available  for  exciting  plasma  oscillations, 
whereas,  in  the  narrow  pulse  case  the  plasma  frequency  Is  well  inside 
the  bandwidth  of  the  pulse  so  substantial  energy  is  available.  Excitation 
of  plasma  oscillations  takes  on  much  greater  importance  when  we  deal 
with  pulse  width  measurements  In  the  next  section,  so  It  Is  Important  to 
understand  what  produces  them. 


37 


3.  Square  Pulse  Approximation 


The  static  limit  and  the  impulse  approximation  leave  the  intermediate 

cases  unevaluated  on  a  quantitative  basis.  To  obtain  workable  formulas  for 

these  cases  it  is  necessary  to  relinquish  the  Gaussian  description  of  the  pulse 

in  favor  of  a  simpler  pulse  shape.  Dhong  and  Van  Duzer  (56)  have  shown  that 

analytical  results  for  threshold  conditions  with  arbitrary  pulse  widths  can 

be  obtained  if  the  driving  pulse  has  vertical  sides  and  square  corners.  For 

describing  this  pulse  we  use  the  same  notation  as  for  the  Gaussian  pulse,  i.e., 

it  has  a  maximum  current,  normalized  to  I  .  of  1  ,  and  a  width  W.  In  this  cal- 

o  pi 

culation  the  system  is  Initially  assumed  to  be  dissipationless  but  a  correction 
is  made  for  that  later.  Using  the  washboard  model,  described  earlier,  and 
energy  conservation,  the  phase  which  must  be  reached  by  the  end  of 

the  square  pulse  in  order  to  have  continuing  movement  down  the  washboard  is: 


_.nib  +  (21b  +  1pi)*b  +  2cos*b 

,  rf 


(36) 


Using  a  linearized  s1n$  term,  expanded  about  $b,  the  differential  equation 
(6),  with  I(t)  ■  Ip(t)  +  Ij,»  can  be  solved  for  the  time  dependence  of 
the  phase  during  the  pulse: 


(  37) 


where  w'  5  u  /cos$b.  («p  is  a  general izatlon  of  the  plasma  frequency  up  explicitly 
includes  the  phase  dependence  of  the  effective  Inductance  of  a  Josephson  junction). 
Solving  Eq.  (37)  for  t  and  then  setting  t  ■  W  so  that  $w  yields: 


Combining  the  threshold  relationship  Eq.  (36)  with  Eq.  (38)  allows  the 
calculation  of  the  pulse  width  W  at  threshold  for  a  given  bias  and  pulse 
height,  assuming  no  damping. 

Including  a  correction  for  damping,  Dhong  and  Van  Duzer  $fi«)  obtained 
the  following  expression  for  the  pulse  width  at  threshold: 


W'  ■  W  + 


W-VsinV  ■ 


(39) 


This  relationship  has  been  evaluated  for  a  number  of  values  of  pulse  width 
with  fixed  pulse  height  1^  *  0.5.  The  results  are  displayed  in  Fig.  13 
along  with  several  data  points  from  an  evaluation  of  the  impulse  approximation 
formula,  Eq.  (34).  In  the  figure  these  points  are  to  be  compared  with  the 
results  of  circuit  simulations,  represented  by  the  solid’llne. 

The  square  pulse  approximation  predicts  the  threshold  current  to  an 
accuracy  of  102  or  better,  as  compared  with  the  simulations,  over  the  pulse 
width  range  0.09  xg  <  W  <  0.5  xg.  This  approximation  diverges  from  the 
simulations  for  small  pulse  widths  because  the  correction  term  for  dis¬ 
sipation  becomes  comparable  to  W.  For  large  pulse  widths  It  diverges  from 
the  simulations  because  the  linearization  of  s1n<|>  becomes  inaccurate.  The 
Impulse  approximation  Is  comparably  accurate  for  W  <  .25  x^.  One  caution 
that  should  be  observed  Is  that  the  RCSJ  model  Is  not  expected  to  be  very 
accurate  for  pulse  widths  of  the  order  of  1  ps  or  less,  the  regime  where  the 
microscopic  frequency  dependence  of  the  superconducting  response  becomes 
significant.  In  a  related  problem  (41 )  for  which  comparisons  with  the 
microscopic  model  were  made  for  an  0.8  ps  pulse  width,  the  differences 
between  the  models  were  about  352. 


The  most  Important  general  conclusion  from  Fig.  13  Is  the  observation 
that  the  static  limit  Is  a  good  approximation  for  pulse  widths  greater  than 


or  equal  to  one  plasma  period,  l.e. ,  for  M  >_  tp  s  2ir/wp  In  this  region 
the  Josephson  pulse  detector  Is  linear,  by  which  we  mean  that  at  threshold 
a  change  In  ip1  can  be  compensated  for  by  an  equal  but  opposite  change  in 
1b  following  the  relationship  (1fa  +  1pl)Tft*  1-  1°  the  next  section  we 

argue  that  t  ,  and  therefore  the  minimum  pulse  width  for  linear  detection, 
can  be  as  small  as  2.1  ps.  For  narrower  pulses  nonlinear  threshold  conditions 
apply  as  described  by  Eq.  (39)  for  intermediate  pulse  widths  and  Eq.  (34)  for 
the  narrowest  pulses. 


B.  Pulse  Width  Determination 

The  Josephson  pulse  detector  described  in  the  preceding  section  Is 
a  very  simple  device.  Given  such  a  device  one  can  measure  Its  critical 
current  density  and  from  that  infer  Its  R^Cj  and  from  Fig.  7.  However, 
If  this  detector  is  used  to  detect  pulses  of  unknown  width  and  amplitude, 
what  can  be  learned  about  the  pulses?  About  all  one  can  measure  Is  the 
bias  threshold  for  detection  of  the  pulses,  but  that  tells  very  little 
about  them  since  the  switching  may  be  due  to  wide  pulses  of  relatively 
low  amplitude  or  from  narrower  pulses  of  higher  amplitude.  Knowledge  of 
the  detector  response  time  does- not  change  the  situation.  However  If  the 
pulse  width  could  be  measured,  the  pulse  amplitude  could  be  inferred. 

Before  describing  a  technique  for  pulse  width  measurement.  It 
Is  Instructive  to  consider  a  hypothetical  situation  In  which  two  tri¬ 
angular  shaped  current  pulses  are  applied  to  an  Ideal  threshold  detector. 

For  simplicity  the  widths  W  and  the  amplitudes  I  the  two  pulses  are 

pl 

equal.  Figure  14  (a)  Is  a  plot  of  current  In  the  detector  as  a  function 
of  time  for  a  particular  delay  t  between  the  centers  of  the  pulses.  The 
Idealized  detector  Is  a  device  which  Instantly  switches  the  state  of  Its 
output  terminals  (and  remains  switched}  if  an  adjustable  current  level 
I4  ts  exceeded  at  Its  Input  terminals.  Consider  the  Input  to  the  detector 
to  be  repeatable  and  the  detector  to  be  resettable  so  that  one  can 
experimentally  determine,  for  a  given  t,  the  maximum  value  of  Id  for 
which  threshold  Is  crossed  by  the  Input  signal. 


The  result  of  such  an  Idealized  experiment  Is  shown  In  Fig.  14  (b).  If 
the  pulses  are  well  separated  In  time  the  maximum  value  of  1^  for  detection, 
Id(Th),  Is  equal  to  Ip^ .  As  t  Is  decreased  the  pulses  begin  to  overlap  but 
I^CTh)  does  not  change  until  one  pulse  begins  to  overlap  the  center  of  the  other 
pulse.  The  maximum  value  of  Id(Th)  is  2Ip^  and  occurs  for  t  *  0.  The  result 
of  this  Idealized  experiment  is  an  Image  of  one  of  the  pulses  in  the  Id (Th )  vs. 
t  plot.  The  Image  pulse  has  an  offset  baseline  but  the  pulse  shape  Is  unchanged. 

The  scheme  for  using  Josephson  junctions  for  pulse  width  measurements  Is 
similar  to  the  Idealized  case  just  described  but  there  are  significant  dif¬ 
ferences.  The  idealized  device  Is  linear  In  that  a  change  in  the  bias 
Id  produces  an  exactly  equal  change  In  the  current  required  to  reach  thres¬ 
hold.  For  a  Josephson  junction, that  Is  true  only  In  the  static  limit. 

In  both  the  square  pulse  approximation  Eq.  (39)  and  the  Impulse  approximation 
Eq.  (34)  the  relationships  between  Ifa  and  1^  at  threshold  are  nonlinear. 

With  th'e  Idealized  detector  the  fidelity  of  the  image  of  the  driving 
pulse  is  a  result  of  the  linearity  of  the  detector.  With  a  nonlinear 
detector  distortion  In  the  pulse  Image  can  be  expected. 

Our  Investigations  of  this  procedure  are  limited  to  a  circuit  like 
the  pulse  detection  circuit  of  Fig.  13  with  the  addition  of  a  second 
Gaussian  pulse  source.  We  confine  our  attention  to  calculating  the 
switching  threshold  of  the  detector  as  a  function  of  time  delay  between 
two  pulses  for  delays  of  up  to  a  few  times  the  pulse  width.  In  essence 
this  Is  a  study  of  the  Interference  pattern  of  two  pulses  applied  In 
coincidence,  or  near  coincidence,  to  a  nonlinear  detector.  This  Is  a 
special  case  of  the  waveform  measurement  technique  shown  In  Fig.  3  In 
which  the  signal  Is  Identical  to  the  sampling  pulse. 


For  simplicity  the  two  Gaussian  pulses  have  the  same  amplitudes, 

*pl  *  *p2  *  0.33,  and  the  same  widths  W^«W2»W.  Our  first  simulations  of 
this  type  of  experiment.  Fig.  15  (b),  produced  a  surprisingly  broadly- 
structured  Interference  pattern.  In  this  plot  of  Ib(Th)  versus  delay 
time  between  first  and  second  pulses,  the  image  of  the  pulse  Is  a 
negatively  directed  peak  at  zero  delay.  Only  the  pattern  for  positive 
delay  times  is  displayed  since  the  pattern  is  symmetric  about  zero 
delay.  This  case  Is  for  W  *4.72  ps,  rp  *  11.5  ps,  and  R^Cj  »  243  ps. 

Note  that  the  width  of  the  central  minimum  in  Fig.  15  (b)  is  WQ  »  9.6  ps. 

For  the  present  we  Ignore  the  broad  structure  and  ask  how  does  WQ  change 
as  the  width  of  the  applied  current  pulses  is  changed?  In  Fig.  15  (c) 

W  has  been  approximately  doubled  to  9.54  ps  and  It  Is  found  that  the 
central  minimum  did  Increase  In  width  but  only  by  about  382  Instead  of 
the  anticipated  factor  of  2.  Thus  It  Is  clear  that  the  detector  Is  not 
responding  in  a  simple  way,  and  that  the  broad  structure  In  the  interference 
pattern  cannot  be  Ignored. 

Evidently  this  structure  Is  due  to  plasma  oscillations.  These  oscillations 
are  less  damped  here  than  In  Fig.  12  since  P,  Is  now  larger.  In  Fig.  15  (a) 
the  time  dependence  of  $  Is  shown  with  only  the  first  of  the  two  pulses 
applied.  It  Is  Immediately  obvious  that  the  ringing  behavior  Is  present  for 
the  time  delays  of  Fig.  15  (b)»  and  Is  of  the  correct  frequency  to  cause 
the  observed  structure.  Thus  our  conclusion  Is  that  the  plasma  oscillations 
are  so  strong,  particularly  In  Fig.  15  (b],  that  the  system  Is  not  very 
sensitive  to  the  applied  pulse  width. 
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In  Fig.  15  (d)  the  width  of  the  pulse  Is  Increased  to  19.02  ps  and 
the  plasma  oscillations  are  clearly  less  important  than  in  the  preceding 
two  figures.  It  is  revealing  to  calculate  the  bandwidth  of  this  pulse. 

Using  Eq.  (27)  the  result  is  3  23  GHz.  This  Is  low  compared  with 
the  plasma  frequency  of  the  detector  which  Is  about  87  GHz.  Consequently 
there  is  very  little  energy  available  to  excite  the  plasma  oscillations. 

Thus  it  becomes  clear  why  the  extended  structure  in  Fig.  15  (b)  diminishes 
as  wider  pulse  widths  are  used. 

The  oscillatory  behavior  of  the  interference  pattern  for  large  delay 
times  is  much  reduced  in  Fig.  15  (d)  as  compared  with  the  graph  immediately 
above  it.  At  the  same  time  we  note  that  WQ  changes  essentially  by  a  factor 
of  2  between  these  two  figures.  Thus  the  width  of  the  pulse  image  now 
changes  in  proportion  to  changes  in  the  driving  pulse  width,  which  Is  the 
desirable  result.  This  behavior  is  correlated  with  the  pulse  width  W  becoming 
comparable  to  and  greater  than  the  period  of  a  plasma  oscillation  Tp. 

To  further  substantiate  this  point  of  view  additional  simulations 
were  done  for  the  narrow  pulse  width,  4.72  ps,  used  for  Fig.  15  (b)  but 
now  applied  to  a  junction  with  a  shorter  plasma  period,  t  »  4.76  ps. 

The  result  Is  the  considerably  lower  level  of  plasma  oscillations  evident 
in  the  Interference  pattern  Fig.  15  (e)  as  compared  with  Fig.  15  (b). 
Furthermore  It  Is  now  found  that  when  the  applied  pulse  width  is  increased 
6y  a  factor  of  2,  to  9.43  ps,  the  width  of  the  pulse  image  also  Increases 
by  a  factor  of  2,  as  is  seen  by  comparing  Figs.  15  (e)  and  (f). 
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Thus  the  interference  pattern  provides  a  sensitive  means  of  measuring 
pulse  width  if  W  >_  t  .  Me  note  in  passing  that  our  empirically  derived 
criterion  for  sensitive  pulse  width  measurements  can  also  be  expressed  in 
the  frequency  domain.  Using  Eq.  (38)  we  obtain  for  Gaussian  pulses,  fp^ 


1/2 


>.  2.27,  where  f  Is  the  plasma  frequency  in  Hz.  This  result  is  similar  to 
our  earlier  observation  that  if  the  plasma  frequency  is  well  above  the 
spectrum  of  the  applied  pulse,  then  the  plasma  oscillations  are  not  highly 
excited  and  simple  behavior  results.  However  the  width  of  the  pulse  image 
Mq  is  not  necessarily  numerically  equal  to  the  pulse  width  M.  The  widths 


of  the  interference 
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patterns  in  Figs.  15  (d)  and  ( f )  are  about  1.3  times  the  driving  pulse 
widths.  It  is  expected  that  for  larger  pulse  widths  the  pulse  image  will 
more  nearly  approximate  the  driving  pulse.  Without  further  theoretical 
development,  detailed  simulations  are  required  to  determine  the  driving 
pulse  width  from  the  Interference  pattern.  An  experiment  of  this  type  has 
been  done  by  Paris  in  which  a  pulse  width  of  26  ps  was  measured  with  little 
distortion  (28). 

We  are  now  in  a  position  to  consider  the  minimum  measureable  pulse 

width  for  Josephson  detectors  driven  by  current  sources.  The  previously 

2 

discussed  maximum  current  density  of  28,000  A/cm  corresponds,  in  Fig.  7, 
to  a  plasma  period  of  2.1  ps.  Therefore  the  minimum  measureable  pulse 
width  is  also  2.1  ps  if  the  change  in  the  width  of  the  pulse  Image  must  be 
proportional  to  the  change  in  the  width  of  the  pulse.  Narrower  pulses 
can  be  measured  If  a  less  sensitive  relationship  between  image  width  and 
pulse  width  Is  acceptable.  Whether  or  not  that  might  be  the  case  depends  on  the 
experimental  signal -to-nolse  ratio  which  Is  beyond  the  scope  of  the  present 
article. 

VI.  Superconducting  Striplines  at  High  Frequencies 

In  previous  sections  we  have  detailed  the  speed  limitations  of 
the  Josephson  junction  In  various  applications.  A  limitation  which 
remains  to  be  considered  Is  that  imposed  by  the  Interconnections 
between  junctions  In  circuits. 

As  was  realized  early  on  (57),  the  operation  of  Josephson  circuits 
at  high  speeds  requires  that  the  junctions  be  Interconnected  using 
transmission  lines  with  matched  loads.  This  arrangement  prevents  multiple 
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reflections  of  the  signals  and,  assuming  an  ideal  transmission  line, 
allows  full  utilization  of  the  devices'  switching  speed.  Swihart  has 
shown  (58)  that  at  frequencies  well  below  the  superconducting  energy-gap 
frequency  the  striplines  used  to  interconnect  superconducting  microcircuits 
are  nearly  ideal,  having  low  loss  and  low  dispersion.  However,  at 
frequencies  approaching  the  energy  gap,  both  attenuation  and  dispersion 
can  be  significant  for  striplines  having  lengths  typical  of  microcircuit 
dimensions  (59). 

The  speed  limitations  Imposed  by  superconducting  striplines  were 
explored  theoretically  by  Kautz  (59)  using  the  theory  of  Mattis  and 
Bardeen  for  the  conductivity  of  a  superconductor.  Figure  16  shows  the 
attenuation  and  phase  velocity  of  an  example  stripline  chosen  to  be 
typical  of  those  used  in  superconducting  microcircuits.  Results  are 
shown  for  temperatures  of  4.2  K  at  which  the  line  Is  superconducting  and 
9.2  K  at  which  the  line  is  normal.  Above  the  energy  gap  frequencies  of 
the  Nb  and  Pb-alloy  conductors,  720  and  610  SHz  respectively  at  4.2  K,  the 
attenuation  and  phase  velocity  are  nearly  the  same  for  the  superconducting 
and  normal  states.  Well  below  the  gap  frequency,  the  superconducting 
line  shows  Its  superiority,  having  an  attenuation  orders  of  magnitude 
lower  than  the  normal-state  line  and  virtually  no  dispersion.  Below 
about  10  GHz  the  attenuation  and  dispersion  of  the  superconducting  line 
are  low  enough  that  lines  1  cm  or  less  In  length  can  be  considered 
Ideal.  At  frequencies  approaching  the  gap,  however,  both  attenuation 
and  dispersion  become  significant  even  for  a  1  cm  length  of  line. 


45 


/ 


The  effects  of  attenuation  and  dispersion  at  these  high  frequencies 
on  pulses  Is  shown  In  Fig.  17.  Here  we  see  the  evolution  of  two  pulses 
which  are  Initially  Gaussian,  with  widths  W  »  2.4  and  1.2  ps,  as  they 
propagate  along  the  example  stripline.  For  these  pulses,  virtually  all 
of  the  energy  Is  contained  In  frequency  components  below  the  energy  gap, 

99%  of  the  energy  being  below  290  and  580  GHz  for  the  2.4  and  1.2  ps 
pulses,  respectively.  As  Fig.  17  shows,  the  2.4  ps  pulse  can  travel  no 
further  than  about  0.1  cm  without  significant  distortion,  while  the 
limit  for  the  1.2  ps  pulse  is  about  0.01  cm.  Thus,  for  microcircuit 
applications,  the  quality  of  the  example  line  is  excellent  for  pulse 
widths  down  to  a  few  picoseconds  but  for  shorter  pulses  becomes  rapidly 
unsatisfactory. 

Finally,  we  examine  the  extent  to  which  the  high-frequency  limit  of 
the  example  stripline  can  be  Increased  by  changing  Its  materials,  tempera¬ 
ture,  or  dimensions.  Use  of  a  superconductor  with  a  larger  energy  gap, 
for  example,  leads  to  a  direct  Increase  In  the  high  frequency  limit. 
Unfortunately,  such  large-gap  superconductors  are  refractory  and  difficult 
to  work  with.  Lowering  the  operating  temperature  leads  to  a  dramatic 
decrease  In  attenuation  below  the  gap  frequency  but  reduces  the  dispersion 

only  slightly.  If  the  example  strlpllne  Is  operated  at  1.0  K  then  the 
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sub-gap  attenuation  does  not  exceed  10  dB/m  but  dispersion  at  frequencies 
just  below  the  gap  still  prevents  propagation  of  sub-picosecond  pulses 
over  significant  distances  (59).  Finally,  because  both  the  attenuation 
and  the  frequency  dependent  terms  of  the  phase  velocity  are  Inversely 
proportional  to  the  dielectric  thickness,  use  of  a  thicker  dielectric 
yields  an  overall  Improvement  In  pulse  propagation  (60).  Increasing  the 


dielectric  thickness  by  a  factor  of  10,  for  example,  provides  a  10-fold 
Increase  In  the  distance  a  given  pulse  can  propagate  before  It  becomes 
severly  distorted.  However,  for  the  level  of  miniaturization  required 
In  a  complex  microcircuit,  the  dielectric  can  probably  be  no  thicker 
than  1  ym  and  sub-picosecond  pulses  will  propagate  no  more  than  a  cm 
without  distortion.  Thus,  the  energy  gap  frequency  represents  a  firm 
upper  limit  for  pulse  propagation  In  Superconducting  microcircuits. 
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VII.  Discussion 

The  behavior  of  superconducting  circuits  In  the  picosecond  regime 
has  been  studied  In  only  a  few  expelments.  Consequently,  In  assessing 
picosecond  applications  of  these  circuits,  we  must  rely  heavily  on 
circuit  simulations.  There  are  two  models  available  for  simulating 
the  active  part  of  the  circuit:  the  RCSJ  model  and  the  microscopic 
model . 

The  RCSJ  model  Is  a  simplification  of  the  microscopic  model  and 
Is  used  almost  exclusively  In  the  work  reported  here,  partly  because 
It  requires  much  less  computer  time,  but  also  because  It  is  technically 
adequate.  To  establish  Its  adequacy  we  compared  it  with  frequency 
domain  experiments  In  Section  II. A.,  and  we  compared  It  with  the 
microscopic  model  In  two  areas:  the  evaluation  of  the  Josephson 
Inductance  in  Section  III.B  and  the  calculation  of  picosecond  wave¬ 
forms  In  Section  V. A.  We  summarize  these  results  as  follows:  for 
electrical  responses  that  occur  in  time  intervals  of  a  few  picoseconds 
or  longer,  the  RCSJ  model  Is  adequate  for  practical  purposes.  However, 
as  shorter  time  periods  are  considered  It  should  be  recognized  that  the 
RCSJ  model  Is  Increasingly  Inaccurate,  with  waveform  errors  as  large 
as  approximately  35%  at  0.8  ps.  Consequently  we  conclude  that  all  of 
the  time-domain  experiments  that  we  discussed  are  justlfably  analyzed 
with  the  RCSJ  model. 

Within  the  limitations  of  the  RCSJ  model  we  examined  the  step 

function  response  of  a  junction  and  derived  a  simplified  formula  for  the 

turn-on  delay  t  associated  with  switching  from  the  zero  voltage  state 
D  ' 

to  the  energy  gap  voltage.  The  theory  of  turn*>on  delay  was  also 


extended  to  Include  arbitrary  Initial  bias  levels  Instead  of  only  the 
single  case  of  Initial  bias  equal  to  the  critical  current.  It  was  shown 
that  the  total  step  height  of  the  bias  current,  rather  than  simply  the 
amount  of  overdrive.  Is  Important  in  determining  xQ.  A  simple  expression 
for  the  delay  until  the  first  appearance  of  a  voltage  pulse  In  a  single 
junction  SQUID  was  also  developed. 

Detailed  consideration  was  then  given  to  the  theory  for  measurements 
of  the  height  and  width  of  pulses  using  a  latching  Josephson  detector. 
Formulas  were  derived  for  the  switching  threshold  condition  for  three 
ranges  of  pulse  widths.  We  found  the  threshold  condition  to  be  especially 
simple  for  pulse  widths  W  >  x  .  where  t  Is  the  period  of  a  plasma 

—  p  p 

oscillation.  The  simple  threshold  condition  Is  a  linear  relationship, 

(1b  +  fpiJ-rj,  *  1.  With  the  RCSJ  model  any  width  of  pulse  will  switch 
the  detector  If  It  has  sufficient  amplitude  but  the  threshold  condtttons 
are  nonlinear  for  W  <  Xp.  The  transmission  lines  limit  the  minimum  pulse 
width  to  about  1  ps. 

Simulations  show  that  pulse  width  measurements  can  be  made  by 
studying  the  detector  response  for  pulses  arriving  In  near  coincidence 
at  the  detector.  It  was  found  that  this  technique  gives  a  sensitive 
measure  of  pulse  width  as  long  as  K  >  tp,  the  same  criterion  as  for  the 
detector  to  be  linear. 

For  pulse  widths  less  than  Xp  the  technique  no  longer  proportionately 
measures  the  width  of  the  pulse  and  becomes  Increasingly  Insensitive  to 
pulse  width  as  the  width  Is  decreased. 


These  results  suggest  a  model  In  which  the  detector  junction  Is  an  Ideal 
Instantaneously  responding  device  In  series  with  a  low-pass  filter  which  Is 
resonant  at  the  plasma  frequency.  The  Q  of  the  filter  Is  determined  by  the 
loading  across  the  junction.  Thus  most  of  the  ringing  reflected  In  Fig.  15  (b) 
can  be  suppressed  by  adding  a  resistor,  but  frequencies  above  <i»p  cannot  be 
readily  observed. 

Thus  we  conclude  that  the  plasma  period  rp  Is  an  approximate  measure 
of  the  highest  time  resolution  that  can  be  achieved  In  time-domain  measure¬ 
ments.  (Obviously  a  more  complete  understanding  must  Include  consideration  of 
signal -to-nolse  ratio.)  Because  of  self-heating  effects  tp  cannot  be  made 
smaller  than  about  2  ps  with  present  fabrication  procedures. 
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Appendix  A. 


Time  Domain  Formulation  of  the  Microscopic  Theory 

The  microscopic  theory  seems  to  be  not  well  known  In  detail,  probably 
because  the  phenomenological  model  has  consistently  produced  accurate 
descriptions  of  experimental  circuits.  Two  things  may  chnage  this  situation 
tn  the  future.  First  the  theory  has  been  recast  Into  a  time  domain  formulation 
which  makes  numerical  computation  simpler  (61).  Second,  the  Increasing  speed 
of  practical  circuits  may  require  simulations  of  circuits  for  which  the 
phenomenological  model  Is  too  Inaccurate.  For  these  reasons  we  Include  in 
this  Appendix  a  brief  description  of  the  microscopic  theory. 

The  time  domain  formulation  of  the  microscopic  theory  carries  present 
and  past  Information  about  the  junction  voltage  In  a  function  U(t): 


U(t)  »  exp  (Hs(t)) 


exp  ( 


'Hr  J  v(t' 


)dt'), 


CA.l) 


The  function  U(.t)  Is  very  much  like  a  single-particle  wave  function.  It  Is 
a  solution  of  the  time-dependent  Schroedlnger  equation  having  Hamiltonian  -eV 
where  V  is  the  voltage  across  the  junction. 

The  theory  describes  the  junction  current  in  terms  of  U(t): 


I(t) 


Iqp(t-t')  U(t')  dt' 


+ 


Ij(t-t')  U(t')  dt'}. 


(A. 2) 
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Here  the  kernels  Iqp(t)  and  I jCt )  produce  the  details  of  the  time-dependent 

★ 

response  of  the  junction.  The  U  denotes  the  complex  conjugate  of  U. 

We  now  examine  the  form  of  the  two  kernels  which  are  plotted  in  Fig.  18. 

First  Ij  which  describes  the  supercurrent  reveals  a  logarithmic  singularity 
(negative)  just  after  time  zero.  This  Integrable  singularity  represents  a 
very  fast  but  not  Instantaneous  response.  If  the  time  scale  were  compressed, 
by  letting  the  energy  gap  2a  become  large,  the  singularity,  combined  with 
Eq.  CA.2),  transforms  into  the  Instantaneous  current  IQs1n<f>  of  the  phenomenological 
model.  Somewhat  after  time  zero,  Ij  reveals  oscillations  at  the  energy  gap 
frequency  2a/h.  These  correspond  to  the  resonance  known  as  the  Riedel  peak 
shown  In  Fig.  5. 

It  Is  also  Interesting  to  observe  that  the  second  term  of  Eq.  (A. 2)  Involving 
I.  Includes  a  product  of  the  form  U(t)U(t').  This  form  tends  to  double  the  phase 
$s  making  It  Into  the  pair  phase  of  Eq.  (2).  In  .contrast  the  first  term  of 
Eq.  (A. 2)  contains  the  product  U*(t)U(t')  which  does  not  double  the  phase  <j>s . 

Thus  the  second  term  corresponds  to  the  flow  of  paired  electrons  (the  super¬ 
current},  while  the  first  term  describes  the  flow  of  single  (normal  or  quasi¬ 
particle)  electrons. 

The  kernel  I  is  slightly  more  Involved.  This  kernel  contains  all  the 

qp 

Information  found  In  the  current-voltage  characteristic  of  Fig.  5.  First 
we  note  that  the  characteristic  becomes  a  linear  resistance  RN  for  voltages 
above  the  energy  gap.  Since  a  pure  resistance  describes  an  Instantaneous 
response.  Is  appears  In  Iqp  as  a  delta  function.  To  simplify  Iqp  we  there¬ 
fore  separate  out  the  Instantaneous  response. 

Iqp(t)  -  -2  0<eRN)  «'Ct)  +  IqpQ(t)  .  (A. 3) 
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The  portion  of  the  kernel  I  Q  Is  shown  In  Fig.  18.  It  also  contains  a 
peak  just  past  time  0.  This  peak,  however.  Is  not  as  sharp  as  that  In  Ij 
and  thus  corresponds  to  a  somewhat  slower  response.  The  oscillations  at 
the  gap  frequency  In  Ij  also  appear  In  and  correspond  to  the  abrupt 

rise  In  current  of  the  voltage-dependent  resistance  at  the  energy  gap. 

Although  the  time-domain  formulation  of  the  microscopic  theory  Is 
conceptually  equivalent  to  that  In  the  frequency-domain.  It  Is  particularly 
useful  In  numerical  simulations  (61].  As  devices  are  fabricated  which  are 
faster  than  those  commonly  used  today,  simulating  them  may  necessarily 
Involve  this  formulation  of  the  theory. 

The  time-domain  formulation  could  be  significantly  Improved  in  one 
respect.  The  energy  gap  frequency  oscillations  In  both  Ij  and  IqpQ  maintain 
significant  amplitude  to  quite  long  times.  This  Is  a  result  of  the  sharp 
gap  structure  seen  In  the  current-voltage  characteristics.  In  practice,  this 
structure  Is  smeared  over  a  range  of  voltages.  As  a  result  the  kernels  Ij  and 
IqpO  for  Pract-,cal  Junctions  should  correspondingly  exhibit  a  faster  fall  off 
In  the  amplitude  of  the  energy  gap  oscillations.  The  development  of  a  physical 
model  for  this  fall -off  would  make  It  possible  to  truncate  the  Integrals  of 
Eq.  (4)  after  a  suitable  time.  This  would  speed  numerical  computations  and 
reduce  memory  requirements. 
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Appendix  B. 


Turn-On  Delay  in  a  Single- Junction  SQUID 

Although  not  at  the  moment  of  apparent  practical  Importance,  turn-on 
delay  In  a  single- junction  SQUID  Is  interesting  in  its  own  right.  Me  shall 
sketch  briefly  how  one  can  estimate  tg  for  this  case.  Many  mathematical 
details  for  flux  entry  Into  a  SQUID  may  be  found  in  reference  62.  Consider 
the  case  in  which  the  junction  of  Fig.  4  Is  shunted  by  a  superconducting 
loop  of  Inductance  L.  This  circuit  is  described  by  Eq.  (6)  if  we 
add  to  the  right  side.  The  washboard  potential  is  now  given  a 

positive  net  curvature,  becoming  a  parabola  with  ripples.  The  potential 
energy  describing  this  case  Is 

V  *  cost  +  (t-&b)^i»  (B.l) 

where  *  LIk^o*  APart  from  an  unimportant  constant,  Eq.  (B.l)  goes 
to  Eq.  (17)  as  L  -*■  «.  It  Is  clear  that  this  SQUID  will  have  less  of  a 
chance  for  switching  to  Vg  because  of  the  curvature  of  the  parabola. 

The  time  delay  before  one  or  more  fluxoids  enter  the  SQUID  after  ap 
plication  of  a  step  Increase  In  1^  is  of  potential  Interest,  however. 

As  Ifa  Is  raised  slowly  from  zero,  t  will  advance  from  zero,  always 
staying  In  the  first  local  minimum  given  by  dV/dt*  0.  Thus  t  Is  that 
value  which  solves 


♦  +  S0  sint  *  p5. 


(B.2) 
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When  d  1/d  $  *  0  at  the  minimum,  the  particle  will  spill  out,  corresponding 
to  the  entrance  of  one  Cor  more)  fluxolds  Into  the  SQUID.  Thus,  fluxoid 
entry  occurs  when  cos$  ■  -1/SQ  and  Eq.  (B.2)  are  satisfied  simultaneously. 

The  critical  or  threshold  current  for  this  to  occur  is  given  by 

s£  «  d>c  +  SQs1  o«t»c  *  cos"1  (-s’1)  +  (B2  -  l)1^2.  CB.3) 

Suppose  that  when  this  current  is  reached.  It  is  suddenly  stepped  to 
l£  (1+V-  **  we  then  exPan<*  s,,n<f>  about  initial  value  si n«<>c ,  retaining 

only  the  linear  term  in  $  -  $  ,  we  find  from  the  generalization  of  Eq.  (6) 

w 

to  the  SQUID, 

sb  6o  *  LC  *  +  If  *  ‘  (B.*) 

The  term  with  d>  has  dropped  out,  leaving  the  same  form  of  equation  as  In 
Section  IV  for  the  junction  not  shunted  by  Inductance.  If  we  assert,  as 
there,  that  the  turn-on  delay  time  can  be  approximated  by  the  time  required 
for  the  phase  to  advance  by  one-half  radian,  we  get 

l8-5) 

where  f  ■  s£/S0  *  1  +  w/C2bq1  +  (2&2)-1.  The  approximate  equality  here 
results  from  expanding  the  right  side  of  Eq.  CB.3)  in  powers  of  s’*,  and  Is 
an  excellent  approximation  even  at  the  lowest  value  of  SQ  (unity)  for  which 
the  system  particle  can  be  trapped  away  from  the  bottom  of  the  parabola. 

Equation  (B.5)  should  be  compared  with  Eq.  (18),  which  it  generalizes. 


Numerical  simulations  show  that  Eq.  (B.5)  Is  quite  a  good  approximation 

to  the  time  for  the  rising  of  the  first  voltage  pulse.  The  effect  of  £ 

the  shunt  Inductance  Is  to  shorten  the  turn-on  delay  time. 

c 


Appendix  C 


Turn-On  Delay  in  a  Two-Junction  SQUID 

A  two-junction  SQUID,  or  interferometer,  is  depicted  in  Fig.  19  (a). 

It  is  evident  that  such  a  device  is  Inherently  slower  than  a  single 
junction,  principally  because  the  two  junctions  will  not  necessarily  act 
in  concert. 

There  is  a  trivial  exception  to  this  which  we  shall  mention  here, 
because  the  Idea  occurs  again  later.  When  the  Interferometer  is  completely 
symmetric,  and  there  Is  no  control  current  Ic,  the  bias  current  Ifa  divides 
equally  between  the  two  halves  of  the  interferometer.  The  Inductances  L  and 
damping  resistors  draw  no  currents.  When  Ib  is  stepped  to  2Iq (1 +6Q ) ,  one 
half  of  this  passes  through  each  junction.  The  junctions  act  together,  and 
the  turn-on  delay  is  given  by  Eqs.  (18),  (19),  or  (23),  where  wp  still 
refers  to  the  single- junction  plasma  frequency. 

Of  course,  the  reason  for  using  a  multi-junction  Interferometer  is 

f 

because  the  control  current  allows  one  to  set  the  switching  threshold  value 
of  Ib  over  a  range  of  values,  as  shown  in  Fig.  19  (b).  When  Ic  f  0,  the 
analytical  problem  is  much  more  difficult  than  for  a  single  junction.  The 
only  published  treatment  of  turn-on  delay  In  a  two-junction  SQUID  at  the 
moment  is  that  of  E.  P.  Harris  (.52).  We  shall  here  follow  Harris' 
reasoning  and  give  his  results.  . 

If  the  Interferometer  is  heavily  damped  to  suppress  the  LC  resonance 
effects  of  the  circuit,  the  interferometer  Is  described  approximately 
by  the  two  equations 
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+  $_  +  60s1rnj>_cos4>+  *  8  , 


Cc.i) 


2LC  <ji+  +  *+  +  2eosin<j»+cos<(»_ 


(C.2) 


where  $+  *  y  U-j*  $2^’  an<1  the  S's  have  been  defined  following  Eq.  (B.l). 

One  sees  from  Eqs.  (C.l)  and  (C.2)  that  If- 6  is  held  constant  at  or  below 
a  threshold  value,  while  ^  is  stepped  up,  then  $+  will  first  respond;  $ 
will  follow  by  virtue  of  the  coupling  between  the  two  equations.  Similarly, 
if  s  is  fixed  and  sr  is  stepped,  <j>  first  responds.  Harris  considered 
only  the  latter  case. 

The  control  current  is  taken  as  zero  for  t  <  0,  and  1^  is  fixed  (for 
all  t)  at  Ibl,  so  that  $_(t  <  0)  »  0  from  Eq.  (C.l).  At  t  -  0,  I  is 

stepped  to  a  value  *02’  greater  than  the  threshold  value  Xcl*  Equation 

« 

(C.2)  shows  then  that  s1n$+  is  given  approximately  by  Bbl sec<j>_/2eo  until 
the  time  at  which  the  latter  quantity  equals  unity.  Harris  calls  this 
time  Interval  T^nt  and  develops  an  approximate  expression  for  it,  given 
below.  At  that  time,  <j>+  begins  Increasing  rapidly,  that  is,  the  two  junctions 
begin  moving  together  and  the  Interferometer  tends  to  behave  as  a  single 
junction.  The  time  for  $+  to  advance  from  w/2  to  3ir/2  is  designated  t  t, 
and  the  total  turn-on  delay  time  of  the  interferometer  is  taken  to  be 
tD  *  Tint  +  Text*  For  Text’  Harr,s  simP1>  adopts  the  single  junction  re¬ 
sult.  In  which  the  capacitance  Is  twice  the  single  junction  capacitance, 
and  IQ  Is  replaced  by  Im,  the  threshold  value  of  Ib  when  Ic  has  Its  final 
value  IC2  (see  Fig.  19  b).  Harris's  approximate  expression  for  is 
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t  +  t,  where 
o  i 
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Figure  Captions 


Fig.  1 


Fig.  2 


Fig.  3 


Experimental  results  for  the  Induced  dc  In  a  Josephson 
junction  with  applied  radiation  at  891  and  805  GHz. 

Spectral  features  appear  at  1.84  and  1.66  mV,  corresponding 
In  Eq.  (3)  to  Josephson  oscillations  at  the  frequencies  of 
the  applied  radiation.  The  feature  at  0.18  mV  is  not  due  to 
applied  radiation  but  corresponds  to  the  difference  frequency 
of  the  applied  signals. 

(a)  Sampling  circuit  for  measuring  the  risetime  of  a  signal 
I$.  A  sampling  oscilloscope  In  a  servo  loop  Is  used  to 
automatically  measure  the  time  at  which  the  signal  plus  the 
bias  cause  Jl  to  switch.  Is  Is  reproduced  by  plotting  the  , 

bias  versus  switching  time,  (b)  Measured  waveform,  using  the 

I, 

above  technique,  of  the  switching  transition  of  a  Josephson  .* 
junction.  >f 

Sampling  circuit  for  measuring  an  arbitrary  signal  waveform 
I$.  A  sharp  pulse  of  amplitude  Is  obtained  by  dif¬ 
ferentiating,  In  an  RC  network,  the  switching  transition 
waveform  of  junction  Jl .  This  pulse  and  the  unknown  signal 
Is  are  applied  to  a  detector  junction  J2  which,  through  1^,  Is 
then  adjusted  to  the  threshold  for  switching.  The  resulting 
value  of  Ib2  is  the  essential  Information  needed  for  determining 
the  value  of  1$  at  the  time  the  pulse  occurs.  The  entire  signal 
waveform  Is  measured  by  sweeping  the  time  at  which  the  pulse 
occurs. 
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Fig.  4.  Schematic  representation  of  the  RCSJ  model,  (a)  The  "J 


symbol  represents  the  supercurrent  path  which  has  current  IQs1  n<p . 

Cj  is  the  Intrinsic  junction  capacitance  and  R  may  be  either  the 

voltage  dependent  resistance  Rj  or  a  parallel  combination  of  Rj 

with  an  external  shunt  resistance,  (b)  The  voltage  dependent 

resistance  Rj  has  three  segments.  Below  the  energy  gap  voltage 

V  the  resistance  Is  R, .  Above  v  Is  a  low  resistance  region 
g  I  9 

which  is  followed  by  a  higher  resistance  designated  as  RN>  the 
junction  normal  state  resistance.  Rj  Is  an  approximation  to  the 
quasiparticle  I-V  characteristic  of  the  microscopic  model. 

Fig.  5.  Results  of  the  microscopic  theory  of  Josephson  junctions.  The 
illustrated  functions  reveal  the  frequency  dependence  of  the 
supercurrent  and  the  voltage  dependence  of  the  quasi  particle 
current  as  described  In  the  text. 

Fig.  6.  Theoretical  Impulse  response  of  a  tunnel  junction  .without 
capacitance.  The  current  calculated  from  the  microscopic 
model  Is  Illustrated  for  a  voltage  Impulse  applied  at  time 
zero.  The  amplitude  of  the  impulse  was  adjusted  to  produce  a 
phase  shift  of  ir/2. 

Fig.  7.  Macroscopic  properties  of  small  Pb  alloy  Josephson  junctions.  The 
critical  current  density  jc,  the  time  constant  RnCj,  the  plasma 
frequency  «0,  and  the  Josephson  penetration  depth  Xj  are  shown 
as  functions  of  tunneling  barrier  thickness  d. 


Fig.  8.  The  potential  energy  washboard,  (a)  The  potential  energy  of 

2 

Eq.  (17),  normalized  to  aip,  Is  plotted  against  $/2ir  for  Ib  =  IQ. 

(b)  Comparison  of  potential  energies  for  three  values  of  I b/ IQ : 
top  curve,  0.9;  middle  curve,  1.0;  bottom  curve,  1.1.  The  circles 
on  the  upper  and  middle  curves  shows  the  Initial  position  and 
potential  energy  of  the  system  "particle".  The  circle  on  the 
lower  curve  lies  at  1/2  radian  beyond  ir/2. 

Fig.  9.  Switching  curve  (from  simulations)  for  a  current  biased  junction, 

stepped  from  Ib  *  IQ  to  1.01  IQ,  Illustrating  both  the  turn-on  delay 
and  rise  time  portions.  IQ  *  0.0796  mA,  Cj  =  1.0  pF,  Vg  *  2.5  mV, 

R-!  *  200  n,  Rl  ■  «.  The  method  used  for  determining  -tp  from  the 
switching  curve  Is  shown.  The  10-90%  voltage  levels  on  the 
risetime  portion  are  also  indicated.  For  weak  coupling  R^  * 
irV^/4IQ  -  24.7  fl,  so  RnCj  »  24.7  ps.  Compare  with  tr  (10-90%). 

Fig.  10.  Turn-on  delay  as  determined  from  Eq.  (18)  (upper  curve),  Eq.  (19) 
(lower  curve),  and  from  simulations  (circles).  The  junction  para¬ 
meters  are  those  of  Fig.  9.  1^  Is  stepped  from  IQ. 

Fig.  11.  One  period  of  the  voltage  across  the  load  resistance  of  the  microwave 
driven  pulser  with  an  applied  frequency  of  10  GHz.  The  microscopic 
theory  of  junctions  was  used  to  calculate  this  waveform.  The 
resulting  pulse  width  Is  2.2  ps  with  the  following  circuit  parameters, 
RnCj  »  3.24  ps, 


( 
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Fig.  12.  Calculated  detector  responses  for  four  different  applied  pulse 


widths.  The  applied  current  pulses  are  shown  at  the  top  of  each 
column  of  Illustration  with  the  corresponding  detector  responses 
arranged  below.  The  detector  properties  and  the  applied  pulse 
heights  (1p^  ■  0.5)  are  the  same  In  all  cases.  The  response 
of  the  junction  voltage  V  and  the  phase  $  (radians)  are  illustrated 
for  bias  levels  just  above  threshold  (labeled  A)  and  for  just 
below  threshold  (labeled  B).  For  the  four  columns,  reading 
from  left  to  right,  the  pulse  widths  are:  324,  32.4,  3.24,  and 
0.324  ps;  and  the  bias  levels  above  and  below  threshold  are: 
0.500215,  0.500214;  0.502139,  0.502138;  0.53905,  0.53904; 

0.95528,  0.95527  mA.  The  detector  properties  are  I  *  1  mA, 

Rn  *  1  .578  Q,  R1  /Rjy  -  5,  Cj  «  2.05  pF. 

Fig.  13.  (a)  Circuit  diagram  for  simulations  of  detector  response.  I 

Is  a  current  source  for  a  Gaussian  pulse  and  Ip  Is  the  constant 
bias  supply.  The  X  symbol  represents  the  Josephson  junction, 
which  Is  described  by  the  RCSJ  model,  (b)  Calculations  of 
detector  threshold  bias  1p(Th)  as  a  function  of  pulse  width  VI. 

For  all  points  in  this  graph  the  pulse  height  is  1pl  •  0.5. 

The  circles,  connected  by  a  solid  line,  represent  data  from 
the  RCSJ  model,  the  squares  are  calculations  from  the  impulse 
response  Eq.  (34),  and  the  triangles  are  calculations  using 
the  square  pulse  approximations  Eq.  (39).  The  maximum  threshold 
current  for  ip^  •  0.5  In  1.5,  the  value  closely  approached  by 
the  Impulse  approximation  at  the  narrowest  pulse  widths.  The 
detector  parameters  are  the  same  as  for  Fig.  12. 
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Fig.  14.  Response  of  an  idealized  linear  detector,  (a)  I  (t)  is  the 

current  to  the  detector.  This  current  has  the  form  of  two  Identical 
triangular  pulses  with  delay  time  r  between  them,  (b)  Threshold 
bias  of  the  detector  Id(Th)  plotted  versus  x,  showing  the  "Image" 
of  the  pulses  of  (a). 

Fig.  15.  Pulse  width  measurement  simulations,  (a)  Time  dependence  of  a 
lightly  damped  plasma  oscillation  excited  by  a  single  pulse  of 
width  W  a  4.72  ps,  with  Its  center  arriving  at  the  detector  at  time 
zero,  (b)  -  (f)  Interference  patterns.  These  patterns  are 
graphs  of  the  threshold  bias  as  a  function  of  delay  time  between 
the  arrival  of  Identical  pulses  at  the  detector.  In  all  cases 
the  pulse  amplitude  1  ^  *  0.33.  The  broadly  structured  Interference 
pattern  of  (b)  is  a  reflection  of  the  plasma  oscillations  In  (a). 

For  illustrations  (b)  -  (d)  the  detector  characteristics  are 
unchanged;  the  applied  pulse  widths  have  the  values  4.72,  9.54,  and 
19.0  ps,  respectively. 

The  Illustrated  values  of  1/2WQ  are  the  half  widths  of  the  "Image" 
of  the  applied  pulses.  The  junction  parameters  are  IQ  3  1  mA, 

Rn  -  1.412  n,  R1/Rn  «  17,  Cj  -  10.1  pF,  and  «  -  0.548  x  1012s_1. 

For  the  Illustrations  (e)  and  (f)  the  plasma  frequency  of  the 
detector  has  been  Increased.  The  applied  pulse  widths  are 
4.72  and  9.43  ps,  respectively.  The  junction  parameters  are 
IQ  ■  1  mA,  R|,j  ■  1.4.2  n,  Rj/Rn  ■  17,  Cj  ■  1.63  pF,  and  Mp  « 

1.37  x  101  V1. 
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Fig.  16.  Theoretical  attenuation  (solid  line)  and  phase  velocity  (dashed 
line)  for  an  example  superconducting  stripline  at  4.2  K  and  for 
the  same  stripline  In  the  normal  state  at  9.2  K.  The  example 
stripline  consists  of  a  0.4  ym  thick  Nb  ground  plane,  a  0.1  ym 
layer  of  NbjOg  dielectric  (assumed  lossless),  and  a  0.4  ym  thick 
Pb  alloy  strip  (assumed  wide  compared  to  the  dielectric  thickness). 

Fig.  17.  Simulated  pulse  propagation  on  the  example  stripline  described  in 
Fig.  16.  The  histories  of  two  pulses  are  shown  In  (a)  and  (b), 
respectively,  by  plotting  the  pulse  shape  at  various  distances  L 
along  the  line.  The  pulses  are  originally  Gaussian  shaped,  at 
L  *  0,  with  W  *  2.4  ps  In  (a)  and  1.2  ps  In  (b). 

Fig.  18.  The  two  kernels  Ij  and  I^q  In  the  time  domain  formulation  of 
the  microscopic  theory. 

Fig.  19.  (a)  Two-Junction  Interferometer,  showing  damping  resistors  Rd 

and  a  load  resistor  R^.  Each  X  represents  the  junction  depicted 
In  Fig.  4  (a),  (b)  Sketch  of  a  typical  threshold  curve  for  a 
two- junction  Interferometer,  showing  control  stepped  from  0  to 
Ic2»  at  bias  current  1^. 
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We  report  the  observation  of  heterodyne  mixing  in  superconductor-insulator-superconductor 
tunnel  junctions  with  significant  conversion  gain  and  with  a  noise  temperature  comparable  to  the 
photon  noise  limit.  (Double  sideband  L  ~ 1  =  1 .40  ±  0. 14,  TM  £  1 .5  K.  at  36  GHz.)  The  mixing 
arises  from  the  strong  nonlinearity  in  the  quasiparticle  tunneling  currents  at  voltages  comparable 
to  the  full  energy  gap.  Gain  is  observed  when  the  onset  of  quasiparticle  current  is  sufficiently 
sharp  that  quantum  effects  are  important.  The  observed  gain  is  in  good  quantitative  agreement 
with  calculations  based  on  the  work  of  Tucker.  It  should  make  possible  the  construction  of 
photon-noise-limited  microwave  heterodyne  receivers. 

PACS  numbers:  07.62.  +  s,  74.50.  +  r 


Applications  such  as  radio  astronomy  and  space  com¬ 
munications  require  coherent  millimeter-wave  receivers 
with  noise  temperatures  that  approach  the  photon  noise  lim¬ 
it  T„  =  hv/k.  The  largest  double  sideband  (DSB)  conver¬ 
sion  efficiency  possible1  from  a  classical  resistive  mixer  is 
L~'  —  1.0.  Therefore  both  a  very  quiet  mixer  and  a  large 
improvement  in  intermediate  frequency  (IF)  amplifier  noise 
would  be  required  to  achieve  this  goal  with  a  classical  mixer. 
In  this  letter  we  report  the  first  observation  of  quasiparticle 
mixing  with  significant  conversion  gain  and  the  first  micro- 
wave  mixing  experiment  to  reach  the  photon  noise  limit. 

Mixing  due  to  the  quasiparticle  currents  in  a  supercon¬ 
ductor-insulator-superconductor  (SIS)  tunnel  junction2  3  is 
not  correctly  described  by  classical  theory  when  the  corner 
on  the  quasiparticle  /-  V  curve  occurs  within  the  voltage 
range  hv/e.  Evidence  for  photon-assisted  tunneling,  ob¬ 
served  in  the  form  of  bumps  separated  by  hv/e  on  the  bias 
voltage  dependence  of  the  IF  signal,2  shows  the  necessity  of  a 
quantum-mechanical  treatment.  Tucker  has  recently  used 
the  full  quantum  tunneling  theory4  to  calculate  the  perfor¬ 
mance  of  an  ideal  SIS  junction  operated  as  a  three-port  mix¬ 
er.3  He  predicts  a  new  physical  phenomenon  of  great  techno¬ 
logical  importance:  quiet  mixing  (TM  ~hv/k )  with  large 
conversion  gain.  This  prediction  suggests  that  SIS  quasipar¬ 
ticle  heterodyne  receivers  can  reach  the  photon  noise  limit. 

Lead  alloy  SIS  junctions  are  convenient  for  mixing  ex¬ 
periments  because  of  their  immunity  to  thermal  cycling  and 
the  availability  of  lithographic  techniques  for  producing 
junctions  with  small  areas.  The  onset  of  quasiparticle  tunnel¬ 
ing  in  these  junctions,  however,  is  usually  smeared  signifi¬ 
cantly  over  the  voltage  range  of  0. 1 5  mV  which  is  equal  to 
hv/e  at  our  operating  frequency  of  36  GHz.  For  our  experi¬ 
ments  we  have  used  Pb-Bi  (29  wt.  %  Bi)  junctions  processed 
to  minimize  this  effect.6  The  resulting  /•  V curve  for  a  typical 
3 X  3/um  junction  with  a  normal  resistance  of  21  fl  is  shown 
in  Ftg.  1(a). 

Mixer  experiments  were  carried  out  using  a  measure¬ 
ment  system  which  was  designed  for  the  quantitative  evalua¬ 
tion  of  mixer  parameters,2  rather  than  for  use  as  a  low-noise 
receiver.  The  junction  was  deposited  on  a  0.28  X  l.Ox  19 


mm  Si  strip  which  was  placed  in  the  E-field  direction  across 
a  full-height  XX -  band  microwave  waveguide.  An  adjustable 
stub  J/i,  in  front  of  the  junction  and  a  plunger  behind  it 
formed  the  low-temperature  if  impedance  matching  struc- 


FIG.  1.  Static /-^curves  are  shown  in  the  upper  plot  for  a  junction  at  l.SK 
(a)  without  PLO  and  (b)  with  PLO .  Plots  of  IF  output  power  in  the  frequency 
ran|e  from  30  to  80  MHz  are  shown  in  the  lower  portion  of  the  flgurt  as  a 
function  of  junction  bias  voltage.  Point  (c)  was  obtained  with  a  50-/7  1 ,5-K 
load  in  place  of  the  mixer,  curve  (d)  with  a  matched  t.S-K  load  in  front  of 
the  mixer,  and  curve  (e)  with  a  calibrated  36-GHz  signal  applied  to  the 
mixer  from  an  oscillator.  Values  of  mixer  noise  temperature  were  deduced 
from  (c)  and  (d),  and  conversion  efficiency  from  (e).  Curve  (f)  shows  the 
oscillations  in  the  mixer  output  calculated  from  quantum  three-port  Y- 
mixcr  theory  with  no  free  parameter.  Below  —2.1  mV  Joeephson  (pair 
tunneling)  currents  degrade  mixer  performance.  In  this  region  the  mixer 
has  high  conversion  efficiency,  but  also  large  none. 
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ture.  Radio  frequency  leakage  was  minimized  by  the  use  of 
choke  structures  on  the  plunger,  stub  adjustment,  and  IF 
ports.  Signal  and  LO  power  from  carefully  calibrated  36 
GHz  oscillators  were  combined  in  a  10-dB  directional  cou¬ 
pler  and  introduced  into  the  cryostat  through  a  section  of 
stainless-steel  waveguide.  Noise  from  the  oscillators  and 
from  room-temperature  objects  was  reduced  by  the  use  of  a 
19.4-dB  cold  attenuator  at  the  entrance  to  the  mixer.  The 
combined  loss  of  the  waveguide  and  the  attenuator  was  de¬ 
duced  from  measurements  of  their  transmittance  as  a  func¬ 
tion  of  the  height  of  the  helium  in  the  cryostat.  Three  identi¬ 
cal  50-/7  stainless-steel  IF  cables  were  installed  in  the 
cryostat.  One  was  terminated  with  a  short  circuit  to  measure 
IF  cable  loss,  one  was  terminated  in  a  50-/7  cold  load  to 
provide  a  noise  source  for  calibrating  the  IF  amplifier  train, 
and  one  was  connected  to  the  mixer.  A  directional  coupler 
was  used  to  inject  a  signal  into  the  mixer  output  to  evaluate 
its  coupling  to  the  IF  system.  Measurements  were  made  in 
the  IF  frequency  band  of  30-80  MHz  using  a  string  of  tran¬ 
sistor  amplifiers  with  Ttf  —  50  K,  followed  by  a  spectrum 
analyzer. 

The  static  /-  V  curve  of  a  typical  Pb-Bi  junction  is 
shown  without  LO  drive  in  Fig.  1(a)  and  with  the  LO  power, 
stub,  and  plunger  settings  which  maximize  the  conversion 
efficiency  in  1(b).  Photon-assisted  tunneling  steps  are  weak¬ 
ly  seen  on  the  I-  V curve  measured  with  LO  drive.  These  steps 
appear  as  very  pronounced  bumps  on  the  plot  of  mixer  signal 
versus  bias  voltage  in  Fig.  1(e).  Measured  parameters  for 
mixers  nade  from  very  similar  junctions  are  given  in  Table  I. 
The  stub  and  plunger  settings  and  the  bias  levels  were  chosen 
to  maximize  the  mixer  signal  (and  therefore  the  conversion 
efficiency).  When  the  stub  and  plunger  were  adjusted  in  such 
a  way  that  the  conversion  efficiency  was  the  same  in  both 
sidebands,  the  overall  single-sideband  (SSB)  conversion  effi¬ 
ciency  referred  to  the  entrance  of  the  mixer  block  was  ob¬ 
served  to  be  as  large  as  Z.  “'  =  0.55  ±  0.05.  When  corrected 
for  the  IF  cable  transmittance  (0.89),  and  the  IF  mismatch 
(0.90  ±  0.02),  the  best  mixer  conversion  efficiency  was  esti¬ 
mated  to  be  L  =  0.70  ±  0.07  SSB  or  1.40  ±  0. 14  DSB. 
The  observed  values  of  L  ~x  are  a  factor  —4  larger  than 
those  previously  measured  in  the  same  mixer  mount  with  a 
junction  that  was  operating  close  to  the  classical  limit.2  The 
obsevation  of  a  DSB  conversion  efficiency  greater  than  1 .0  is 


TABLE  I.  Properties  of  two  SIS  mixers  with  normal  state  resistances  of  22 
ti.  operated  at  1.3  K.  The  noise  temperatures  and  conversion  efficiencies  are 
referred  to  a  single  sideband.  Values  of  conversion  efficiency  are  given  sepa¬ 
rately  for  upper  (above)  and  lower  (below)  sidebands.  The  second  listed 
mixer  has  been  optimized  for  SSB  performance.  The  experimental  errors 
are  given  as  30%  confidence  limits. 
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unambiguous  evidence  for  conversion  gain.  Despite  the  low 
IF  frequency,  some  image  rejection  could  be  obtained  by 
operating  with  the  stub  close  to  the  position  at  which  its 
impedance  is  singular.  Mixer  conversion  efficiencies  as  large 
as  0.91  (SSB)  were  observed  in  this  mode  of  operation,  as 
shown  in  Table  I. 

Because  the  nonlinearity  of  the  /-  V  curve  extends  over  a 
narrow  range,  small  values  of  LO  power  are  used,  and  the 
mixer  can  be  saturated  by  relatively  small  amounts  of  signal 
power.  For  a  signal  of  3x  10 W  (~0.3%  of  PLO )  the  IF 
power  was  3%  less  than  the  value  expected  if  there  were  no 
saturation.  Measurements  of  mixer  parameters  were  made 
with  low  enough  power  that  saturation  was  not  important. 
In  a  typical  heterodyne  application  with  a  1-GHz  signal 
bandwidth  B,  the  room-temperature  noise  P  —  kTB 
—  4x  10~12  W  would  not  be  sufficient  to  saturate  the  mix¬ 
er.  Rectification  (direct  detection)  of  this  noise  would  cause 
a  very  small  ( 1  part  in  104)  change  in  the  dc  bias  current.7 

Because  the  onset  of  quasiparticle  tunneling  is  still  sig¬ 
nificantly  rounded  in  the  junctions  used  for  our  experiments, 
Tucker’s  idealized  calculation  of  mixer  properties  is  not  di¬ 
rectly  applicable.  We  have  extended  the  calculations  so  that 
the  imaginary  and  real  parts  of  the  junction  response  are 
obtained  from  the  measured  static  I-  V  curve  and  its 
Kramers-Kronig  transform.1  Predictions  of  conversion  effi¬ 
ciency  and  noise  as  a  function  of  bias  voltage  were  then  com¬ 
puted  for  a  three-port  T-mixer  model  which  assumes  that  all 
harmonic  frequencies  are  shorted  by  the  junction  capaci¬ 
tance.  One  useful  feature  of  our  experiments  is  that  the  shape 
of  the  I-  V  curve  measured  with  various  values  of  LO  power  is 
quite  sensitive  to  the  setting  of  the  stub  and  plunger.  This 
effect  was  utilized  to  deduce  the  rf  source  impedance  Z, .  The 
junction  capacitance  was  regared  as  part  of  this  source  im¬ 
pedance.  The  value  of  Z,  =  14.1  +  jlA  il,  determined  by 
fitting  the  model  to  I-  V curves  measured  with  constant  avail¬ 
able  LO  power,  was  in  good  agreement  with  standard  cavity, 
theory  and  was  used  in  the  calculations  of  mixer  perfor¬ 
mance.  Some  evidence  for  inadequacy  of  the  three-port 
model  was  provided  by  this  procedure.  The  quality  of  the  fit 
was  less  good  for  the  junctions  used  here  than  for  other  junc¬ 
tions  with  larger  capacitance.  The  IF  load  impedance  for  our 
matched  FET  amplifier  was  50 12  (real). 

The  predictions  of  quantum  three-port  T-mixer  theory 
with  no  free  parameters  are  shown  in  Fig.  1.  The  calculated 
mixer  signal  as  a  function  of  bias  voltage  with  constant  rf 
source  impedance  shows  the  same  oscillations  that  appear  in 
the  experimental  data.  To  avoid  confusion,  only  parts  of 
each  oscillation  are  plotted  as  Fig.  1(0-  The  agreement  with 
the  experimental  mixer  output  shown  in  Fig.  1(e)  is  very 
good  where  the  experimental  output  is  high,  but  diverges 
rapidly  as  the  bias  voltage  approaches  3  mV.  We  can  under¬ 
stand  this  striking  result  in  terms  of  the  rapidly  increasing 
first-harmonic  content  of  the  LO  waveform  as  the  bias  is 
increased  toward  3  mV.  In  three-port  T-mixer  theory  these 
harmonics  are  rigorously  shorted,  and  so  do  not  degrade  the 
conversion  efficiency.  In  the  experiment,  the  junction  ca¬ 
pacitance  is  not  completely  effective  in  shorting  the  harmon¬ 
ics,  so  a  significant  fraction  of  the  signal  power  is  converted 
to  harmonic  frequencies. 
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This  interpretation  is  confirmed  by  preliminary  five- 
port  Y-mixer  model  calculations  in  which  the  junction  is 
assumed  to  be  terminated  by  its  own  capacitance  at  the  first- 
harmonic  frequency.  The  value  of  the  reactive  element, 

=  —  j5  Cl,  was  determined  by  a  second-order  fit¬ 
ting  to  the  experimental  /-  V curve  with  LO  power.  It  is  with¬ 
in  a  factor  ~2  of  the  value  estimated  from  the  Josephson 
critical  current  density.  Preliminary  calculations  using  the 
five-port  K-mixer  model  indicate  that  the  conversion  effi¬ 
ciency  for  fundamental  mixing  is  reduced  from  the  predic¬ 
tion  of  the  three-port  model.  The  divergence  in  the  conver¬ 
sion  efficiency  shown  at  ~2.8  mV  in  Fig.  1(f)  is  no  longer 
predicted.  The  predictions  of  the  five-port  model  are  within 
a  factor  —  2  of  the  experimental  coversion  efficiency  at  all 
voltages.  It  appears  likely  from  these  results  that  the  SIS 
junction  will  be  a  very  efficient  first-harmonic  mixer  if  it  is 
biased  near  the  point  of  approximate  symmetry  of  the  I-  V 
curve  near  3  mV. 

We  have  every  reason  to  believe  that  the  predicted  gain 
mechansim  has  been  observed  and  that  larger  gain  can  be 
obtained  from  junctions  with  sharper  corners  on  their  /-  V 
curves,  and/or  from  junctions  operated  at  higher  frequen¬ 
cies.  Some  insight  into  the  conversion  gain  can  be  obtained  if 
we  observe  that  the  conversion  efficiency  is  proportional  to 
the  dynamic  resistance  of  the  junction  with  LO  applied.  Di¬ 
rect  calculations  using  our  three-port  mixer  model  show  that 
the  large  oscillations  in  L  arise  primarily  from  the  effects 
of  the  photon  assisted  tunneling  steps  on  the  dynamic  resis¬ 
tance.  In  particular,  gain  occurs  where  the  /-  V  curve  is  near¬ 
ly  horizontal  so  that  this  resistance  is  large.  In  the  extreme 
quantum  limit,  regions  of  negative  resistance  are  predicted 
which  give  infinite  gain  or  oscillations.  This  association  of 
conversion  gain  with  rf-induccd  steps  on  the  I-V curve  is 


very  similar  to  the  situation  previously  observed  in  point- 
contact  Josephson  effect  mixers.9 

Because  of  the  increase  in  conversion  efficiency,  it  has 
been  possible  to  improve  on  our  previous  measurements  of 
SIS  mixer  noise  temperature.  The  values  given  in  Table  I  are 
much  lower  than  the  upper  limits  obtained  previously  using 
the  same  experimental  techniques.2  The  lowest  observed 
mixer  noise  T„  <  1 .5  K  (DSB)  is  equal  to  the  photon  noise 
limit  hv/k  =  1.7  K  within  the  experimental  uncertainty  of 
approximately  ±  2  K.  The  contribution  to  this  noise  from 
the  LO  has  been  computed  from  the  quantum  theory.4  5  We 
obtain  the  value  1 .4  K  from  the  three-port  Y-mixer  model  for 
the  mixer  shown  in  Fig.  1. 

We  are  grateful  to  J.R.  Tucker  for  showing  us  his  calcu¬ 
lations  prior  to  publication.  A  portion  of  this  research  was 
supported  by  the  U.S.  Office  of  Naval  Research. 
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Large  growth-induced  anisotropy  to  preferential  occupation  of  the  iron  sites 
in  garnets 

L.  C.  Luther,  R.  C.  LeCraw,  E.  M.  Gyorgy,  and  S.  L.  Blank 
Bell  Laboratories,  Murray  Hill,  New  Jersey  07974 

(Received  26  December  1979;  accepted  for  publication  26  February  1980) 

Large  growth-induced  anisotropies  {K, )  arising  from  the  preferential  occupation  of  Ir4*  on 
octahedral  sites  have  been  measured.  For  YIG  with  approximately  0.2  Ir  atoms  per  garnet 
formula  unit  and  Mg2,  as  a  charge  compensator,  values  of  /T„  =  +4.5  X 105  ergs/cm3  have  been 
obtained.  The  induced  anisotropy  is  proportional  to  the  Ir4*  present  as  deduced  from  the  optical 
absorption. 

PACS  numbers:  61.50.Cj 


Preferential  occupation  of  the  rare  earths  on  dodecahe¬ 
dral  sites  is  well  known  to  produce  a  uniaxial  (actually  orth¬ 
orhombic)  anisotropy  in  iron  garnets.  In  fact,  most  of  the 
bubble  technology  depends  on  this  mechanism.  In  this  letter 
we  report  on  the  preferential  occupation  of  4 4  and  5 4  ions  on 
the  iron  (octahedral  and  tetrahedral)  sites.  With  Ir4*  and 


Mg2*  substituted  for  Fe3*  in  YIG  we  have  obtained  a  growth- 
induced  uniaxial  anisotropy  as  large  as  4.5  X 103  ergs/cm3 
for  (1 1 1)  films.  This  may  be  compared  with  3  X 105 
ergs/cm3,  the  largest  growth-induced  Ku  reported  thus  far 
for  the  garnets. 1  The  composition  of  this  garnet  was 
Sm,.1Lu,|Fe,Ol2. 
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A  superconducting  sampler  for  Josephson  logic  circuits 

C.  A.  Hamilton,  F.  L  Uoyd,  R.  L  Peterson,  and  J.  R.  Andrews 

Electromagnetic  Technology  Division.  National  Bureau  of  Standards.  Boulder,  Colorado  80303 

(Received  9  July  1979;  accepted  for  publication  20  August  t979) 

A  method  is  described  for  automating  a  technique  which  is  used  to  sample  transition  duration 
(rise  time)  in  superconducting  logic  circuits.  The  method  is  based  on  measuring  the  time  at  which 
a  biased  Josephson  junction  switches  under  the  influence  of  an  applied  signal.  The  system 
transition  duration  is  limited  primarily  by  time  jitter  which  is  estimated  to  be  7  ps.  Transition 
durations  of  as  little  as  9  ps  have  been  observed. 

PACS  numbers:  74.50.  +  r 


The  characterization  of  Josephson  logic  circuits  often 
requires  the  measurement  of  very  fast  low-level  signals.  Such 
measurements  are  difficult  to  perform  because  of  the  loss  of 
bandwidth  for  signal  transmission  to  room  temperature. 
Also,  standard  sampling  oscilloscopes  have  transition  dura¬ 
tions  (rise  time)  of  about  20  ps  whereas  high  current  density 
Josephson  devices  switch  in  times  less  than  10  ps.  It  is  there¬ 
fore  desirable  to  provide  for  high-speed  measurements  di¬ 
rectly  on  the  superconducting  logic  chip.  This  may  be  ac¬ 
complished  by  using  a  biased  Josephson  detector  as  an 
amplitude  discriminator  for  the  signal  to  be  measured.' 1 
This  detector  will  switch  to  the  voltage  state  when  the  signal 
plus  the  bias  exceed  the  detector  critical  current.  Since  the 
bias  and  the  critical  current  are  known,  the  signal  level  at  the 
time  the  detector  switches  can  be  determined.  Thus,  by  mea¬ 
suring  the  delay  in  the  detector  output  as  a  function  of  bias, 
rising  edges  on  the  signal  may  be  mapped  out.  This  has  com¬ 
monly  been  done*  •*  by  recording  traces  of  the  detector  output 
for  a  number  of  values  of  bias.  The  signal  may  then  be  com¬ 
puted  point  by  point.  This  method  is  somewhat  tedious  and 
may  yield  erroneous  results  if  the  time  base  used  drifts  dur¬ 
ing  the  measurements.  Since  time-base  drift  may  be  as  much 
as  several  ps/min,  it  is  a  serious  problem  for  ultrafast 
measurements. 

Typical  sampling  oscilloscopes  used  to  record  the  de¬ 
tector  waveform  have  an  external  time-base  input  which 
may  be  used  to  specify  the  sampling  time.  By  driving  this 
input  from  the  sampled  voltage,  the  sampling  time  can  be 
locked  to  the  steepest  portion  of  the  detector  output  wave¬ 
form.  The  time-base  input,  therefore,  is  a  measure  of  the 
detector  switching  time.  The  signal  waveform  is  thus  repro¬ 
duced  by  sweeping  the  bias  and  recording  the  bias  versus  the 
time-base  input.  This  method  yields  the  signal  waveform  di¬ 
rectly  and  in  a  time  short  enough  to  minimize  drifl  problems. 

Figure  1  shows  in  more  detail  how  this  is  accomplished. 
The  superconducting  integrated  circuit  is  shown  in  the 
dashed  box.  It  consists  of  a  first  junction  J 1  which  serves  to 
produce  a  very  fast  signal  transition  which  is  then  sampled 
by  the  second  detector  junction  /  2.  The  sequence  is  as  fol¬ 
lows:  a  tunnel  diode  step  generator  operating  at  60  kHz  ap¬ 
plies  a  repetitive  200-ps  transition  duration  step  to  junction 
J 1.  Whenitscriticalcurrent/C,  is  exceeded,/ 1  switches  toits 
energy  gap  (2.5  mV)  in  a  time  determined  primarily  by  its 
critical  current  and  capacitance.  This  applies  a  signal  cur¬ 
rent  /,  through  the  8-/2  resistor  to  junction  /  2.  /  2  will  switch 


when  the  current  /,  plus  the  bias  current  Ib  exceeds  Icl .  The 
output  of  J  2  is  transmitted  to  room  temperature  and  sam¬ 
pled  at  time  t.  The  sampled  voltage  is  compared  with  a  refer¬ 
ence  level,  typically  1  mV,  by  applying  the  reference  as  an 
offset  voltage  to  the  sampling  head  amplifier.  The  difference 
signal  from  the  sampling  head  amplifier  is  integrated  by  in¬ 
verting  amplifier  A1  and  used  to  drive  the  sampling  time 
base  input.  Thus,  if  the  sampled  voltage  is  below  the  refer¬ 
ence  level,  the  integrator  output  increases,  advancing  the 
sampling  time  to  the  point  where  the  signal  equals  the  refer¬ 
ence.  If  the  sampled  voltage  is  greater  than  the  reference,  the 
sampling  time  is  decreased  in  the  same  way.  This  servo  loop 
thus  causes  the  sampling  time  to  lock  to  the  point  at  which 
the  rising  edge  of  the  /  2  transition  crosses  the  reference  level. 
The  time  base  input  is  therefore  a  measure  of  the  time  at 
which  J  2  switches,  i.e.,  the  time  at  which  /,  equals  Icl  -  /„. 
As  Ib  is  slowly  swept  by  the  ramp  generator  (  -  30-s  sweep 
time),  the  X-  Y  recorder,  which  plots  -  Ib  versus  sampling 
time,  will  trace  out  the  rising  portion  of  the  transition  in  /, 
Since  the  sampling  rate  (60  kHz)  is  much  greater  than  the 
servo  loop  frequency  response,  noise  and  time  jitter  in  the 
sampling  system  are  averaged.  The  600-ps  delay  line  on  the 
IC  chip  provides  isolation  to  separate  the  /  2  transition  from 
cross  talk  between  input  and  output  lines  at  the  chip  mount. 
The  timing  network  and  analog  switch  pull  the  bias  to  zero 
between  samples  to  reset  the  junctions  to  the  zero  voltage 
state. 

The  advantage  of  this  method  lies  in  the  fact  that  it  can 
measure  transition  durations  which  are  substantially  faster 
than  the  transition  duration  of  either  the  chip  mount  (~  50 
ps)  or  the  sampling  scope  ( ~  20  ps).  The  measurable  transi- 


FIG.  1.  Block  diagram  of  the  transition  duration  measurement  system  uti¬ 
lizing  a  Josephson  junction  as  an  amplitude  discriminator. 


718 


Appi.  Ptiya.  Lett.  98(9),  1  November  1979 


0003-895 1 779/21 071 8-02800.50 


718 


FIG.  2.  Simulation  of  the  signal  voltage  waveform  (solid  curve)  and  its 
reconstruction  (dashed  curve)  based  on  the  use  of  a  second  junction  as  a 
sampling  device.  The  time  scale  of  the  reconstructed  curve  has  been  arbi¬ 
trarily  shifted  to  bring  the  two  curves  together.  The  simulation  parameters 
are  / ,  —  2.9  mA,  C,  =  4.8  pF.  ,  =  1.45  mA,  C,  =  2.4  pF.  The  input 
is  a  1 -mV/ps  ramp. 

tion  duration  is  limited  primarily  by  the  jitter  in  the  sampling 
system.  The  accuracy  of  the  recorded  waveform  is  also  sub¬ 
ject  to  the  assumption  that  the  detector  junction  waveform  is 
constant  and  that  its  delay  depends  only  on  the  instant  at 
which  its  critical  current  is  exceeded.  These  assumptions  will 
break  down  at  sufficiently  high  speeds.  The  most  significant 
error  results  from  the  fact  that  the  turn-on  delay1  of/ 2  is  a 
function  of  overdrive.  Thus  for  low  bias  values  (signal  near 
maximum),  the/  2  transition  has  some  additional  delay.  This 
has  the  effect  of  making  the  recorded  waveform  slightly 
slower  than  the  actual  signal. 

Figure  2  is  the  result  of  a  simulation  of  our  circuit.  The 
solid  curve  is  the  computed  signal  voltage  produced  by  / 1. 
The  dashed  curve  is  the  computed  result  which  the  above 
described  method  would  produce  in  the  absence  of  jitter.  The 
error  caused  by  variable  turn-on  delay  is  evident  in  the  in¬ 
creased  separation  of  the  two  curves  near  the  signal  maxi¬ 
mum.  The  fine  structure  caused  by  Josephson  oscillations  is 
also  washed  out.  The  critical  current  density  of  both  junc¬ 
tions  is  4200  A/cm2,  giving  them  an  estimated  RNC  time  of 


FtO.  2.  Photograph  of  the  superconducting  integrated  circuit. 


TIME  IpS) 

FIG.  4.  Recorded  transition  duration  of  the  signal  current  /  for  a  junction 
with  critical  current  density  of  4200  A/cm1. 

2  ps  (R  v  is  the  junction  normal  state  resistance).*  The  effect 
of  jitter  is  roughly  equivalent  to  passing  the  signal  through  a 
low-pass  filter  of  rise  time  equal  to  the  jitter  width. 

Figure  3  is  a  photograph  of  the  integrated  circuit  used  in 
this  experiment.  It  is  a  seven-level  circuit  which  uses  a  lead 
alloy  process  similar  to  that  described  in  Ref.  5.  The  junction 
areas  are  approximately  40  and  80  jim2  as  indicated.  The 
microstripline  showing  at  the  left  is  a  small  portion  of  the 
line  used  to  obtain  600  ps  of  delay  in  the  output  signal.  The 
parameters  used  in  the  simulation  of  Fig.  2  are  taken  from 
measurements  on  this  circuit. 

Figure  4  is  the  experimental  transition  duration  data 
recorded  for  this  circuit.  The  two  traces,  recorded  about  one 
minute  apart,  are  for  increasing  and  decreasing  bias.  They 
are  different  because  the  sampling  time  drift  decreases  the 
apparent  transition  duration  in  one  case  and  increases  it  in 
the  other.  Assuming  constant  drift,  the  correct  transition 
duration  (10%-90%)  will  be  the  average,  or  about  9  ps.  This 
result  represents  the  combination  of  the  signal  transition  du¬ 
ration  estimated  to  be  about  5  ps  and  the  jitter  limited  system 
response  of  about  7  ps. 

A  second  experiment  was  performed  in  which  the 
switching  transition  of  / 1  was  intentionally  rolled  off  by  an 
on-chip  RC  filter.  In  this  case  the  recorded  waveform  was 
the  expected  exponential  rise.  The  time  constant  (3.7  ns) 
matched  the  filter  to  the  accuracy  with  which  the  filter  com¬ 
ponents  were  known  (30%). 

It  is  important  to  note  that  the  very  fast  switching  speed 
reported  here  is  only  one  part  of  the  delay  which  occurs  in 
logic  circuits.  Gheewala*  has  used  a  chain  of  gates  to  mea¬ 
sure  logic  delays  of  as  little  as  13  ps. 

This  work  is  supported  by  the  Office  of  Naval  Research 
under  Contract  No.  N00014-79-F-0020. 
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